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Abstract 


Munition  development  has  always  been  driven  by  the  necessity  of  delivering  enough 
explosives  to  a  targeted  object  to  destroy  it.  Targets  that  are  protected  by  steel  reinforced 
concrete  housings  have  become  increasingly  more  difficult  to  destroy.  Improvements  must  be 
made  in  munitions  engineering  design  to  either  deliver  more  payload  to  the  target  or  to  make  the 
weapon  more  potent.  In  most  cases,  due  to  aircraft  weight  limitations,  the  delivery  of  more 
payload  is  not  an  option.  Therefore,  improving  the  destructive  power  of  a  weapon  of  a  given 
payload  requires  the  use  of  more  powerful  explosives.  However,  when  the  potency  of  an 
explosive  is  increased,  its  sensitivity  to  premature  detonation  also  increases. 

The  characteristics  of  the  metal  casing  containing  the  explosive  contribute  significantly  to 
the  weapon’s  detonation  sensitivity.  The  casing  experiences  significant  heating  during  weapon 
penetration.  This  dynamic  heating  can  cause  the  weapon  to  detonate  before  it  reaches  its  target 
location.  In  the  past,  computer  codes  used  to  model  detonating  weapons  have  not  heating  into 
account  in  their  performance  predictions.  Consequently,  the  theoretical  models  and  the  actual 
field  tests  are  not  in  agreement.  New  models,  that  include  temperature  information,  are  currently 
being  developed  which  are  based  on  work  done  in  the  area  of  computational  fluid  dynamics. 

In  this  research,  a  remotely  located,  high-speed,  infrared  (IR)  camera  is  used  to  obtain 
detailed  measurements  of  the  passive  radiation  from  an  object  in  an  energetic  environment.  This 
radiation  information  is  used  to  determine  both  the  emissivity  and  the  temperature  of  the  surface 
of  an  object.  However,  before  the  temperature  or  emissivity  was  determined,  the  functional  form 
of  the  emissivity  was  assumed  to  be  an  Mh  degree  polynomial  with  respect  to  wavelength 
dependence.  Using  a  commercial  Fourier  transform  spectrometer,  this  functional  form  was 
verified  to  accuracy's  on  the  order  of  0. 1  %,  to  be  a  polynomial  of  the  second  degree. 

With  the  advent  of  large  (greater  than  128  X  128  pixels),  high-speed,  IR  detector  arrays,  it 
has  now  become  possible  to  realize  IR  imaging  spectrometers  that  have  very  high  spatial 
resolution.  The  IR  spectrometer  system  developed  in  this  research  utilized  a  large  detector  array 
to  allow  multiple  spectral  images  to  be  formed  simultaneously  on  the  image  plane,  with  the 
information  from  each  spectral  image  being  concurrently  collected.  In  conjunction  with  the 
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correct  emissivity  model,  this  imaging  IR  spectrometer  can  determine  the  temperature  of  the 
surface  of  an  object  within  an  energetic  environment  to  within  +  5  degrees  Celsius.  These 
experimentally  verified  temperature/emissivity  maps  can  then  be  integrated  into  the  newly 
developed  computer  models.  This  additional  information  will  result  in  more  accurate  computer 
codes  for  modeling  the  energetic  environment.  In  turn,  this  will  allow  the  weapon  designer  to 
accurately  optimize  weapon  performance  with  respect  to  different  materials,  geometries  and 
kinetics. 
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Introduction 


In  munition  development  there  has  been  a  great  need  to  accurately  characterize  the 
materials  found  in  exploding  and  deflagrating  reactions  to  determine  when  a  weapon  may  fail. 
The  primary  area  of  concern  within  munition  development  is  penetrator  and  explosively  formed 
projectile  type  weapons.  With  explosively  formed  projectiles,  the  first  concern  is  to  understand 
how  the  projectile  material  (typically  a  copper  alloy)  responds  to  the  explosive  charge.  After 
formation,  the  concern  is  to  understand  how  the  projectile  will  propagate  toward  the  target. 
Finally,  understanding  how  the  projectile  interacts  with  the  target  is  critical  for  effective  weapon 
design.  The  operational  scenario  of  the  explosively  formed  projectile  is  shown  in  Fig.  1  - 1 . 


Figure  1-1:  Operational  scenario  for  explosively  formed  projectiles. 
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Determining  the  temperatures  for  each  of  the  three  stages  in  the  scenario  can  become  quite 
involved.  However,  the  factors  that  largely  determine  instrumentation  needs  are;  time  of  event, 
image  resolution,  and  the  number  of  spectral  images  needed  for  accurate  temperature 
reconstruction. 

For  penetrators,  the  primary  concern  is  understanding  how  the  casing  material  (typically 
an  iron  alloy)  responds  to  the  stresses  of  penetrating  a  hardened  target  as  shown  in  Fig.  1-2. 


Figure  1-2:  In  penetrating  a  hard  target,  the  weapon  experiences  stress  induced 
by  the  multi  layers  of  steel  and  concrete. 


During  hard-target  penetration,  excessive  heating  of  the  weapon  casing  can  cause  massive 
warping  and/or  premature  detonation  of  the  bomb.  Thus,  there  is  a  need  to  characterize  how  the 
explosives  inside  the  casing  are  affected  by  these  temperature  changes.  This  is  an  important 
aspect  of  understanding  if  the  bomb  will  hold  together  until  it  reaches  its  target  and  subsequently 
detonates  in  the  designed  manner.  Determination  of  the  dynamically  changing  temperature  of  a 
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weapon’s  casing  material  is  essential  to  characterizing  its  plastic  deformation  characteristics  and 
failure  points. 

Computer  code  simulations  have  been  developed  as  input  to  munition  design  and 
development.  In  recent  years,  the  munitions  community  has  recognized  that  the  inclusion  of 
temperature  effects  in  these  high-speed  deformation  models  is  becoming  more  critical.  The  code 
developers  have  shown  that  when  temperature  characterizations  are  included  in  their  models, 
significantly  more  accurate  predictions  are  made.  However,  for  the  accurate  temperature  analysis 
of  the  next  generation  of  weapons,  the  development  of  new  innovative  pyrometry 
instrumentation  is  required.  There  are  two  main  forces  that  drive  the  level  of  difficulty 
associated  with  making  temperature  measurements  in  an  energetic  environment.  Measuring 
temperature  in  these  environments  is  difficult  due  to  high  rates  of  object  motion  (on  the  order  of 
300  meters  per  second  or  Mach  1.0)  and  the  unknown  emissivity.  In  addition,  the  temperature 
may  also  be  changing  at  rates  that  are  on  the  order  of  hundreds  of  degrees  per  millisecond.  To 
develop  the  proper  hardware  that  can  optically  measure  temperature  it  is  important  to  break  down 
the  overall  requirement  into  its  fundamental  instrumentation  needs.  The  temporal  and  spatial 
resolution  factors  of  the  instrument  design  are  somewhat  straightforward  and  easy  to  determine 
from  basic  radiometry.  Conversely,  before  a  pyrometer  is  designed,  the  number  of  necessary 
spectral  images  needed  for  accurate  temperature  prediction  must  be  determined.  Here  the 
functional  form  of  the  target’s  IR  spectrum,  in  terms  of  its  emissivity  and  temperature,  will 
determine  the  necessary  number  of  wavelengths  or  spectral  bandwidths. 

IR  radiative  characteristics  can  be  separated  into  an  emissivity  contribution  and  a 
temperature  contribution.  A  material’s  spectral  emissivity  is  defined  as  the  wavelength 
dependent  coefficient  that  characterizes  the  optical  radiation  of  a  surface.  To  find  the  emissivity 
of  a  material  measurement  of  the  material’s  radiant  power  at  each  discrete  wavelength  and  a 
blackbody’s  radiant  power  at  the  same  discrete  wavelength  are  obtained.  Temperature,  area,  and 
solid  angle  of  emission  are  kept  constant  during  these  measurements.  The  emissivity  for  each 
wavelength  is  then  found  by  taking  the  ratio  of  these  two  values.1  The  functional  dependence  of 
the  resulting  emissivity  curve  can  often  be  accurately  modeled  as  an  A/^-degree  polynomial.  The 
degree  of  this  polynomial,  and  the  functional  form  of  the  temperature  contribution  (typically 
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following  the  Planck’s  radiation  law),  directly  dictate  the  number  of  wavelengths  required  to 
accurately  determine  temperature. 

This  research  discusses  the  design,  development  and  test  of  a  four-color,  IR,  high-speed, 
imaging  spectrometer.  This  spectrometer,  in  conjunction  with  accurate  emissivity  models  and 
temperature  radiation  model,  is  used  to  produce  an  optical  pyrometer  system  that  solves  some  of 
the  temperature  measurement  needs  of  the  munition  development  community.  Specifically,  the 
characterization  of  dynamic  temperature  changes  of  casing  materials  is  investigated. 

The  goal  of  this  research  is  to  develop  an  accurate  theoretical  model  for  emissivity  to  be 
used  with  the  four-color,  IR,  high-speed,  imaging  spectrometer  to  determine  temperature.  This 
optical  pyrometer  significantly  reduces  the  instrumentation  inaccuracies  found  in  previous 
systems.  This  research  will  first  discuss  a  brief  history  of  pyrometry,  laying  the  foundation  for 
the  importance  of  this  work.  This  will  lead  to  a  discussion  of  a  temperature  measurement 
solution  to  the  problem  of  material  deformation.  In  addition  to  highly  accurate,  high-speed 
instrumentation,  the  accurate  measurement  of  temperature  requires  the  determination  of  an 
accurate  emissivity  model.  Thus,  a  theoretical  model  for  emissivity  of  metals  will  be  presented. 
The  accuracy  of  this  model  will  be  verified  using  a  commercial  Fourier  transform  spectrometer. 
The  details  of  the  design  and  test  of  a  high-speed,  IR,  imaging  spectrometer  will  follow. 
Combining  the  emissivity  model  with  this  spectrometer  produces  an  optical  pyrometer  system 
that  is  shown  to  be  capable  of  measuring  temperature  to  within  accuracies  of  ±  5  C°. 
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1.  Background 


When  small  amounts  of  force  are  applied  to  a  material,  it  undergoes  elastic  deformation, 
but  will  return  to  its  original  condition  after  the  force  is  released.  Applying  more  and  more  force 
to  the  material  will  eventually  cause  it  to  reach  its  “yield  point,"  where  it  will  first  begin  to 
plastically  deform  and  eventually  fracture.  Beyond  this  yield  point,  irreversible  structural 
deformation  change  occurs.  During  plastic  deformation,  changes  occur  to  the  material  resulting 
in  crystal  realignment  and  the  generation  of  heat  energy.  Most  of  the  deformation  related  energy 
is  converted  into  heat.  This  heat  energy  manifests  itself  as  a  change  in  the  material  s 
temperature.  Using  a  correctly  determined  emissivity  model  and  temperature  radiation  model  for 
the  casing  material,  a  spectrometer  can  be  used  to  determine  the  temperature  within  an  energetic 
environment.  This  chapter  will  discuss  how  a  heated  object  generates  an  IR  spectrum  and  the 
manner  in  which  the  temperature  of  this  object  is  determined. 


1.1.  Radiation  Law  and  Emissivity  Model 

A  certain  level  of  temperature  dependent  molecular  vibration  is  inherent  in  all  objects  that 
are  at  a  temperature  above  absolute  zero.  For  materials  that  are  both  in  thermal  equilibrium  and 
are  assumed  to  have  a  continuum  of  energy  states,  the  statistics  of  these  energy  modes  follow  a 
Maxwell-Boltzman  distribution.2  Although  the  actual  radiation  of  optical  energy  takes  on  this 
statistical  functional  form,  Planck  showed  that  energy  states  do  not  occur  as  a  continuum,  but 
rather,  they  occur  in  discrete,  quantized  levels.  Using  this  knowledge,  Planck  derived  the  ideal 
radiant  exitance  (where  emissivity  is  equal  to  1  at  all  wavelengths).  The  foim  of  this  ideal 
radiation  law  for  materials  in  thermal  equilibrium  is  known  as  the  Planck  distribution  law.3  This 
distribution  dictates  the  amount  of  energy  a  given  material  will  radiate  at  a  given  wavelength.  If 
the  material  is  a  perfect  radiator,  freely  radiating  and  absorbing  energy  to  and  from  free  space 
without  any  losses,  it  is  referred  to  as  a  blackbody.  The  actual  Planck  distribution  is  modified  by 
the  material’s  total  absorption.  The  total  absorption  of  a  material  combines  the  effects  due  to  a 
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material  s  inherent  bulk  absorption  with  the  losses  due  to  the  manner  in  which  energy  couples 
through  the  material’s  surface.  In  addition,  if  a  material  allows  no  light  transmission  and  it  is  in 
thermal  equilibrium,  then  the  material’s  total  absorption  will  equal  its  emission  (emissivity). 

Incorrect  assumptions  cause  errors  in  optical  pyrometers.  Often  the  largest  contributor  to 
this  error  is  uncertainty  in  the  functional  form  of  the  emissivity.  Emissivity  can  be  characterized 
with  a  Mth  degree  polynomial  in  wavelength  with 

=  aM+\^  JraM^  +  + . +  fli?  (1-1) 

where  the  a’s  are  unknown  constants  that  describe  emissivity  behavior.4  If,  for  example,  a  given 
material  has  an  emissivity  that  does  not  change  with  wavelength  (a  graybody),  we  would  only 
need  to  determine  a, ,  because  all  the  other  a’s  would  be  zero.  Consequently,  only  two  spectral 
measurements  would  be  necessary  to  form  two  equations  for  the  two  unknowns:  the  constant,  al , 

and  the  temperature,  T.  If  a  Mlh  -degree  polynomial  provides  an  accurate  fit  to  a  given  material’s 
emissivity  profile,  there  will  be  M  coefficients,  aM+v..a2 ,  of  wavelength  that  will  be  unknown. 
With  the  unknown  constant  emissivity  coefficient,  a, ,  and  the  unknown  temperature,  T.  a 
minimum  of  M+2  spectral  measurements  are  required. 

However,  prior  to  this  research,  instrumentation  system  errors  have  been  shown  to 
increase  strongly  as  the  number  of  spectral  measurements,  M,  increased.  This  has  made  the 
practice  of  taking  more  than  two  or  three  spectral  measurements  ineffective.5  These 
instrumentation  errors  are  primarily  attributable  to  using  separate  detector  arrays,  each  with  its 
own  support  electronics  and  optics,  to  collect  the  individual  spectral  images.  In  addition,  as 
higher  degree  polynomial  models  are  used,  the  error  in  calculating  the  M'h  -coefficient  is 
weighted  by  a  successively  larger  and  larger  wavelength  term  (i.e. 
A  =  5 fMn:  Pc'  =  125,  A4  =  625,  Pi  =  3125)  causing  a  wavelength  to  the  M‘h  error  in  the  calculated 

emissivity. 

In  the  past,  there  have  been  many  different  methods  used  in  attempting  to  measure  the 
spectral  characteristics  of  emissivity  and  temperature  for  a  variety  of  materials.  Depending  on 
the  assumptions  made,  and  a  priori  knowledge  of  the  material’s  characteristics,  various 


temperature  accuracies  have  been  achieved.  Historically  the  most  widely  used  algorithms  were 
initially  developed  for  the  use  in  metal  processing  and  later  were  applied  to  the  problem  of  casing 
material  deformation.  However,  these  algorithms  were  developed  using  various  assumptions  for 
emissivity  that  are  invalid  for  deforming  casing  materials.  These  assumptions  were  made  to  limit 
the  number  of  necessary  spectral  images  to  two.  This  was  done  for  the  same  reasons  stated 
earlier;  i.e.  prior  to  this  research,  instrumentation  errors  made  additional  spectral  measurements 
ineffective.  The  theoretical  results  suggest  that  if  the  instrumentation  was  redesigned  so  that 
system  errors  affected  each  of  the  spectral  images  equally,  more  spectral  measurements  could  be 
utilized  and  thus  a  highly  accurate  tool  would  be  produced. 

If  a  single  detector  array  can  be  used  to  collect  all  of  the  spectral  information  the  common 
set  of  optics  and  support  electronics  will  allow  all  spectral  images  to  vary  with  noise  in  the  same 
manner.  Consequently,  most  of  this  instrumentation  error  can  be  eliminated.  While  computer 
rounding  errors  can  also  cause  a  similar  problem  of  weighting  an  erroneous  coefficient  with  a 
larger  and  larger  value,  if  the  model  used  for  the  emissivity  is  less  than  a  sixth  degree 
polynomial,  the  round-off  error  will  not  be  significant.  This  research  will  show  that  temperature 
emissivity  separation  can  be  performed  accurately  with  a  second  degree  polynomial  model  for 
the  emissivity. 


1.2.  Separating  Temperature  and  Emissivity  from  Radiant  Power 
Measurements 

The  main  problem  behind  utilizing  thermal  radiation  in  making  IR  pyrometer 
measurements,6,7  is  in  attaining  accurate  knowledge  of:  the  spectral  emissivity  of  the  target,  the 
effects  of  atmospheric  absorption,  the  transmittance  of  the  optics,  and  the  spectral  response  of  the 
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detector.  Spurious  signals,  coming  from  other  emitting  bodies  that  lie  out  of  the  field  of  view  of 
the  detector,  make  the  measurement  even  more  difficult.  The  path  radiance  term  is  generally 
small  and  can  thus  be  ignored.  Including  all  these  terms,  the  basic  radiometric  equation  that 
maps  power  to  a  single  detector  element  is: 


<S>  =  LAQ, 


(1-2) 


where  L  is  the  radiance  of  the  object  for  a  specific  temperature,  A  is  the  area  of  the  object,  and  Q 
is  the  solid  angle  subtended  by  the  optics.  Taking  into  account  the  transmission  of  the  optics  and 
the  atmosphere  the  power  on  the  detector  equation  can  be  rewritten  in  differential  form  as: 

(1-3) 

r„(; l)  T'{X)Aad9dX 


dO(A,T>T„0)  = 


A  eXT  -  1 


s(X,T,0)- 


Sl  ) 

A5 

r-) 

(l  -  £(A.,T,0))  e,(Z,r„0) 


where: 

c,=  1.28x  104Wp4/cm2 
c2=  1.44  x  104n  K 

T—  Temperature  of  target 

Ts  =  Equivalent  temperature  of  the  background 

X  =  Wavelength  in  microns 
s(X,  T,  9)  =  Emissivity  of  target 

£s(A,  Ts,0)=  Emissivity  of  background 
ra  =  Transmission  of  atmosphere 
t0  Transmission  of  optics 

d9- Differential  solid  angle 
dA  =  Differential  wavelength 
A0  =  Area  of  the  object 
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The  second  term  inside  the  brackets  is  due  to  radiation  coming  from  objects  out  of  the  field-of- 
view  of  the  detector  but  reflecting  off  the  object  of  interest.  When  care  is  taken  in  setting  up  an 
experiment  and  the  object  of  interest  is  hot  compared  to  other  objects  in  the  surrounding  area,  the 
reflected  radiation  will  be  very  small  and  can  be  left  out  of  the  expression.8  Further  simplifying 
equation  1-3  can  be  done  when  emissivity  variations  across  the  solid  angle  are  negligible.9  The 
emissivity  will  vary  as  one  minus  the  reflectivity  and  it  is  known  that  the  reflectivity  for  metals 
will  see  very  little  variation  over  a  wide  angle,  0.10  The  emissivity  dependence  with  solid  angle 
(0)  can,  therefore,  be  dropped  and  equation  1-3  can  be  rewritten  as: 


d<5>(X,T)  = 


A5 


fL 

,xr 


£(A,T)ra(X)ro(A)Aod0dA . 


(1-4) 


After  calibration,  the  transmission  and  detector  response  terms  of  equation  1-3  arc 
known  and  the  basic  radiometric  equation  simplifies  to  two  unknowns.  The  problem  of 
determining  the  true  temperature  from  passive  IR  radiation  is  now  coupled  only  to  the 
uncertainty  in  predicting  emissivity.  It  is  the  attempt  of  this  section  to  give  some  understanding 
of  the  models  for  emissivity  that  were  developed  in  the  past  and  substantiate  their  usage.  The 
discussion  will  start  with  a  simple  model  that  uses  only  one  color  measurement  and  the  section 
will  conclude  with  a  new  model  developed  for  the  work  of  this  dissertation. 


1.2.1.  Single-Color  Pyrometry 


The  primary  assumption  behind  single-color  pyrometers  is  that  the  emissivity  for  the 
object  must  be  known  accurately  within  the  bandwidth  of  the  sensor.  After  the  emissivity  is 
known,  the  signal  at  a  single  detector  for  an  object  that  fills  the  instantaneous  field  of  view,  can 
be  represented  by. 
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era(X)To(A)A0QdX. 


(1-5) 
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The  limits  of  the  integral  have  been  set  to  take  into  account  all  spectral  energy.  However,  in 
reality  the  transmission  term  will  allow  only  a  certain  number  of  spectral  components  to 
contribute  to  the  power  at  a  given  detector.  The  temperature  can  be  solved  for  after  the  power  at 
a  detector  is  known.  First,  Equation  1-5  needs  some  simplification.  At  high  temperatures  (above 
100C°)  the  -1  of  the  denominator  is  small  compared  to  the  eXT  term  and  can  be  dropped.  The 
emissivity's  dependence  on  temperature  and  wavelength  has  also  been  removed  because  this 
technique  depends  on  an  accurate  model  of  emissivity  at  the  temperature  of  the  object  and  the 
wavelength  of  the  pyrometer.  Equation  1-5  can  then  be  solved  explicitly  for  T.  Depending  on 
the  accuracy  of  measuring  the  other  terms  in  equation  1-4,  reasonably  accurate  solutions  are 
found  for  temperature.  Since  temperature  variations  affect  the  incident  power  as  an  exponential 
and  the  other  terms  have  more  linear  effects, 11  maximum  sensitivity  is  obtained  at  the  steepest 
slope  of  the  blackbody  distribution  curve  that  is  to  the  right  of  its  peak.  Additional  sensitivity  is 
obtained  by  separating  the  IR  energy  into  two  spectral  regions  and  this  is  the  topic  to  the  next 
section. 


1.2.2.  Dual-Color  Pyrometry 

Dual-color  pyrometers  are  used  when  emissivity  is  not  known,  but  a  known  relationship 
exists  between  the  emissivities  at  any  two  wavelengths  (i.e.  ex ,  =  f[sx  2)  or 

ln(f/i,i)~  /(ln(^A,2))  )•  Surface  characteristics  of  emissivity  vary  with  many  factors  (oxidation 

layer  thickness,  surface  roughness,  etc.)  and  two  color  pyrometers  combine  two  spectral  radiance 
measurements  to  arrive  at  a  more  accurate  temperature  estimation.  The  most  commonly  used 
two-wavelength  algorithm  today  is  the  ratio  temperature  (RAT)  method.  The  derivation  of  the 
RAT  algorithm  begins  with  solving  for  the  emissivity  at  two  points,  i.e., 
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(1-6) 
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Where  the  M’s  are  the  measured  values  for  radiant  exitance  and  the  integral  has  been  dropped  for 
simplicity.  Assuming  the  ratio  of  the  two  emissivities  is  a  known  constant, 
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Taking  the  log  of  each  side  and  rearranging  terms,  equation  1-9  becomes 
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The  RAT  temperature  equation  is,  finally, 
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(1-13) 
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and  Tn  and  Tn  are  the  brightness  temperatures  that  would  be  inferred  if  the  emissivity  for  the 

target  was  1 .0.  In  practice,  there  are  some  errors  in  measurement  and  compensation  must  be 
made  for  calculated  temperature.  This  relationship  is 

1  A  B  „ 

'  Al  a2 

where  A  and  B  are  dependent  on  the  specific  algorithm  used12  and  C  represents  the  error  term. 
For  the  RAT  algorithm  A=Ar/X\  and  B=-Ar/X2-  This  technique  works  well  when  the  emissivity 
ratio  is  constant.  When  this  is  not  the  case,  other  algorithms  are  needed  for  accurate  results. 

There  are  two  classifications  for  dual-wavelength  algorithms.  The  first  type  can  be 
referred  to  as  the  exact  method.  These  algorithms  require  knowledge  of  how  emissivity  at  one 
wavelength  varies  with  the  emissivity  at  the  second  wavelength.  The  RAT  algorithm  is  an 
example  of  this  technique.  A  set  of  two  simultaneous  equations,  one  for  each  measurement  at  a 
given  wavelength,  can  be  solved  to  produce  an  exact  answer  if  the  emissivity  relationship,  such 
as  in  equation  1-13,  is  known  (given  no  measurement  error).  More  generally 
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(1-16) 
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dMXi  =  -(  y — ^^2’  (M7) 

can  be  solved  if  the  exl=f(&x  2)  relationship  is  known  (i.e.,  two  equations  and  two 

unknowns).  If  the  relationship  between  the  two  emissivities  is  known  but  fluctuates  over  time 
then  a  second  technique  can  be  employed.  A  least-squares  method  implemented  using  multiple 
measurements  over  time  to  reduce  the  effects  of  noise.  The  following  table  1-1  gives  a  listing  of 
the  commonly  used  dual-wavelength  algorithms  and  the  models/assumptions  used  for  the 
emissivity.  The  models/assumptions  have  been  made  with  a  great  deal  ot  knowledge  ol  the 
surface  roughness. 

In  the  measurement  of  the  temperature  of  casing  material  the  surface  roughness  of  the 
material  is  changing  due  to  deformation,  changing  the  relationship  between  the  two  emissivities. 
Consequently,  dual  color  algorithms  do  not  predict  temperature  accurately.  A  new  model  needs 
to  be  developed  for  emissivity  that  is  independent  of  surface  roughness,  and  only  needs 
wavelength  information  from  an  arbitrary  metal  surface. 
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1.2.3.  Multi-Color  Pyrometry 


When  the  models  for  emissivity  used  to  realize  a  two-color  pyrometer  are  not  accurate,  a 
new  model  needs  to  be  generated  and  this  will  require  additional  colors  (spectral  regions)  to  be 
collected  by  the  IR  instrumentation.  Knowing  the  response  of  each  of  these  spectral  regions,  a 
system  of  equations  can  be  generated.  Solving  this  system  of  equations  will  separate  the  coupled 
response  of  temperature  and  emissivity.  A  system  of  equations  is  constructed  from  the  energy 
residing  in  each  of  the  spectral  regions  of  the  pyrometer.  These  equations  are, 

oo 

<£>.  =  J — -7-“ — z(X)  ra{/ l)  T0[X)AfPdA,  for  i=l  to  n,  (1-18) 

0  A5L^-1 

where  n  is  the  number  of  spectral  regions.  The  area  of  the  object  and  solid  angle  is  known 
quantities  after  system  calibration,  and  can  be  represented  as  a  single  function,  i?,(/ 1).  In 

addition,  after  calibration,  the  transmission  terms  of  equation  1-17  can  also  be  incorporated  into 
the  this  same  function,  further  simplifying  the  expression.  The  power  on  the  detector  Equation 
(1-17)  can  be  rewritten  as 


O,  =  j — y~ - T-e{^Xi^Ri{X)dX ,  fori=Tton.  (1-19) 

0  A5  je^-lj 

The  first  equation  (i=l)  has  two  unknowns,  temperature  and  emissivity.  Each  additional 
equation  added  for  each  of  the  spectral  regions,  i=2,3,...,  introduces  one  more  unknown 
emissivity  with  temperature  remaining  the  same  in  all  equations.  Consequently  the  system  of 
equations  is  always  under  determined  (more  unknowns  than  equations)  and  some  assumptions 
need  to  be  made  about  the  behavior  of  emissivity.  The  typical  assumption  is  that  the  emissivity 
follows  a  Mlh  degree  polynomial  with  wavelength, 
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s(X)  ~aM+iAM  +aMXM  1  +aM_iAM  2  + . +ax. 


(1-20) 


There  are  M+l  unknown  coefficients  of  X  and  with  the  temperature  undefined  leaves  M+2 
unknown.  Consequently,  information  regarding  the  amount  of  optical  energy  in  M+2  spectral 
regions  is  needed  to  form  M+2  equations.  There  are  a  number  of  mathematical  techniques  that 
can  be  to  used  solve  a  system  of  equations  where  the  number  of  equations  is  equal  to  or  greater 
than  the  number  of  unknowns.  Several  techniques  will  be  compared  in  chapter  2,  Emissivity 
Model  Development,  Test  and  Analysis. 


1.3.  Prior  Pyrometer  Instrumentation 

Up  to  this  point,  the  discussion  has  been  limited  to  the  models/algorithms  used  to 
separate  temperature  effects  from  emissivity.  Now  the  instrumentation  used  in  the  past  to  collect 
the  spectral  data  will  be  detailed.  These  types  of  instrumentation,  in  addition  to  being  inaccurate, 
are  unable  to  acquire  the  required  spectral  information  quickly  enough  to  meet  the  needs  of 
energetic  materials  pyrometry.  In  the  past  several  years  there  have  been  many  attempts  to 
develop  fast  imaging  spectrometers.14,  15,  16  While  each  technique  has  had  its  strengths,  none 
showed  promise  to  meet  the  requirements  of  munition  development.  Consequently,  a  large  part 
of  this  dissertation  will  discuss  the  development  of  a  new  imaging  spectrometer  with  no  moving 
parts. 

Spectrometers  requiring  moving  parts  range  from  a  rotating  filter  wheel  to  the 
displacement  of  a  mirror  in  an  interferometer.  Most  of  these  systems  are  not  fast  enough  to  be 
used  in  actual  munition  testing,  but  can  be  used  to  validate  high-speed  instrumentation,  as  well  as 
a  variety  of  lower  speed  applications.  It  is  important  to  understand  these  systems  in  order  to 
substantiate  the  need  for  alternative  instrumentation.  There  are  three  techniques  used  in  moving 
part  spectrometers:  (1)  Moving  a  series  of  filters  in  front  of  a  fixed  detector  (moving  filter 
wheel);  (2)  Moving  the  image  over  a  series  of  detectors  each  with  its  own  filter;  and  (3) 
Changing  the  optical  path  difference  of  an  interferometer.  These  three  techniques  will  be 
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discussed  to  justify  why  a  new  technology  needed  to  be  developed  to  meet  the  requirements  in 
munition  testing. 

The  moving  filter  wheel  spectrometer  is  probably  the  most  basic  and  easiest  to  implement 
of  all  of  the  spectrometer  systems.17  The  idea  is  that  a  filter  can  be  placed  in  front  of  a  detector 
or  detector  array  for  which  the  image  of  some  object  is  formed.  After  data  is  collected  for  the 
one  filter,  a  second  filter  takes  its  place.  Eventually  .^enough  spectral  information  is  collected  to 
characterize  a  given  event.  The  strength  of  this  approach,  as  mentioned,  is  ease  of 
implementation.  In  addition,  the  filter  responses  can  be  characterized  providing  accurate  filter 
calibration  that  does  not  change  over  time.  This  is  a  good  system  to  use  for  general 
characterization  for  slowly  changing  events.  However,  it  suffers  from  the  mechanical  limitations 
of  rotating  a  filter  wheel  or  some  other  mechanical  device.  After  a  given  detector  integration,  the 
filter  needs  to  accelerate  and  then  stop  so  that  the  detector  can  accomplish  subsequent 
integrations.  This  can  take  a  substantial  amount  of  time,  and  when  the  spectral  radiance  of  the 
event  is  changing,  the  instantaneous  radiant  exitance  characteristics  can’t  be  evaluated. 

Some  of  this  limitation  can  be  reduced  by  replacing  the  discrete  filter  wheel  by  a  circular 
variable  filter  wheel  and  rotating  with  a  constant  angular  velocity18.  If  the  filter  is  in  collimated 
space,  then  the  overall  shift  of  the  pass  band  will  move  a  reasonable  amount  between  each 
subsequent  detector  readout.  The  limitation  changes  from  how  fast  the  filter  wheel  can  be 
rotated,  stopped,  rotated,  stopped,  to  a  limitation  in  the  speed  of  the  readout  for  the  detector. 
This  can  mean  substantial  increase  to  the  overall  speed  of  the  system  and,  consequently,  the 
speed  of  collection  of  multiple  spectral  regions.  If  the  spectral  separation  is  not  significant  then 
the  system  of  equations  will  be  highly  singular  mandating  the  collection  of  more  spectral  regions 
than  the  system  of  equations  indicates  (i.e.,  M+2  equations  M+2  unknowns).  Inherently,  in  these 
types  of  systems  the  uncertainty  in  the  variable  filter  wheel  position,  and  consequently,  the  filter 
pass  band,  creates  an  uncertainty  in  the  amount  of  collected  optical  energy.  In  high-speed  events 
this  uncertainty  becomes  unreasonable  for  the  accurate  reconstruction  of  the  temperature  of  the 
subject  event.  Even  with  a  moderately  accurate  emissivity  model,  the  reconstruction  of 
temperature  inevitably  becomes  a  limitation  of  the  instrumentation  collecting  the  spectral  energy. 
Hence,  the  filter  wheel  approach  is  not  the  instrumentation  of  choice  for  high  speed  events. 


26 


Instead  of  moving  the  filter  as  described  above,  the  image  may  be  moved  across  a  series 

of  detectors,  each  with  it  own  fixed  pass  band,  support  electronics  and  imaging  optics.  This 

system  is  typically  referred  to  as  a  framing  camera.19  The  approach  invariably  requires  the  use  of 

a  rotating  mirror  in  a  pupil  plane  of  the  optical  system.  Scanning  in  the  pupil  plane  avoids  the 

overlap  of  spectral  energy  from  detector  array  to  detector  array  and  singularity  of  the  system  of 

equations  can  be  circumvented.  This  type  of  device  is  inherently  faster  than  the  aforementioned 

technique  because  the  scanning  mechanism  is  lighter  and  can  be  rotated  faster.  While  this  system 

is  faster  than  the  rotating  filter  approach,  the  geometrical  setup  (separation  between  elements) 

limits  the  effective  F/#  to  a  minimum  of  about  20  or  higher.  With  the  knowledge  that  the 
radiometry  follows  a  dependence  proportional  to  ^-it  can  be  seen  that  the  efficiency  of  this  type 

of  system  is  at  least  two  orders  of  magnitude  less  than  a  typical  F/2  system  that  could  be 
implemented  in  a  single  color  pyrometer.  Consequently,  in  order  to  get  the  same  amount  of 
optical  energy  onto  the  detector  array  as  for  the  single  color  system,  the  integration  time  for  the 
framing  camera  must  increased  by  two  orders  of  magnitude.  This  makes  the  system  too  slow  for 
use  in  the  energetic  environment. 

Instead  of  moving  the  filter  or  the  image  across  a  series  of  detector  arrays,  a  Fabry-Perot 
can  be  scanned  through  wavelength.  A  Fabry-Perot  type  spectrometer  uses  two  reflective 
surfaces  to  form  a  resonator  cavity20.  Frequencies  of  light  that  are  tuned  to  that  cavity  are 
allowed  to  pass  more  freely  due  to  interference  losses  (at  each  of  the  two  surfaces)  from  the 
multiple  reflections  of  a  common  beam  of  light.  When  the  frequency  of  light  experiences  an 
even  number  of  wavelengths  in  the  cavity,  the  losses  will  be  minimized.  The  pass  band  of  the 
cavity  can  be  scanned  by  moving  one  of  the  mirrors.  A  detector  is  then  used  to  obtain  the  spectra 
for  the  light  entering  the  system.  However,  either  this  type  of  system  is  still  slowed  by  the 
readout  of  the  detector  or  the  scan  of  the  mirror  since  the  radiance  information  for  each  of  the 
spectral  regions  is  not  collected  concurrently. 

Instead  of  using  an  interferometer  to  pass  only  a  small  portion  of  the  wavefront,  one  may 
choose  to  delay  a  portion  of  the  optical  energy  and  interfere  it  with  the  original  signal21.  With  no 
optical  path  difference  between  the  two  arms  of  an  interferometer,  an  unaltered  signal  would  be 
acquired  at  the  detector.  As  one  of  the  mirrors  is  scanned,  the  detector  will  respond  in 
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accordance  to  the  partial  coherence  of  the  incoming  light.  After  scan  of  the  mirror,  the  Van 
Cittert-Zemike22  theorem  can  be  employed  to  reconstruct  spectral  content.  Simply  stated,  the 
detector  data  will  be  a  Fourier  transform  of  the  spectrum.  This  technique  is  referred  to  as  a 
Fourier  transform  spectrometer  (FTS).23  The  strength  of  this  method  is  that  it  has  high  spectral 
resolution.  However,  FTS  suffers  from  the  same  problem  as  the  Fabry-Perot  interferometer 
taking  a  significant  amount  of  tinre  to  move  the  scan  mirror  or  readout  the  detector. 

Since  the  energetic  environment  requires  the  acquisition  of  data  in  a  few  microseconds 
the  above  techniques  are  all  too  slow  to  be  used.  It  becomes  apparent  that  a  snapshot 
spectrometer  is  required.  There  are  two  ways  to  construct  an  imaging  snapshot  spectrometer. 
The  image  can  be  separated  into  multiple  images  on  a  focal  plane,  each  at  a  different  color,  or  the 
image  can  be  diffracted  over  the  focal  plane.  The  multiple  diffractive  orders  can  be  used  to 
reconstruct  the  image  at  some  later  point.24  While  both  of  these  techniques  have  shown  promise 
to  solve  the  energetic  materials  pyrometry  problem,  the  reconstruction  algorithm  needed  when 
the  image  is  diffracted  over  the  focal  plane  is  computationally  cumbersome,  and  takes 
tremendous  amounts  of  time  to  process  even  small  images  (10  pixels  by  10  pixels).  In  addition, 
since  the  diffracted  orders  will  not  fill  the  detector  array  evenly  the  focal  plane  between  the 
diffractive  orders  is  wasted.  Consequently,  the  method  of  choice  is  to  separate  spectral  energy 
into  multiple  images  on  a  focal  plane  concurrently.  After  looking  at  the  types  of  systems  that 
separate  images  onto  a  single  detector  array,  or  multiple  arrays,  it  quickly  became  obvious  that 
advantages  in  new  focal  plane  technology  could  be  exploited.  A  new  instrument  was  developed 
for  use  on  the  casing  material's  found  in  munition  testing. 
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2.  Emissivity  Model  development,  Test  and  Analysis 


This  chapter  will  discuses  the  development  of  a  model  for  emissivity  of  metals  that 
depends  only  on  the  spectral  form  of  IR  radiation.  This  model  will  be  validated  for  accuracy 
using  special  FTS  instrumentation.  In  addition,  the  algorithms  used  to  find  the  temperature  and 
emissivity  from  multi-spectral  measurements  will  be  discussed.  However,  before  the  model  for 
emissivity  is  developed,  the  thermal  equilibrium  assumption  needs  to  be  validated  to  show  that 
even  if  the  temperature  is  changing,  the  spectral  behavior  for  the  radiation  will  still  follow  a 
blackbody  distribution. 


2.1.  Validation  of  Thermal  Equilibrium  Assumption 

Light  emitted  from  solids  under  conditions  of  thermal  equilibrium  and  without  any  other 
external  sources,  is  known  as  thermal  light.25  The  Planck  blackbody  radiation  law  uses  thermal 
equilibrium  statistics  to  derive  the  ideal  wavelength  dependent  IR  spectrum.  However,  the 
surface  temperature  of  energetic  materials  may  be  highly  dynamic  and  the  assumption  of  thermal 
equilibrium  may  not  be  correct.  An  evaluation  is  needed  to  show  how  much  the  temperature  can 
change  during  a  given  time  interval  before  the  IR  radiation  spectrum  deviates  from  the  ideal 
blackbody  radiation  law.  This  evaluation  can  start  by  analyzing  a  material  that  has  two  energy 
states. 

Consider  a  cavity  of  unit  volume  whose  walls  contain  many  atoms  of  two  energy  states 
separated  by  an  energy  difference,  hv.  The  cavity  will  support  broadband  radiation26.  Now  let 
A,  and  N2  represent  the  number  of  atoms,  per  unit  volume,  occupying  energy  levels  one  and 

two.  Spontaneous  emission  creates  radiation  in  the  cavity,  assuming  that  at  least  one  atom  is  at 
energy  level  two.  Subsequently,  absorption  and  stimulated  emission  can  occur.  These  three 
processes  exist  in  equilibrium  when  steady  state  is  reached.  This  derivation  will  first  consider 
the  spontaneous  emission  rate  and  then  move  on  to  absorption  and  stimulated  emission. 
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The  probability  that  an  atom  in  the  upper  energy  state  undergoes  spontaneous  emission 
into  any  of  the  modes  within  the  time  interval  At  is  PspNt  =  A tjtsp  ,  where  tsp  is  the  spontaneous 

emission  time.  There  are  N2  such  atoms  and  the  average  number  of  photons  emitted  in  At 
seconds  is  therefore  N2At/tsp  .  This  is  also  the  number  of  atoms  that  leave  level  2  during  the 

time  interval.  Consequently,  the  rate  that  atoms  increase  in  level  2,  due  to  spontaneous  emission, 
is  given  by, 

dN  2  _  N2  r?_n 


The  solution  of  this  differential  equation  is  an  exponential  decaying  function  and  a  function  of 
initial  state  N02  at  time  zero; 


N2  =  N0  e  v  . 


(2-2) 


In  the  presence  of  radiation,  absorption  and  stimulated  emission  will  occur.  This  changes  the 
populations  of  iV,  and  N2 .  Since  there  are  N]  atoms  in  level  one  capable  of  absorbing,  the  rate 

of  increase  of  the  population  of  atoms  in  the  upper  level,  due  to  absorption,  is. 

dN 2  _  Njt2  (2-3) 

dt  tsp 

where  n  is  the  average  number  of  photons  per  radiation  mode.  In  a  similar  way,  the  stimulated 
emission  can  be  written  as, 


dN 2  _  N2n 
dt  tsp 


(2-4) 
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Equation  2-4  shows  that  the  absorption  and  spontaneous  emission  rates  are  proportional  to  the 
average  number  of  photons  in  each  mode.  Summing  the  three  processes  we  can  now  write  the 
rate  equation  for  the  overall  process  as. 


dN2  N2 

dt  fsp 


Nxn 


sp 


N2n 


(2-5) 


Notice  that  equation  2-5  does  not  include  transitions  to  other  than  the  two  levels.  However. 

adding  more  levels  is  a  trivial  extension  when  the  system  is  in  thermal  equilibrium.27  Thermal 

equilibrium  will  first  be  assumed  and  then  the  temperature  will  be  allowed  to  change  to 

determine  if  the  results  change.  In  steady  state,  the  rate  of  change  between  the  levels  will  be 
dN7 

zero,  — —  =  0 .  Consequently,  the  ratio  of  the  populations  of  the  two  states  is 
dt 


_  n_  1 
n 


(2-6) 


If  the  system  is  in  thermal  equilibrium,  the  populations  obey  a  Boltzman  distribution,  i.e., 


(2-7) 


and  the  radiation  mode  distribution  is 

_1 _ 

hv\ 

kBT)  _  j 


(2-8) 


The  relation  above  was  derived  under  the  assumption  of  thermal  equilibrium,  but  is  this 
assumption  a  good  one?  If  the  temperature  changes  from  Ij  to  T2  during  the  time  interval.  At, 

the  change  in  the  populations  will  be  related  as: 
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A N2  =  - AN , . 


(2-9) 


If  ten  atoms  move  to  level  two,  then,  as  a  consequence,  level  one  looses  ten  atoms.  The  change 
in  the  populations  can  also  be  written  as: 


A N2  = 


n 


1 


I+i 

n 


N, 


Jr. 


(2-10) 


n  is  much  smaller  than  one  and  the  one  in  the  denominator  can,  therefore,  be  dropped.  Equation 
2-10  can  be  rewritten  as: 


^  <2-") 

Equation  2-1 1  cannot  be  solved  without  knowledge  of  the  population  of  level  one;  however,  a  an 
approximation  for  A  AT  can  be  made.  In  equation  2-11  the  first  term  is  much  largei  than  the 

second  term  and  they  are  both  positive.  Dropping  the  second  term, 

(2-12) 

We  can  also  write  the  population  of  level  two  in  terms  of  level  one  as 

and  (2-13) 

A N2  *1  0N2. 

Consequently,  if  the  observation  time  is  one  microsecond  and  the  relaxation  time  is  10 
seconds  then  the  JAT  /dt  term  is  negligible,  i.e., 
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(2-14) 


dN2  AN2 
dt  At 


10  N2 
At 


«10  7  TV,. 


Equation  2-5  can  be  rewritten  in  terms  of  the  ratio  of  the  two  populations  as: 

TV2  n 

N[  t^XO1 +\  +  n' 


(2-15) 


^ sp  1 0  term  is  associated  with  the  changing  population  term  and  is  very  small  compared  to  1 . 
Therefore, 

dN2 

X”0’  <2-16> 

and  the  steady  state  assumption  holds  true  even  when  temperature  is  changing.  It  appears  that  the 
blackbody  radiation  law  will  provide  a  good  model  for  the  temperature  dependent  spectral 
radiation  until  the  integration  time  approaches  the  relaxation  time,  at  least  for  the  given  1 00° 
temperature  change. 


2.2.  Emissivity  Behavior  For  Ideal  Metal 

The  manner  by  which  that  an  ideal  solid  generates  an  IR  spectrum  has  been  discussed  and 
now  attention  is  diverted  to  the  characterization  of  a  metal's  absorption  spectrum  and,  its  ideal 
emissivity.  The  phenomenon  of  the  bulk  material’s  IR  absorption  and  emission  spectra  needs 
more  detail.  This  is  done  by  modeling  the  metal  as  an  electron  plasma  with  characteristic 
relaxation  time  and  effective  mass.  The  next  step  will  be  to  consider  how  the  surface  affects 
radiation  into  free  space.  The  surface  will  alter  the  total  amount  of  energy  at  a  given  wavelength 
that  is  allowed  to  escape  from  the  bulk  material.  Contributing  to  this  effect  is  Fresnel  reflection, 
surface  roughness  causing  diffraction,  surface  roughness  altering  refraction  and  oxidation  or 
other  surface  deposits.  The  free  electron  model  will  be  discussed  first. 
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2.2.1.  Free  Electron  Model 


When  in  equilibrium,  according  to  Kirchoffs  law,  the  absorption  of  a  material  is 
equivalent  to  the  emissivity.  Consequently,  if  absorption  characteristics  of  a  material  can  be 
determined,  the  emissivity  characteristics  will  be  known.  Since  much  of  the  munition 
application  is  considering  the  optical  properties  of  heated  metals,  the  derivation  will  start  with 
obtaining  the  absorption  relationship  for  a  simple  metal. 

The  earliest  attempts  to  predict  the  optical  properties  of  metals  were  made  by  Lorentz  and 
Drude.28  They  assumed  the  metal  to  contain  electrons  that  were  essentially  free  to  move  under 
the  influence  of  the  electric  field  induced  by  the  incident  electromagnetic  wave.  These  free 
electrons  are  the  valence  electrons  in  the  outer  shell  of  the  atoms  constituting  the  metal.  When 
an  electromagnetic  wave  is  incident  upon  its  surface,  an  oscillating  electric  field  of  angular 
frequency  to  is  induced  parallel  to  the  surface  in  the  metal.  The  free  electrons  will  oscillate  under 
the  influence  of  this  field  at  the  frequency  of  the  incident  wave.  Since  the  electrons  are  moved  by 
the  electric  field  they  must  absorb  some  of  its  energy.  The  absorption  of  a  material  can  be 
derived  from  its  dielectric  function.  In  turn,  the  dielectric  function  can  be  derived  from  a 
material’s  polarization  defined  by  the  density  of  dipole  moments.  To  derive  the  dipole  moment 
of  the  free  electrons  it  is  first  assumed  that  the  electrons  respond  to  an  electric  field,  Eft), 
polarized  in  the  x  direction.  The  equations  of  motion  for  each  electron  in  the  field,  using 
Newton’s  second  law,  are  written  as 


d2x(t )  me  dx(t ) 
c  dt 2  r  dt 


=  -eE(t) , 


(2-17) 


where,  x(t)  is  the  electron  displacement  from  equilibrium,  z  is  the  electron  relaxation  time  and 
me  is  the  effective  mass  of  an  electron.  Note  that  the  motion  of  the  nucleus  has  been  neglected  in 

Equation  2-17.  This  is  known  as  the  dipole  approximation29.  The  second  term  of  equation  2-17 
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represents  the  damping  due  to  collisions  with  ions.  Note  that,  because  the  electrons  in  a  metal 
are  free  to  conduct,  there  is  no  restoring  force  in  Equation  2-17. 

Assuming  a  single  wavelength  interacting  at  frequency  co ,  the  time  dependent  field  is 

E(t)  =  E0e-io* .  (2-18) 

r 

Accordingly,  the  position  of  the  electron  will  follow  the  field  oscillation  and  can  be  written  as: 

x(t)  =  x0e~'*.  (2-19) 

Substituting  Eq.  2-18  and  Eq.  2-19  into  Eq.  2-17  and  dividing  both  sides  by  e  10)1 ,  the  equation  of 
motion  becomes 


or 


2  1(0  TP 

-m„co  xn  -  m„  —  xn  =  - eE n . 


ex 


m„co 


cox  -i 
co2x2  + 1 


(2-20) 


(2-21) 


This  is  the  solution  to  the  differential  equation  of  motion.  The  polarization  density,  Pc ,  of  the 
conducting  electrons  is  the  dipole  moment,  -ex(t),  multiplied  by  the  electron  density,  ne.  and  is 
given  by 


Pc  =  ~neex{t)  =  - 


ne2x 


mea> 


- E(t ) 


cox  -i 
co2x2  +1 


(2-22) 


The  polarization  due  to  the  electrons  is  accompanied  by  the  polarization  induced  by  the  applied 
electric  field,  Pb,  which  can  be  thought  of  as  the  background  polarization.  Ph  is  due  to  the 

motion  of  bound  particles,  such  as  inner-shelf  electrons  about  their  equilibrium  position.  The 
total  polarization  is  the  summation  of  the  electronic  and  background  terms,  or, 
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p  =  pc+ph. 


The  optical  susceptibility  of  a  material  is  defined  as, 


PifQ) 

E{a>) 


(2-24) 


The  displacement  field,  D,  can  be  written  in  terms  of  the  electric  field  and  the  susceptibility, 


D-  E+  4kP  =  £[1  +  )]  . 


(2-25) 


In  addition,  the  relationship  between  the  displacement  electric  field  and  the  dielectric  function. 


D  =e(co)E, 


(2-26) 


can  be  used  to  obtain  the  relationship  between  susceptibility  and  dielectric  function. 


e  ( co )  =  1  +  4  7tx(co) . 


(2-27) 


Thus,  the  dielectric  function  can  be  represented  in  terms  of  polarization  density  and  electric  field 


:  (co)  =  1  +  4tt—  +  47t— 
EE 


,  .  Ph  A7me2x  coz-i 

=  1  +  4  7C-L - — 5-; - 

E  m„co  co  r  +1 


(2-28) 


=e\co)  +  i  e"(co). 
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The  imaginary  and  real  parts  will  be  used  to  compute  the  index  of  refraction  and  absorption 
respectively.  At  high  frequencies,  the  last  term  containing  Pc  is  negligible30  and  the  remaining 
terms  fold  into  the  high  frequency  term,  e®.  Equation  2-28  can  be  rewritten  as 


4 me2r  cor  -  i 
meco  (02t2  + 1 


(2-29) 


This  formula  is  known  as  the  Drude  oscillator  model  for  the  dielectric  function  and  can  also  be 


written  as 


e  {(o)  =e'm+  e\— 


4me2r2 


1 


m„ 


co2t2  + 1 


■  +  r 


4me2T 


1 


m„co  a>2r2  + 1 


(2-30) 


Now  the  coefficient  for  absorption  in  terms  of  dielectric  function  can  be  solved.  This 

O  1 

relationship  is  derived  in  many  texts  and  is 


where 


a(eo)  -  — — —  e”(<y) , 
cn(a>) 


(2-31) 


n(a>)  = 


2 


g'+V(g')2+(^')2 


(2-32) 


and  c  is  the  speed  of  light  in  a  vacuum.  Substituting  equations  2-30  and  2-32  into  equation  2-3 1 
and  arranging  terms  we  have 


J2co^'~ 


-j2ccAnm2T  1 


a(a>)- 


m„co  co 


V+l 


ATm^x2  1 


m,  co 


y+i  u 


(2-33) 


r  ,  4 1  V  ( 
K*  me  ©V+lJ 


4  me2r  1 

m  co  co2r  +1 
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Given  this  complex  equation  it  seems  rather  a  formidable  task  to  solve  for  the  absorption. 
There  are  several  material  contestants  that,  in  all  likelihood,  will  be  unknown.  There  are  five 
unknown  constants  in  the  system:  .  Since  emissivity  is  equal  to 

absorption,  an  optical  pyrometer  will  be  required  to  resolve  six  spectral  regions.  As  was 
discussed  in  the  last  chapter  a  system  of  equations  can  then  be  solved  to  determine  the  five 
unknown  constants  and  temperature. 

The  six  equations  made  by  the  six-color  pyrometer  are  nonlinear  and  require  a  nonlinear 
optimization  technique,  such  as  that  proposed  by  Levenberg  and  Marquardt.  These  types  of 
algorithms  tend  to  take  a  long  time  to  optimize  and  often  become  caught  in  local  minimums  and 
the  global  solution  is  not  always  found.  However,  some  simplifications  can  be  made  for  our 
application.  The  em  term  is  due  to  the  displacement  of  bound  electrons,  and  for  metals,  can  be 

assumed  to  be  unity.  Therefore, 


42coA7me2T 


a((Q )  = 


mm  o)2r2  + 1 


1- 


4  7me2T2 


m„ 


co2z2  +  1 


■  +  . 


1- 


A7me2T2 


m„ 


co2t2  +  \) 


+ 


A7me2T 


1 


V 


V  m co  co2t2 +  \J 


(2-34) 


4  7me2 

where  the - term  is  known  as  the  plasma  frequency.  Equation  2-34  takes  the  form 

m. 


a(a> )  = 


V2 


CT 


2  (A 


M+iM  +u. 


(2-35) 


where 


A  = 


4  met  1 

me  o)2r2  + 1 


(2-36) 


Conductivity  measurements  have  shown  that  the  relaxation  time  for  metals  is  on  the  order  of 

A 


10  14  seconds32.  Because  our  sensor  is  operating  between  3  and  5  pm,  t2co2  >  36  and 


0)T 


IS 
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small,  as  compared  to  the  (1-A)2  term,  and  can  be  left  out  of  the  expression.  Equation  2-35 
becomes 


(2-37) 


If  the  denominator  is  expanded  as  a  binomial  series,  equation  2-37  becomes 


,.1,1  11  31 

a(co)  A 

CT  4a  2  8  y— 


(2-38) 


For  metals,  the  magnitude  of  the  refractive  index  in  the  denominator  is  large.  As  a  consequence, 
the  higher-order  terms  are  small  and  can  be  left  out  of  Equation  2-38.  The  equation  simplifies  to 

—  A 

a(oo)  =  -rp=r  =  —  yfA  •  (2-39) 

\J-A\  CT 

If,  for  instance,  the  refractive  index  were  to  have  a  value  of  5,  A=26.  This  would  yield  a  value 
for  the  approximated  absorption  of  —5.1,  which  is  within  2%  of  the  unapproximated  value  of 

CT 

—  5.2.  Using  the  fact  that  o)2z2  is  large  compared  to  1,  equation  2-39  can  be  rewritten  as 
cr 


a(co) 


4  me"  1 
c2z1m„  (o 


(2-40) 


Therefore,  a  first-degree  polynomial  should  provide  a  good  model  for  our  calculations  in 
temperature  emissivity  decoupling.  There  are  no  constant  or  squared  terms  in  equation  2-40  and 
experimental  results  in  the  past  have  shown  emissivity  will  depend  on  more  than  just  the  linear 
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wavelength  term.33  This  additional  wavelength  dependency  is  shown  when  the  surface  is 
included  in  the  evaluation  of  how  a  material  generates  an  IR  spectrum. 


2.2.2.  Effects  of  Surface  on  Emissivity 

If  a  piece  of  polished  aluminum,  or  any  metal,  is  roughened  with  sand  paper  the 
reflectivity  decreases;  consequently,  emissivity  increases.  In  addition,  other  surface  features, 
such  as  oxidation  or  oily  residue  left  behind  during  processing  can  also  affect  spectral  emissions. 
In  this  section,  the  contribution  of  the  surface  effect  will  be  discussed,  and  the  emissivity  model 
will  be  modified  to  incorporate  these  considerations. 


2.2.2. l.Fresnel  Reflection  at  Surface  Boundary 

The  characteristic  energy  modes  that  are  available  for  radiation  are  known  from  the 
Planck  radiation  law  and  the  method  by  which  energy  is  reflected  back  into  the  bulk  material  by 
the  boundary  needs  to  be  defined.  The  three  fields  of  concern  are  the  incident  field,  En  the 
reflected  field,  Er ,  and  the  transmitted  field,  Et ,  represented  in  Figure  3.1. 
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Figure  2-1:  Surface  fields  of  a  metal  sample 


The  E  and  H  fields  for  the  system  are34: 
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Maxwell’s  third  equation  can  now  be  used  to  write 
dEz  dEy 
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Canceling  the  terms  that  are  zero,  Equation  2-42  becomes 


(2-41) 


(2-42) 


42 


(2-43) 
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H  =  —  E 

'oO«o  „  i„orl0  ' 


(2-44) 


For  the  reflection  component,  from  Equation  2-45,  there  is  a  sign  change  on  t,  consequently. 


c 


(2-45) 


Now  the  boundary  conditions  need  to  be  incorporated.  The  tangential  component  will  be 
continuous  at  the  boundary,  e.g., 


E,o  +  Er o  =  E, o  and  Hi0  +  Hr0  =  H,0 , 


(2-46) 


CL e  +  CLL  £  -  CLLLLL  £ 

C  C  C 


(2-47) 


{Eio-Er0)  =  ^{E,i0+Er0) 


(2-48) 


the  transmitted  wave  is  in  free  space,  therefore  s0  =  1  and  Eq.  2-52  becomes 


-  ^o(^  +  l)  > 


(2-49) 


Er  _  ye  - 1 
E,  a/g  + 1 


where 


==e’+z  e" . 


(2-50) 


(2-51) 
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However,  since  r2co 2  >  36 ,  the  thermal  energy  lies  in  the  high  frequency  regime.  Consequently, 
the  imaginary  part  of  the  dielectric  function  is  small  compared  to  the  real  part  and  can  be 
neglected35.  The  reflectivity  can  now  be  written  as 


where 


i_J_ 

r  e, ,  v?-i  4^ 
E,0  V5  +  1  1  +  _L’ 
4^ 


(2-52) 


CO 
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(2-53) 


Equation  2-57  can  also  be  written  in  terms  of  the  transmission  of  the  surface.  The  transmission 
defines  the  amount  of  energy  that  escapes  the  surface  and  radiates  into  free  space. 


1  + 


f  =  1- 


Jr 0  _ 


V?  M  M 
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1  + 


1  + 


1  +  - 


(2-54) 


The  co p  refers  to  the  plasma  frequency  introduced  earlier.  The  plasma  frequency  is  typical  in  the 

ultraviolet36;  consequently,  for  IR  radiation  the  real  part  of  the  dielectric  function  is  negative. 
Since  the  IR  frequencies  of  interest  are  well  below  the  plasma  frequency  the  binomial 
approximation  can  be  used  to  simplify  Equation  2-58. 


_  _2_f _ 1_)  _2 _ 2_ 
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(2-55) 


The  second  term  of  this  expression  can  be  dropped  because  it  is  very  small  compared  to  the  first. 
The  transmission  then  becomes 
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(2-56) 


where 


t(co)  = 


con 

— ! -  >  1  and 
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»  1  . 


(2-57) 


Transmission  is  the  square  of  transmissivity,  hence 


T{co)  =  t 2 


(2-58) 


Rewriting  in  terms  of  wavelength 


T{X)  = 


\6k2c2 

V2  ' 


(2-59) 


This  model  needs  to  be  taken  one  step  further  in  the  evaluation  of  other  surface  effects.  For  the 
model  to  be  accurate,  the  effects  of  scratches  and  grain  structure  of  the  material  need  to  be 
included. 


2.2.2.2.Surface  Topography 

The  published  literature  on  the  emissivity  of  metal  surfaces  largely  ignores  the  wavelength 
dependent  effects  of  diffraction  or  scatter.  Diffraction  is  caused,  in  part,  by  surface  topography 
errors  that  are  induced  by  the  manner  in  which  the  surface  was  fabricated  or  prepared.  Surface 
topography,  as  well  as  the  inherent  crystalline  grain  boundary  structure  of  the  metal  surface,  causes 
emitted  energy  to  diffract  through  the  surface  of  a  material.  For  many  years,  the  published 
literature  on  the  diffraction  theory  of  scatter  has  treated  surface  errors  by  representing  them  as  a 
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summation  of  sinusoidal  gratings  of  varying  spatial  frequencies.  Each  component  sinusoidal 
grating  causes  a  single  wavelength  to  diffract  into  an  angle  that  is  not  predicted  by  Snell  s  Law  or 
the  Fresnel  equations.  This  angular  dependence,  however,  can  be  predicted  by  using  the  grating 
equation: 


sin  9  -  sin  6i  =  —j~  (2-60)  ■ 

a 

When  using  an  optical  system  to  measure  emissivity,  the  solid  angle  subtended  by  the 
optical  system  dictates  which  diffracted  angles  can  be  collected.  Since  the  desired  lesult  is 
emissivity  as  a  function  of  wavelength,  and  not  as  function  of  angle,  it  is  more  appropiiate  to  use 
the  grating  equation  to  see  what  wavelengths  will  be  collected  for  a  fixed  range  of  scattering  angles. 
Thus,  the  angular  dependence  of  the  diffracted  energy  can  be  rewritten  with  a  wavelength 
dependency. 

It  has  been  shown  by  DeWitt37  that  when  a  metal  surface  is  systematically  roughened  with 
sandpaper,  the  spectral  characteristics  of  its  measured  emissivity  will  change.  This  emissivity 
change  is  a  direct  representation  of  the  change  in  the  nature  of  the  energy  that  escapes  from  the 
body  of  the  material  and  into  free-space.  One  possible  explanation  for  this  spectral  change  in 
emissivity  can  be  attributed  to  the  characteristics  of  the  surface,  where  optical  energy  is  coupling 
through  the  surface  by  diffraction.  The  results  of  Dewitt’s  wavelength  dependent  emissivity  can  be 
fit  by  a  multiple  order  polynomial: 

e(X)  =  AXm  +  BAm~'  +CAm~2+....  (2-61) 

With  a  detailed  analysis  of  the  features  present  on  the  surface  of  the  material,  the  polynomial 
coefficients  (A,  B,  C,  ...)  can  be  found. 

A  detailed  analysis  of  the  diffraction  related  emissivity  is  beyond  the  scope  of  this 
dissertation.  However,  the  general  approach  should  be  outlined.  First,  the  statistics  of  the  surface 
topography  must  be  determined.  The  topography  can  then  be  broken  down  into  its  fundamental 
sinusoidal  components.  Each  of  the  components  can  then  be  treated  separately  as  a  single 
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diffraction  grating,  and  the  effects  of  all  the  sinusoidal  gratings  can  then  be  summed  to  determine 
the  total  diffraction  related  emissivity.  This  effect  can  then  be  added  to  the  Fresnel  and  bulk- 
material  effects  to  determine  the  total  emissivity. 


2.2.2.3.Surface  Chemistry 

As  mentioned  earlier,  a  real  metal  surface  unavoidably  carries  some  film  of  one  type  or 
another  ,and  this  affects  how  the  IR  energy  escapes  the  bulk  material.  In  the  machining  of  an 
aluminum  part,  for  instance,  oils  and  other  deposits  are  left  behind.  Often,  the  deposits  are 
oxides  of  the  base  material,  and  the  characterization  is  difficult.  The  layer  is  a  mixture  of  several 
different  chemical  species  (metal  atoms,  oxygen  ions,  and  one  or  more  metal  ions),  making  the 
interface  of  the  base  metal  rough.  In  addition,  oxides  form  on  the  surface  in  accordance  with  the 
atmospheric  conditions,  surface  topography,  temperature,  and  grain  structure.  The  growth  of 
oxides  on  a  surface  is  not  fully  understood,  however;  the  general  scenario  is  a  roughening  of  the 
surface,  the  effect  of  which  was  presented  in  the  last  section.  The  assumption  made  here  is  that 
the  bound  states  associated  with  the  oxides  do  not  contribute  significantly  to  IR  radiation.  ’8  This 
is  a  good  assumption  as  long  as  the  bound  electrons  are  held  tightly,  as  in  high-energy 
interactions. 

The  coatings  of  oxide  films  as  well  as  other  material  films  can  have  a  significant  effect  on 
the  thermal  emissions  from  a  surface.  If  the  consideration  can  be  restricted  to  oxides,  the 
evaluation  becomes  simplified  to  a  changing  in  the  roughness  of  the  surface.  The  objects  used  in 
the  testing  of  the  above  theoretical  results  were  prepared  carefully  to  limit  the  surface  chemistry 
to  only  oxidation  formed  on  the  surface.  In  this  way  the  effects  of  emissivity  could  be  validated 
and  tied  directly  to  the  free  electrons,  Fresnel  reflections,  and  surface  topography  of  the  test 
object. 


47 


2.3.  The  Overall  Emissivity  Model 


It  has  been  shown  that  the  bulk  material  and  Fresnel  reflections  at  the  boundary  have  a 
wavelength  dependent  effect  on  the  emissivity.  However,  the  true  behavior  of  emissivity  can 
only  be  theoretically  determined  if  the  surface  features  causing  diffraction  are  known.  The  total 
emissivity  is  a  compilation  of  Equations  2-40,  2-59  and  2-60.  Multiplying  these  three  equations 
we  have. 


t v  3  16;t2c2 

e(X)  =  Vr  ~ - -jT  (AXm+BXm-'  +  CXm~2  +...), 

C  T  e  im  COpX 


(2-62) 


If  the  emissivity  model  for  diffraction  follows  a  second  order  polynomial  we  have 


s{X)  =  - 


X  1 60zV  , 
- *'+• 


C  T  e  2  TTC  CO„ 


X  16  Attc- 

2  3  T  X  +  i 

C  T  e  2  TIC  COp 


<*>’.  X  16  Brrc1 


C2T3t.  CO 


(2-63) 


which  is  in  the  form  of 


s(X)  =  DX-]  +EX  +  F .  (2-64) 

The  X'1  can  be  represented  as  a  polynomial  by  using  a  Taylor  expansion.  Expanding  around  X=4 
microns  we  have 

A"=r(i-4)^+(A-4)J^-(1-4)15k-' 

»  0.754 - 0.1 9X  +  0.0 1 6X2  -6.1jc10"5/13 


The  terms  larger  than  the  X2  term  are  small  and  can  be  ignored;  thus,  Eq.  2-63  the  model 
becomes, 


e(k)  =  D0.016X2  +  ( DE  -  £>0.1 9)T  +  £>£0.754 .  (2-66) 
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We  are  left  with  a  second  order  polynomial  as  the  fit  for  the  emissivity  profile. 

While  the  explicit  diffraction  effects  have  not  been  included  in  this  analysis,  the  general 
approach  has  been  discussed.  The  second  order  polynomial  in  equation  2-64  has.  in  the  past, 
provided  a  good  characterization  of  the  emissivity  profiles  for  metals39.  Consequently,  the  next 
section  will  determine  the  emissivity  model  using  real  data. 


2.4.  Verification  Of  Model 


To  validate  the  emissivity  model  proposed  in  this  dissertation,  an  experiment  was  set  up 
using  heated  plates  of  aluminum  and  copper.  The  samples  were  sanded  to  remove  any  oxidation, 
attached  to  an  extended  blackbody,  and  heated  to  various  temperatures  from  100  to  250  K.  A 
Fourier-transform  IR  spectrometer  manufactured  by  Bomem,  Inc.,  was  used  to  collect  spectral 
information,  and  a  system  of  equations  was  constructed.  As  shown  in  Equation  2-1,  the 
emissivity  for  a  material  can  be  expanded  to  as  many  orders  as  there  are  spectral  regions  in  our 
spectrometer.  However,  due  to  computational  rounding  errors,  the  calculations  of  coefficients 
beyond  the  seventh-order  are  erroneous.  Consequently,  the  experiment  gathered  only  seven 
spectral  values.  This  was  sufficient  for  the  validation,  because  the  series,  in  our  model,  should 
converge  after  the  third  term.  The  seven  coefficients  for  the  emissivity  model  were  calculated 
using  experimental  data  collected  for  aluminum  and  copper  samples.  The  system  of  equations 
was  set  up  as  follows: 


ADU^  J 


T  1  +<33/ 1  +#5/l  +#7) 


XT  - 

e  -1 


^(A)  dX  ,  for  i=l,2,..,7,  (2-67) 


where: 

ADU^=  The  number  of  Analog  Digital  Units  (ADU)  output  for  each  detector 
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T=  Temperature  of  target 
A  =  Wavelength  in  microns 

R  (2)=  Response  of  detector/filter  and  transmission  of  optics  in  each  of  the 

spectral  regions  (ADU/Watt) 
dA  =  Differential  wavelength 


The  R,(A)  term  defines  the  limits  of  the  integral.  Note  that  detector  gain,  offset,  solid  angle,  and 
atmospheric  transmission  terms  have  been  left  out  of  equation  2-67.  The  system  did  not  have  a 
direct  means  to  calibrate  out  these  effects  and,  consequently,  the  unknown  constant  terms  were 
lumped  into  the  unknown  constants  for  emissivity,  n,  through  a7 .  Rearranging  the  terms, 

equation  2-77  becomes 


The  system  of  equations  is  non-linear,  but,  when  the  temperature  is  known  (as  for  the 
validation  experiment)  Equation  2-68  becomes  linear.  The  -1  in  the  denominator  is  small 

compared  to  the  first  term,  ex,r ,  and  can  be  left  out  of  the  expression.  The  temperature 

c2 

dependent  constants,  ex,T ,  as  well  as  other  constants  of  the  system,  are  lumped  into  an  overall 
temperature  dependent  constant,  cn  ( T) .  The  equations  can  be  rewritten  as: 
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for  i=l,2,..,7. 


Each  of  the  integrals  can  be  evaluated.  Rewriting,  Eq.  2-66  becomes 

ADU \  =  a, 7,  (7)  +  a2/2  (7)  +  a3/3  (7) 

+  a4/4j  (7)  +  a5/5.  (7)  +  a6/6>  (7)  +  a7/7  (7),  1 

or  can  be  rewritten  in  linear  algebra  form  as 
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• 

ADUh 

4(7)  4(7)  4(7)  4(7)  4(7)  4(7)  4(7) 

ADU , 

4(4  4(4  4(4  4.(4  4(4  4(7)  4(7) 

«6 

ADUXi 

.4(4  4(4  4(7)  4(7)  4(7)  4(7)  4(7) 

a7 

The  solution  for  the  unknown  coefficients  a,,  a2  ,  a3...  is 
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(2-69) 


(2-70) 


(2-71) 


(2-72) 
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The  data  from  the  spectrometer  were  substituted  into  equation  2-72  and  the  solution 
attained.  The  terms  are  higher  than  second  order  were  effectively  zero.  Consequently,  these 
terms  are  omitted.  Since  there  were  no  calibrated  emissivity  plots  for  the  aluminum  and  copper 
test  items  only  the  general  trend  in  data  could  be  calculated.  The  general  setup  for  the  test  is 
shown  in  Figure  2-2. 


Figure  2-2:  Validation  experiment.  The  components  for  this  experiment  were 


delivered  on  a  contract  from  the  government  to  Bomem,  Inc. 

This  experiment  did  show  that  the  emissivity  for  the  two  metals  was  modeled  very  well 
by  a  second-degree  polynomial.  Figures  2-3  and  2-4  show  a  plots  of  the  fit  for  emissivity.  As 
higher  order  fits  were  accomplished,  no  additional  accuracy  was  gained. 
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Figure  2-4:  Plot  of  cooper  emissivity  profile.  The  solid  line  represents  second 
degree  polynomial. 


This  validation  experiment  gave  credence  to  the  approach  of  building  a  four-color  pyrometer. 
Before  moving  on  to  the  technical  discussion  of  this  new  instrumentation  the  algorithms  used  to 


reconstruct  the  validation  experiment  data  should  be  detailed.  These  algorithms  were  also  used 
in  the  four-color  pyrometer  system. 


2.5.  Algorithm  design 

As  mentioned  earlier,  the  system  of  equations  generated  from  the  multi-spectral 
measurements  is  nonlinear  due  to  the  terms  associated  with  the  Planck  distribution; 
consequently,  a  nonlinear  regression  would  be  required  to  attain  a  solution.  These  equations  can 
be  made  linear  by  holding  the  temperature  constant.  After  a  solution  is  attained,  an  error  function 
establishes  the  accuracy  of  the  predicted  temperature.  A  regression  through  temperature  will  find 
the  true  temperature  at  the  minimum  error  function.  In  this  section,  two  approaches  to  finding 
the  linear  solution  will  be  presented,  and  the  section  will  conclude  with  a  nonlinear  technique. 

In  the  last  section  it  was  shown  that  a  second  order  polynomial  provides  a  accurate  fit  lor 
the  emissivity  profile.  Consequently,  only  spectral  information  is  needed  about  four  regions  ot 
the  emissivity  profile.  Simulated  multi-spectral  data  was  used  to  check  each  algorithm  for 
performance  against  varying  emissivity  profiles  and  amounts  of  noise  in  the  system.  All  three 
techniques  were  used  to  reconstruct  the  temperature  using  the  new  imaging  pyrometer,  and  a 
comparison  is  presented  in  the  results  section. 


2.5.1.  Algorithm  Design  And  Simulation  (Gaussian  Elimination) 

The  first  linear  solution  technique  uses  Gaussian  elimination  as  the  basis  to  find  a 
solution  to  the  system  of  four  equations  constructed  from  four  spectral  measurements.  With  a 
second  degree  polynomial  as  the  model  for  emissivity  the  four  equations  setup  by  the  four  color 
pyrometer  are 
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ADU^  =  f"  t  *,(2)  dl  for  i=I,  2,  3  &  4 
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(2-73) 
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Rearranging  the  terms,  equation  2-73  becomes 
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If  the  temperature  is  held  constant,  then  any  three  of  the  four  equations  has  a  linear  solution.  The 
first  three  equations  can  be  written  as 


A1*  i?(/l)  (R  f1™ 

q(7)-^— +^j  q(7) 


4 

/•7mK 

x. 


a  * 


j  <XJ) 


7“ 


for  i=l,  2,  &  3 


(2-75) 


Each  of  the  integrals  can  be  evaluated, 


ADUXi  =  axIh  (T)  +  a2I2i  (T)  +  a,Iy  (T) 


for  i=l,  2,  &  3, 


(2-76) 


or  can  be  rewritten  in  linear  algebra  form  as 
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The  solution  for  the  unknown  coefficients  a,,  a2,  &  a3  is 
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There  many  techniques  for  inverting  a  matrix  and  their  usage  depends  on  the  specific 
application.  The  Gaussian  elimination  method  was  chosen  because  it  provides  the  greatest 
versatility  and  it  is  easily  implemented.  The  code  implementing  matrix  inversion  was  taken  from 
Numerical  Recipes40  and  is  documented  in  appendix  A. 

In  addition,  in  appendix  A  is  an  algorithm  that  was  used  to  analyze  thef.viability  of  the 
linear  solution  discussed  above.  The  first  step  is  to  simulate  real  data.  Figure' 3.5  shows  the 
emissivity  profile  that  was  generated  in  the  VARCAL  subroutine  by  fitting  a  second-degree 
polynomial  to  three  points,  chosen  arbitrarily.  The  emissivity  was  used  with  the  spectral  filter 
characteristic  to  simulate  the  energy  seen  at  each  of  the  four  detectors  in  the  DATUM  subroutine. 
This  generated  the  simulated  detector  responses  for  the  imaging  spectrometer.  The  filter 
characteristics,  also  plotted,  were  generated  in  the  FILTERS  subroutine.  A  Gaussian  filter  shape 
was  used  for  simplicity.  In  the  real  system  the  filter  shapes  were  replaced  with  the  actual  filter 
characteristics  and  a  solution  for  temperature  found  in  the  same  way  as  for  the  simulation. 

To  find  the  solution  for  temperature,  the  coefficients  generated  from  Gaussian  elimination 
are  used  in  the  fourth  equation.  If  there  is  an  error  in  this  fourth  equation,  the  assumed 
temperature  is  increased  or  decreased,  and  the  above  process  repeated.  In  addition,  a  random 
error  (RND(l)  *  ER)  was  added  to  simulate  noise.  The  FITDATUM  subroutine  finds  the 
coefficients  to  the  polynomial  data  fit,  where  temperature  is  constant,  and  uses  the  fourth 
equation  to  calculate  the  error  in  that  fit.  The  main  routine  loops  through  all  values  of 
temperature  and  the  error  is  output  to  the  screen  for  each  temperature.  The  solution  for 
temperature  is  found  when  the  error  in  the  fourth  equation  is  minimized. 
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performance 


This  algorithm  has  performed  quite  well  on  the  simulated  data.  With  no  noise  added  to 
the  data  (i.e.,  ER=0),  the  error  at  the  correct  temperature  is  essentially  zero  when  neglecting 
computational  rounding  errors  (eighth  decimal  place,  or  so).  As  error  is  added,  the  algorithm 
deviates  from  the  ideal.  However,  if  the  noises  are  within  the  Poisson  photon  noise  limited,  seen 
in  a  low  noise  detector  system  ,  then  the  system  will  be  able  to  attain  much  better  than  one- 
degree  accuracy. 


2.5.2.  Singular  Value  Decomposition 

In  the  last  subsection,  a  straight  forward  solution  to  the  system  of  equations  was 
presented,  where  a  simple  matrix  inversion  was  utilized.  If  this  matrix  was  either  singular  or 
numerically  close  to  singular,  inversion  using  Gaussian  elimination  would  fail  due  to 
computational  rounding  errors.  However,  Singular  Value  Decomposition  (SVD)  can  be 
successfully  used  to  invert  such  a  matrix.  In  addition,  SVD  is  a  powerful  method  for  solving 
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most  linear  least-squares  problems.  For  the  imaging  spectrometer  of  this  dissertation  the  nature 
of  singularity  will  only  be  experienced  when  the  images  have  significant  overlap  between  the 
adjacent  spectral  regions.  Overlap  can  be  essentially  eliminated  with  careful  design  and, 
consequently,  SVD  will  be  used  as  a  method  to  accomplish  a  least-squares  fit  of  the  data. 

SVD  methods  are  based  on  the  following  theorem  of  linear  algebra,  the  proof  of  which  is 
beyond  the  scope  of  this  dissertation:  Any  matrix  A  who’s  number  of  rows  M  is  greater  than  or 
equal  to  its  number  of  columns  N,  can  be  written  as  the  product  of  an  M  X  N  column-orthogonal 
matrix,  U,  an  N  X  N  diagonal  matrix,  W,  with  positive  or  zero  elements  (the  singular  values), 
and  the  transpose  of  the  N  X  N  orthogonal  matrix,  V. 
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The  matrices  U  and  V  are  each  orthogonal  in  the  sense  that  their  columns  are  orthonormal. 
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Since  V  is  square,  it  is  also  row-orthogonal,  V  ■  V1  =  1 . 

The  intent  of  this  method  is  to  solve  the  same  system  of  equations  as  was  solved  for  in  the 
last  section.  Again  these  equations  are  non-linear  but  can  be  made  linear  by  assuming  a 
temperature.  This  will  give  us  four  equations,  one  from  each  of  the  four  colors  of  the  imaging 
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spectrometer,  and  three  unknowns  for  the  coefficients  of  emissivity.  This  is  represented  in  linear 
algebra  form  in  equation  2-91  and  2-92. 
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Writing  in  the  form  of  the  decomposed  matrix,  we  have 
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After  the  SVD  calculates  the  coefficients  for  emissivity  they  can  be  substituted  back  into  the  four 
equations  and  an  error  calculated  for  each.  The  temperature  can  then  be  incremented  and,  as  for 
the  Gaussian  elimination  technique,  the  solution  for  temperature  is  reached  when  this  error  is 
minimized. 

The  implementation  of  an  SVD  algorithm  developed  in  Numerical  Recipes  provided  a 
solution  to  equation  2-93  and  is  also  listed  in  the  appendix  A.41  A  function  call  from  the  main 
program  to  the  SVD  subroutines  is  made  in  the  same  manner  as  was  done  for  the  Gaussian 
elimination  regression  techniques.  This  algorithm  uses  the  same  simulated  data  (emissivity 
profile,  filter  responses,  error,  etc.)  to  verify  the  accuracy  of  the  fit  to  the  data  in  the  presence  of 
noise.  While  this  application  of  multi-spectral  pyrometry  generates  a  matrix  that  is  non-singular, 
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the  SVD  does  provide  a  least-squares  fit  of  the  data  that  is  accurate,  and  nearly  as  accurate  as  the 
Gaussian  elimination  technique. 

2.5.3.  Levenberg  Marquardt  Nonlinear  Optimization 

The  previous  two  techniques  took  advantage  of  the  fact  that  only  one  non-linear  term, 
temperature,  is  associated  with  the  system  of  equations.  While  holding  the  temperature  constant 
and  solving  the  residual  linear  set  of  equations  is  a  valid  mathematical  technique,  if  more  teims 
were  associated  with  the  non-linearity  then  the  solution  to  the  system  of  equations  would  almost 
surely  require  a  non-linear  approach.  Modeling  emissivity  as  a  second-degree  polynomial  is  one 
method.  However,  if  bound  surface  states  or  other  material  properties  become  more  significant 
in  the  overall  IR  spectral  radiation,  the  likelihood  of  a  second-degree  polynomial  serving  as  a 
good  fit  for  emissivity  is  diminished.  Although  higher-degree  polynomials  would  provide 
adequate  fits.  However,  instrumentation  inaccuracies  increase  strongly  as  the  number  of  required 
wavelengths  increase.  Consequently,  there  is  a  high  degree  of  likelihood  that  a  non-linear  system 
of  equations  would  result  in  the  best  fitting  of  the  emissivity. 

The  Levenberg-Marquardt  is  an  elegant  method  to  accomplish  a  nonlinear  least-squares 
fit  of  the  data  to  the  model.  This  technique  can  be  used  when  the  merit  function  varies  smoothly 
between  extremes  of  the  inverse  Hessian  method.42  The  Hessian  method  defines  the  direction  of 
the  optimization  profile  (curvature).  When  the  Hessian  method  nears  the  solution  the  steepest 
ascent  method  is  used  to  find  the  global  solution.  The  Levenberg-Marquardt  method  reduces  the 
merit  function  of  a  fit  between  a  set  of  data  points  x(l...ndata),  y(l...ndata)  with  individual 
standard  deviations.  The  emissivity  coefficients  and  temperature  require  seeding  before  the 
algorithm  can  find  a  solution.  If  the  estimate  of  the  magnitude  is  too  far  from  the  true 
coefficients  or  if  the  model  is  inaccurate,  there  is  a  high  probability  that  the  algorithm  will 
become  trapped  in  a  local  minimum. 

The  simulation  showed  that  for  some  emissivity  profiles  this  algorithm  found  solutions 
for  emissivity  and  temperature  that  were  far  off  the  simulated  values  even  with  little  as  no  noise 
added  to  the  system.  For  most  of  emissivity  profiles  that  would  be  realistically  encountered  from 
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metals,  the  algorithm  performed  fairly  well  but  not  nearly  as  well  as  the  two  linear  techniques. 
As  for  the  two  linear  methods,  the  algorithm  to  accomplish  the  Levenberg-Marquardt  was 
obtained  from  Numerical  Recipes.  The  implementation  of  this  algorithm  was  a  function  call  to 
subroutines  from  the  main  program  implementing  the  three  regression  techniques  and  this  code  is 
again  documented  in  appendix  A. 


2.6.  Noise  Contribution  to  Temperature  Error 

The  contribution  of  noise  to  the  system  accuracy  can  be  determined  through  an  analysis  of 
the  extreme  data  points.  In  other  words,  if  there  are  two  equations  each  with  a  ±  error,  then  the 
maximum  error  will  occur  in  one  of  four  the  fits  to  the  two  erroneous  data  values.  For  the  system 
of  four  equations  there  will  be  error  in  optical  power  collected  for  each  of  the  spectral  regions. 
The  matrix  below  describes  the  error  associated  with  our  system  if  the  noise  is  limited  by  the 
background  and  the  read  noise. 

4]  km  4(0  Mr/P' 

ai  =  Mr)  4(0  Ih(T)  ■  i.qAi  ±  +  CT4  (2-84) 

W  ih(D  4 (O  MR I 

In  the  calculation  of  the  argument  coefficients,  an's,  one  combination  of  the  six  error 
values,  associated  with  the  ±  error  of  the  noise  will  be  the  maximum  error  in  the  calculated 
temperature.  The  algorithm  in  appendix  A  calculates  the  temperature  error  associated  with  a 
given  percentage  error  in  the  collected  optical  energy.  It  was  found  during  the  simulation  portion 
of  this  research  that  when  using  either  of  the  linear  techniques,  an  error  of  less  than  5%  (S/N  of 
20)  would  contribute  no  more  than  a  ±  5  C°  to  the  error  in  predicting  temperature.  To  determine 
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The  Levenberg-Marquardt  technique  did  not  perform  as  well  for  all  emissivity  profiles  but  nearly 
as  well  for  some,  approximately  +  30  C°.  Since  the  signal  to  noise  ratio  will  not  exceed  an  error 
of  then  5%,  as  will  be  discussed  in  the  radiometry  section  of  this  dissertation,  the  true  limitation 
of  this  technique  lies  in  the  model  for  emissivity  and  calibration  procedure  for  the  electronics. 
These  topics  will  be  detailed  in  the  technical  portion  of  this  dissertation. 


3.  Imaging  Pyrometer  Design 


As  mentioned  earlier,  multi-spectral  imaging  pyrometers  of  the  past  have  been  largely 
inaccurate.  The  primary  problem  here  is  keeping  good  calibration  as  the  numbei  of  spectial 
regions  increases.  More  than  two  or  three  spectral  measurements  have  been  shown  to  be  useless. 
Accurate  temperature  prediction  for  high-speed  deformation  application  requires  four  spectral 
measurements  in  the  same  IFOV.  It  is  possible  to  form  these  four  spectral  images  onto  a  single 
focal  plane  (Figure  3-1)  and  the  common  set  of  optics  and  electronics  will  substantially  reduce 
system  noise. 
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One  set  of  optics  is  used  for  a  single  detector  array  and  a  snapshot  camera  gathers  all  four  images 
simultaneously.  In  addition,  all  components  are  stationary  and,  consequently,  the  system 
calibrates  very  accurately  and  remains  calibrated  for  long  periods  of  time.  The  general  design  of 
the  components  of  the  system  will  be  detailed  in  this  chapter.  Discussion  will  then  move  on  to 
the  lens  and  interference  filter  designs.  A  general  optical  design  is  first  covered  in  a  step  by  step 
manner  to  show  the  evolution  of  the  high  speed  imaging  pyrometer  produced  from  this 
dissertation. 


3.1.  General  Design 

The  optics  consist  of  a  single  interference  filter  that  is  used  with  a  single  mirror  to 
separate  a  single  beam  of  light  into  four  spectral  regions.  The  concept  of  the  interference  filter  is 
that  it  changes  its  transmission  characteristics  with  the  angle  of  the  incident  optical  energy. 
Figure  4.2  illustrates  this  filter  process. 


Figure  3-2:  The  property  of  an  interference  filter  changing  its  pass  band  as  the 
angle  of  incidence  on  the  filter  changes 

The  spectral  energy  not  in  the  transmission  band  of  the  filter  will  see  a  near  100  per-cent 
reflection  at  the  surface  of  the  filter.  As  the  field  moves  to  angles  of  higher  incidence,  the  filter 
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responds  to  shorter  wavelengths.  Consequently,  as  the  filter  is  tilted  the  light  that  is  passed  for 
each  of  the  fields  shifts  to  shorter  wavelengths  (Figure  3-3). 


Filter 


Figure  3-3:  Interference  filter  behavior  when  tilted 

It  can  be  seen  that  the  on-axis  field  moves  to  pass  green  light  where  the  off-axis  field  moves  to 
pass  blue  light.  In  an  image  forming  system  with  the  filter  in  a  collimated  space,  the  filtei  would 
respond  as  indicated  Figure  3-4. 
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If  the  object  is  a  white  light  source,  the  image  color  will  vary  from  long  wavelength  (green  light) 
on  axis,  to  shorter  wavelengths  (blue  light)  off  axis.  The  next  step  is  to  create  multiple  images 
each  with  a  unique  spectral  content.  This  can  be  done  by  adding  a  mirror  to  the  setup  to  the  left 
of  the  tilted  interference  filter,  figure  3-5. 


67 


Figure  3-5:  Lens  Configuration  with  tilted  filter  and  mirror  to  return  reflected  energy 
subsequently  returning  to  the  filter  at  an  incremental  angle 

The  solid  lines  correspond  to  the  straight-through  component  of  light.  The  light  that  is  not 
passed  through  the  filter  is  reflected  at  the  filter  surface,  and  this  light  is  represented  by  the 
dashed  lines.  Subsequently,  this  reflected  light  experiences  a  second  reflection  at  the  mirror. 
The  total  deviation  is  twice  the  separation  angle  between  the  filter  and  the  mirror.  If  this 
separation  is  done  in  collimated  space,  the  angular  change  will  alter  the  position  that  a  field 
focuses  onto  a  given  focal  plane  without  any  defocus.  To  have  contiguous  images,  the  separation 
angle  must  be  one  half  the  maximum  field  angle  subtended  by  the  object.  Additional  images  can 
be  formed  from  additional  reflections  at  the  filter  surface  and  each  image  will  have  a  unique 
spectral  content  from  passing  through  the  interference  filter  at  different  angles.  By  taking 
advantage  of  these  spectral  transmission  and  reflection  variations  with  incident  angle,  the  need 
for  multiple  spectral  filters  is  eliminated.  In  addition,  multiple  images  are  formed  from  the  same 
beam  of  light  making  the  spectrometer  very  fast. 

In  the  past,  imaging  spectrometers  have  required  significant  time  to  scan  the  spectral 
energy  onto  a  charge-coupled-device  (CCD)  or  other  readout  device.  This  would  take  the  form 
of  a  moving  prism,  diffraction  grating  or  by  moving  filters  in  and  out  of  the  optical  path.  The 
application  of  predicting  temperature  for  a  ballistic  projectile  moving  at  1000  meters  per  second 
or  more  with  1  mm  resolution  will  require  the  spectrometer  to  gather  all  data  in  approximately 
one  microsecond.  The  advantage  of  using  this  new  technique  is  that  the  speed  of  the  system  is 
only  limited  by  the  integration  time  of  the  camera.  Spectral  energy  is  essentially  picked  off  the 
optical  path  without  losing  the  energy  in  adjacent  spectral  regions.  This  creates  a  system  that  is 
very  light  efficient  which  is  extremely  important  for  high  speed  applications.  Now  that  the 
general  high-speed  instrument  design  has  been  detailed,  the  filter  and  lens  components  can  be 
designed. 
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3.2.  Interference  Filter  Design 


In  classical  optical  systems,  filters  are  designed  to  have  small  shifts  in  spectral 
characteristics  with  angle  of  incidence.  The  objective  is  to  have  all  components  of  the  cone  of  a 
focusing  beam  see  the  same  pass  band.  The  application  in  this  dissertation  takes  advantage  of 
shifting  the  band  pass,  therefore,  the  classical  filter  are  not  useful.  The  root  of  the  design  is  in 
separating  each  of  the  images,  spectrally,  as  much  as  possible.  It  is  advantageous  to  separate  the 
spectrum  equally  among  the  four  filtered  images  so  that  irradiance  remains  somewhat  constant 
across  the  array.  To  improve  the  interference  filter,  several  phenomena  can  be  taken  advantage 
of  in  the  design.  The  first  design  that  was  looked  at  involved  changing  the  fundamental 
deposition  of  the  multi-layer  coatings.  While  classical  designs  attempt  to  minimize  spectral  shift, 
it  is  also  feasible  to  maximize  this  shift.  After  accomplishing  a  few  preliminary  designs  using 
the  MACLEOD  thin  film  design  software  and  contacting  filter  vendors,  it  became  apparent  that 
the  cost  of  manufacturing  this  interference  filter  was  prohibitive. 

The  cost  of  the  filter  can  be  reduced  by  using  a  linear  variable  filter.  Here,  the  shift  in 
field  position  on  the  filter  is  used  instead  of  relying  solely  on  the  angular  shift  to  achieve  all  the 
spectral  separation  (Figure  3-6).  The  multiple  optical  paths  traverse  the  interference  plane  in 
different  locations  and  with  the  proper  amount  of  shift  (for  this  case  somewhat  non-linear)  all  the 
spectrum  in  the  three  to  five  micron  band  can  be  used.  While  this  design  worked  well  in 
simulation,  it  was  still  too  expensive  to  manufacture  and  other,  less  expensive  designs  had  to  be 
considered. 
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Figure  3-6:  Using  a  linear  variable  filter  to  get  the  proper  amount  of 
wavelength  shifts  out  of  each  spectral  image. 


The  technique  that  was  selected  took  advantage  of  optical  path  shift  on  the  interference 
plane.  Discrete  filters  at  four  locations  of  the  plane  were  used  to  accomplish  spectral  separation. 
It  can  be  seen  in  Figure  3-7  that  each  optical  path  is  separated  in  the  filter  plane.  This  makes  it 
easy  to  control  the  spectral  content  of  each  image  on  the  focal  plane. 


First  surface  of  imaging  lens 


3D  LRYOUT  I 

XMRGTNG  SPECTROMETER/P YROMETR Y 

FRI  TUN  1H  1996 

Figure  3-7:  Discrete  filter  approach  for  effective  filter  separation 
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While  it  is  desirable  to  manufacture  all  filters  on  the  same  substrate,  again  the 
manufacturing  is  cost  prohibitive.  Consequently,  individual  off-the-shelf  filters  were  obtained 
from  various  venders  and  adhered  to  a  common  substrate.  Experimental  results  will  be  presented 
in  the  results  section  of  this  dissertation.  Now  that  the  interference  filter  has  been  designed  the 
optical  design  can  be  completed. 


3.3.  Optical  Lens  Design 

The  work  on  the  lens  design  started  with  the  paraxial  characteristics  of  the  system.  Since 
the  munition  test  setup  is  fairly  compact  the  system  was  designed  for  finite  conjugates  with  the 
nominal  object  distance  at  0.5  meters.  The  order  from  left-to-right  is  the  objective,  collimation. 
and  imaging  lenses.  The  intermediate  image,  just  after  the  objective  lens,  is  the  field  slop,  and 
the  aperture  stop  is  at  the  collimation  lens. 

The  paraxial  characteristics  are  driven  by  the  field-of-view  of  the  system  and  the  number 
of  spectral  images  designed  to  fall  on  the  focal  plane.  In  addition,  since  the  interference  filter 
design  require  the  fields  for  the  spectral  images  to  be  separated  at  the  filter  plane,  the  size  of  the 
imaging  lens  must  to  be  larger.  This  problem  will  not  be  encountered  when  the  filter  is  designed 
to  spectrally  separate  the  images  using  only  angular  shift.  After  this  filter  is  manufactured  the 
paraxial  characteristics  can  be  defined  for  the  new  system. 

To  start  the  design  of  the  system  of  this  dissertation  the  FOV  first  needed  to  be 
determined.  The  Amber  focal  plane  is  8mm  on  a  side  and  a  25mm  focal  length  imaging  lens  was 
chosen  so  that  the  four  images  were  separated  by  5°  on  the  focal  plane.  The  objective  and 
collocation  lens  combination  is  essentially  a  Keplerian  telescope.  There  are  two  reasons  for  this 
design.  First,  the  collimating  lens  needs  to  be  reasonably  small.  If,  for  instance,  the  collimating 
lens  were  25  mm  in  diameter,  the  separation  between  the  filter  and  the  mirror  would  be  about 
30mm  for  an  initial  incident  angle  of  45°.  For  no  vignetting  of  the  off-axis  fields,  the  focal 
length  of  the  imaging  lens  would  need  to  be  on  the  order  of  25  mm.  This  drives  the  F/#  of  the 
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imaging  lens  to  F/l .  This  is  very  fast  and  difficult  to  design.  To  increase  the  F/#,  the  collimating 
lens  was  made  smaller  so  that  the  image  space  ray  bundle  is  more  compact.  In  addition,  the 
telescope  allows  a  reduction  in  the  field  of  view  so  that  there  is  some  stand  off  from  the 
experiment.  For  the  munition  experiments  that  are  performed,  the  objects  are  no  larger  than  25 
mm.  At  the  500  mm  object  distance,  makes  the  FOV,  approximately  1°  in  object  spac'e.  The 
collimating  lens  produces  a  5°  FOV  making  the  magnification  power  5.  In  addition,  the 
telescope  has  an  intermediate  image  plane,  allowing  the  placement  of  a  field  stop.  This  will 
insure  that  there  is  no  overlap  of  adjacent  out-of-field  spectral  images. 

The  collimating  lens  was  chosen  to  operate  at  F/2  since  an  F/2  system  is  reasonable  to 
design.  After  running  through  a  paraxial  ray  trace,  it  was  found  that,  in  order  for  the  imaging 
lens  to  also  operate  at  F/2,  the  collimating  lens  diameter  needed  to  be  approximately  5  mm.  Now 
the  paraxial  characteristics  can  be  computed. 

focal  length  of  objective  =  50  mm 
focal  length  of  collocating  lens  =  1 0  mm 
focal  length  of  imaging  lens  =  25  mm 

Lens  design  codes  do  not  allow  multiple  transmissions  through  the  single  interference 
filter  while  keeping  track  of  the  reflected  component.  Consequently,  this  lens  design  was  done  in 
multiple  configuration  format.  There  are  four  configurations  corresponding  to  the  four  passes 
through  the  filter.  The  reflected  component  from  the  filter  is  accomplished  by  molding  the  filter 
as  a  mirror.  After  the  light  has  gone  through  the  appropriate  number  of  reflections  the  filter  takes 
on  its  true  characteristics.  The  figures  below  show  a  ray  trace  for  the  paraxial  system  discussed 
above. 
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Figure  3-8:  Paraxial  lens  design  for  the  fourth  reflected  ray  through  the 
interference  filter 
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Figure  3-9:  Paraxial  lens  design  for  the  third  reflected  ray  through  the 
interference  filter 
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Figure  3-10:  Paraxial  lens  design  for  the  second  reflected  ray  through  the 
interference  filter 
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Figure  3-11:  Paraxial  lens  design  for  the  first  reflected  ray  through  the 
interference  filter 
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After  the  paraxial  characteristics  of  the  system  were  defined,  some  real  lenses  could  be 
put  together.  The  analysis  was  started  from  an  aberration  point-of-view  as  to  how  each  lens 
should  be  balanced.  At  this  point,  it  was  known  that  off-the-shelf  lenses  would  be  used. 
Consequently,  the  formal  lens  design  was  accomplished  in  order  to  determine  the  quality  of  each 
of  the  lens  components.  Another  reason  for  completing  this  lens  design  is  that,  in  general,  the 
designs  for  off-the-self  lenses  are  proprietary  and  are  not  easily  obtained.  The  following  design 
can  be  simplified  somewhat  if  it  was  designed  as  a  system.  The  cost  of  assembling  these  lens 
components  for  the  prototype  system,  however,  would  be  much  higher  than  if  off-the-shelf 
components  were  used. 

For  the  objective  lens,  it  is  important  that  the  image  that  occurs  at  the  field  stop  be  free 
from  aberrations  that  shift  the  image  location.  The  concern  is  weather  or  not  the  collimating  lens 
truly  collimates  all  of  the  fields.  Each  of  the  four  images  will  experience  different  optical  path 
lengths,  causing  the  aberrations  to  vary  at  the  imaging  lens  in  a  way  not  predicted  by  third-order 
or  higher  aberrations.  At  the  very  least,  the  collimating  lens  needs  to  compensate  for  objective 
aberrations.  The  objective  lens  was  started  with  a  single  lens  of  germanium  and  best  shape 
factor  to  minimize  spherical  aberration  (sph3)  for  a  50mm  focal  length  lens  was  found.  Then,  a 
second  lens  of  lower  index  was  added  and  balanced  for  sph3  and  longitudinal  color  by  varying 
only  the  surface  curvatures.  Subsequently,  adding  two  more  elements  made  a  symmetric  system 
/surface  curvatures  resulting  in  a  near  zero  sph3,  with  color,  coma,  astigmatism  aberrations, 
petzvall  curvature,  and  distortion  minimized.  The  lens  optimized  to  the  50  mm  focal  length  very 
well,  however,  the  higher  order  aberrations  degraded  the  spot  diagram  performance.  The  default 
merit  function  was  used  to  optimize  the  RMS  spot  diagram  performance.  Less  than  2-micron 
radius  geometric  precision  was  achieved. 

The  imaging  lens  was  designed  in  the  same  way  as  the  objective  lens.  Using  these  two 
designs  the  collimating  lens  design  was  started  by  first  letting  the  surfaces  of  a  triplet  vary  until  a 
reasonable  solution  was  attained.  Again,  the  same  procedure  was  followed  as  for  the  first  two 
lenses,  first  taking  care  of  the  third  order  aberrations,  and  then  using  the  default  merit  function  to 
minimize  the  RMS  spot  size.  All  the  surfaces  and  thickness  of  the  system  were  allowed  to  vary 
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giving  some  reasonable  bounds.  A  diffraction-limited  system  was  reached,  and  tolerance  and 
estimates  from  venders  were  attained.  The  lens  system  is  shown  in  Figure  3-12. 
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Figure  3-12:  Lens  design  for  the  fourth  reflected  ray  through  the  interference 
filter 


As  was  mentioned  earlier,  to  build  these  lenses  from  scratch  would  cost  a  great  deal  of 
money.  Consequently,  for  this  research,  off-the-shelf  components  were  used  to  accomplish  the 
design.  After  this  system  is  calibrated  as  discussed  in  the  next  chapter,  the  four  colors  are 
collected  using  an  IR  focal  plane  array.  For  high-speed  applications,  it  is  both  important  to  use  a 
sensitive  detector  as  well  as  to  get  as  much  energy  through  the  system  as  possible,  for  the 
detector,  it  is  evident  that  an  Indium  Antinomide  (InSb)  snapshot  infrared  array  would  give  the 
best  performance  for  this  application.  The  next  step  is  to  determine  the  available  optical  energy 
and  analyze  the  limitations  of  radiometric  detection  of  thermal  energy  for  each  filter. 
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3.4.  Radiometry 


To  obtain  the  required  scene  information,  the  signal-to-noise-ratio  must  be  determined. 
This  can  be  done  by  looking  at  one  filter  at  a  time  to  estimate  the  accuracy  of  measuring  the 
temperature  with  an  assumed  surface  emissivity.  After  these  signal  characteristics  are 
determined,  the  limitations  of  this  system  can  be  found.  Starting  with  the  assumption  that  the 
surface  acts  as  a  graybody,  it  will  follow  Planck’s  radiation  law: 

_  2a-, 

q  AV(e'-l)  ’  (3_1) 


in  photons/second-area-solid  angle,  where 


c  =  velocity  of  light  =  2.9979  x  10^  m/s, 

A.  =  wavelength  =  3.0  -3.4,  3. 3-3. 8,  3. 8-4.2,  4.5-5.0  pm, 

x  =  hc/kBAT  =  14388  pm-K  /  AT, 

h  =  Planck's  constant  =  6.624  x  10'34  Js, 

kB  =  Boltzmann's  constant  =  1.38  x  10"23  J/K,  and 

T  =  temperature  =  525  K  (~  250  °C), 

s  =  emissivity  =  0.2. 

Given  these  parameters,  Lq  =  2.8  x  \0\7  photons/sec-cm^-sr.  The  photon  flux  from  the  surface 

of  emissivity,  e,  and  area,  Ap ;  to  the  collection  optics  of  area,  A0 ,  separated  by  distance,  R,  is 


Since  radiance  is  conserved,  it  is  equivalent  to  write 
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(3-3) 
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For  the  F/#=2.3  system  developed  in  this  effort,  with  a  detector  area  of  .000009  cm^  and  a 
nominal  emissivity43  of  0.2  (assuming  no  wavelength-dependence),  Nq  =  3.3  x  10^ 

photons/second.  To  calculate  the  total  signal  at  the  detector,  a  total  system  transmittance,  t  =  0.5, 
an  integration  time,  tint  ~  5  ps,  as  well  as  a  detector  quantum  efficiency  of  p  =  0.80,  (typical  for 

our  InSb  detector)  are  included.  This  gives  a  total  signal  of  6.6  x  10^  e"  per  InSb  detector.  As 
this  situation  is  clearly  background  noise  limited,  the  total  noise  in  electrons  will  be  the  quadratic 
sum  of  the  photon  noise  and  the  read  noise  ( (Tread  is  taken  here  to  be  120  e'): 


=  V 8.44  x  105  + 1202  =  916e" 


(3-4) 


and  the  signal  to  noise  ratio  is  6.6  x  10^/916=720  or  about  0.1%.  This  is  a  good  result  because 
the  algorithm  used  to  fit  temperature  remain  accurate  to  a  few  degrees  even  if  there  is  much 
larger  amounts  of  noise.  The  high-speed  pyrometer  was  built  and  tested  for  accuracy,  the  details 
of  which  is  presented  in  the  next  chapter. 
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4.  High  Speed  Pyrometer  Fabrication 


The  hardware  discussed  in  the  last  chapter  was  assembled  using  primarily  off-the-self 
components.  In  this  chapter,  the  mechanical  assembly  for  the  high-speed  instrumentation  will  be 
discussed.  In  addition,  the  experimental  setup,  calibration  instrumentation,  and  test  procedure 
will  be  detailed.  The  testing  of  this  system  consists  low  speed  validation  tests  to  characterize 
camera  performance. 


4.1.  Multi-spectral  Camera 

The  first  step  to  assembling  the  camera  components  was  to  select  a  focal  plane.  In  the 
munition  test  environment,  the  objects  that  are  being  analyzed  are  moving  at  MOCH-one  or 
greater  and  stop-motion  resolution  is  required  to  about  one-millimeter.  Consequently,  the 
requirement  is  that  the  camera  shutter-down  to  one  microsecond  or  better.  After  a  survey  of  the 
IR  cameras  that  could  potentially  meet  this  requirement,  the  Amber  AE173  snapshot  was 
selected.  The  optics  for  the  front  end  of  the  camera  was  built  as  discussed  in  the  previous 
chapter,  and  is  shown  in  Figure  4-1 . 
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Imaging  lens 


Figure  4-1:  High  speed  pyrometer  camera  configuration.  The  lens  assembly  is  at  a  slight 
angle  from  the  front  face  of  the  camera,  which  is  done  to  incorporate  the 
unique  lens  design 


The  lens  assembly  is  mounted  to  the  Amber  Camera  with  four  bolts  and  can  be  removed 
easily  for  component  inspection.  All  components  hard  mount  and  the  only  part  that  moves  is  the 
objective  lens  for  focusing.  The  objective  lens  mounts  in  the  same  clamp  as  the  front  end  of  the 
Amber  camera  making  it  simple  to  interchange  lenses  for  alternative  applications.  After 
detaching  the  filter  assembly  from  the  camera  face,  the  imaging  lens  can  be  seen.  Figure  4-2 
shows  how  the  imaging  lens  fits  in  the  standard  Amber  front  lens  mount.  The  filter  assembly 
mounts  in  existing  lens  mounting  holes,  thus  avoiding  modification  of  the  camera. 
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Figure  4-2:  Disconnection  of  filter  assembly  from  main  camera  body.  The  assembly  is  held 
to  the  front  face  of  the  camera  via  four  mounting  bolts.  The  imaging  lens  fits  in 
the  standard  Amber  mount  in  the  front  of  the  camera.  Both  the  objective  and 
imaging  lenses  were  manufactured  by  DIOP. 

After  removing  the  filter  assembly  from  the  camera,  inspection  of  the  interference  filters 
is  made  by  looking  through  the  back  end.  Just  on  the  other  side  of  the  filters  the  mirror  can  also 
be  inspected.  The  filter/mirror  combination  is  secured  to  the  tubular  housing  with  two  Phillips 
head  screws.  A  picture  of  the  back  end  view  is  shown  in  figure  4-3. 


Figure  4-3:  Interference  filter  and  mirror.  The  mirror  that  folds  the  multiple  beams  back 
towards  the  interference  filters  can  be  seen  through  the  small  hole  directly 
above  the  filters.  On  the  other  side  of  this  assembly  are  the  collimating  lens 
and  the  field  stop. 


The  field  stop  may  be  seen  by  turning  the  filter  assembly  around  and  detaching  the 
objective  lens.  The  field  stop  consists  of  two  razor  blades  fashioned  inside  the  front  end  of  the 
filter  tube  (Figure  4-4).  The  objective  lens  connects  to  a  mounting  flange  on  the  front  of  the  filter 
tube  so  that  an  image  will  form  on  the  field  stop.  Focusing  for  the  entire  system  is  done  by  way 
of  the  “built-in”  focusing  property  of  the  objective  lens. 
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Figure  4-4:  Field-stop  and  collimating  lens.  The  collimating  lens  is  between  the 
interference  filters  and  the  field-stop  and  consequently  cannot  be  seen  in  either 
view.  The  field-stop  is  made  from  two  razor  blades. 


To  validate  the  pyrometer  performance  fairly  accurately,  control  was  needed  over  the 
temperature  of  the  test  metal.  This  could  not  be  done  with  any  accuracy  at  high  speeds, 
consequently,  low  speed  validation  was  implemented.  The  temperature  of  the  sheet  metal  was 
controlled  by  pasting  the  metal  samples  to  an  extended  blackbody  using  heat  sink  compound. 
The  arrangement  is  shown  in  Figure  4-5. 
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camera  calibration  and 
two  point  correction  : 


Figure  4-5:  Experimental  test  item  and  black  body  used  in  system  calibration 


With  temperature  accurately  known,  the  emissivity  could  be  accurately  determined.  This  setup 
was  also  used  to  predict  temperature  for  the  four  spectral  measurements.  The  temperature 
predicted  was  accurately  compared  to  the  true  temperature  reading  of  the  blackbody  controller. 
In  addition,  a  thermocouple  was  mounted  to  the  front  surface  of  the  metal  sample  so  that  there 
was  an  additional  check  for  true  temperature. 
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4.2.  Calibration: 


After  the  pyrometer  camera  components  were  assembled,  the  next  step  was  to  design  a 
calibration  procedure.  The  calibration  is  broken  up  into  two  parts:  calibration  of  the  spectral 
transmission  of  each  of  the  four-color  images  and  calibration  of  the  gain  and  offset  of  the  InSb 
detector.  Several  techniques  for  determining  the  filter  transmission  characteristics  were 
examined  including  imaging  an  extended  blackbody  and  gathering  spectral  data  at  many  different 
temperatures.  As  in  the  problem  of  determining  temperature  from  a  system  of  equations 
generated  from  multiple  spectral  measurements,  a  system  of  equations  from  multiple  temperature 
measurements  can  also  be  solved  to  give  the  spectral  transmission  for  each  of  the  filtered  images. 
It  turned  out  that  the  matrix  describing  the  system  of  equations  is  highly  singular  and  as  a  result 
predicting  the  filter  transmission  is  not  very  accurate.  Consequently,  the  brute  force  technique  of 
sweeping  a  monochromator  through  wavelength  was  chosen.  The  setup  is  shown  in  Figure  4-6. 


An  Acton  Research  Inc.  SpectroPro  500  monochromator  was  controlled  with  the  system 
computer  and  stepped  through  500  values  of  wavelength  from  2700  to  5200  nanometers  in  5 
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nanometer  increments.  An  integrating  sphere  was  used  at  the  output  of  the  monochromator  to 
flat  field  the  detector.  Using  the  monochromatic  data  and  an  extended  blackbody  the  calibration 
of  the  system  filters  was  done  in  two  steps.  The  first  step  is  to  determine  the  spectral  response 
(per-cent  transmission)  of  each  detector  on  the  array.  The  detector  responsivity,  transmission  of 
the  optics,  transmission  of  the  filters,  transmission  of  the  integration  sphere,  transmission  of  the 
atmosphere,  and  the  spectral  characteristics  of  the  monochromator  all,  affect  thev  perceived 
optical  energy,  Eq.  4-1. 

=  Top,  {fy  filters  tsph  M7  atm  {A)Ln,ono  W^de,  (4‘ 1 ) 

Determining  each  of  the  above  system  parameters  is  quite  involved.  However,  since  the 
filter  transmission  can  be  determine  from  the  relative  response  for  each  of  the  four  filters  we  only 
need  to  determine  the  relative  signal  amplitude  before  and  after  the  filters  are  placed  in  the 
system.  Equation  4-2  represents  the  percentage  transmission  for  each  of  the  spectral  regions.  To 
relate  the  energy  on  the  detector  to  the  IR  energy  leaving  the  metal  surface  the  second  step  of 
calibration  needs  to  be  performed,  gain  and  offset. 

^w/ftllersi^)  _  {^)T  filters  sph  {^)T  atm  (2.)Qfldel  _  ^  (4-?) 

In  order  to  get  the  correct  values  for  a  material's  radiant  exitance,  the  gain  and  offset  of 
the  entire  system  needs  to  be  determined.  Since  the  data  from  each  detector  is  digitized,  the 
output  is  in  analog  digital  units  (ADU).  Neglecting  atmospheric  transmission  the  equation  is 

ADU  =  OFFSET +z^(X)Lim(X}nAt.s(X)G.  (4-3) 

where  G  is  the  gain  in  ADU’s  per  watt.  If  the  systems  solid  angle  is  folded  into  the  gain  there  are 
only  two  unknown  constants  in  the  system,  gain  and  offset.  Collecting  data  for  an  extended 
blackbody  at  a  high  and  a  low  temperature  will  form  two  equations  and  the  gain  and  offset  can  be 
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calculated.  Gain  and  offset  can  only  be  done  after  the  filter  transmission  is  characterized. 
Innately  there  was  some  problem  in  this  calibration. 

The  manufacture’s  power  supply  for  the  monochromator  light  source  proved  to  be  too 
noisy  to  accomplish  an  accurate  calibration.  In  addition,  the  power  in  the  laboratory  had  as  much 
as  +10  volts  fluctuation  adding  additional  noise.  To  overcome  these  problems  the  power  supply 
was  replaced,  and  testing  was  done  after  normal  duty  hours,  when  the  AC  power  was  more 
predictable.  Figures  4-7  and  4-8  show  the  fluctuation  before  and  after  the  above  correction  were 
made. 


Figure  4-7:  Error  of  monochromator  from  calibration  run  to  calibration  run.  This 
data  was  taken  during  normal  working  hours  and  with  the  manufacture's 
power  supply. 
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Wavelength  in  microns 

Figure  4-8:  Error  of  monochromator  from  calibration  run  to  calibration  run. 

This  data  was  taken  after  normal  working  hours  and  with  a 
Power  One  Inc.,  F15-15-A  power  supply. 

Along  with  the  power  fluctuations,  the  room  temperature  also  fluctuated.  This  caused  the 
temperature  of  the  optics  to  fluctuate,  and,  consequently,  the  DC  offset  of  the  background  signal 
to  vary.  This  problem  was  overcome  by  masking  off,  at  the  field  stop,  one  line  of  the  array  from 
external  radiation  (Figure  4-9).  These  detectors  could  only  see  the  radiation  of  components  to  the 
right  of  the  field  stop.  Since  all  detectors  on  the  array  will  fluctuate  in  proportion  to  this 
background,  the  masked  detectors  were  monitored  to  compensate  for  the  thermal  variation  of  the 
optics  during  calibration. 
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Figure  4-9:  Detector  masking  approach  to  remove  the  affect  of  thermal 
variation  caused  by  air-conditioning  in  laboratory.  These 
fluctuations  were  as  much  as  2%  over  a  single  calibration  run. 


After  the  fluctuation  in  background  was  normalized  the  two  calibration  runs  were  made,  one  with 
the  filters  and  one  without.  Using  equation  4-2  the  transmission  for  each  of  the  images  could  be 
computed  (figure  4-10). 


Figure  4-10:  The  filter  bandpass  response  for  each  of  the  four  filters.  The 
first  pass  corresponds  to  the  first  filter  on  the  left,  the  second 
pass  is  the  second  filter,  and  so  forth.  Notice  that  each  of  the 
filters  peak  response  is  gradually  decreasing.  This  corresponds 
primarily,  to  out-of-band  absorption. 
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In  the  design  for  the  filter  set  above,  the  goal  was  to  use  as  much  of  the  overall  3  to  5  micron 
spectrum  as  possible  while  staying  away  from  the  C02  notch.  It  can  also  be  seen  in  figure  4-10 

that  for  the  transmission  of  images  2,  3,  and  4  there  is  a  degradation  in  the  peak  transmission. 
This  is  due  to  out  of  band  absorption  that  can  be  removed  with  proper  filter  design. 


4.3.  Experimental  Setup 

After  the  pyrometer  camera  components  were  assembled  and  the  filter  transmission  was 
characterized,  the  next  step  was  to  accomplish  low  speed  validation  testing  to  see  if  the  theory 
agrees  with  reality.  The  setup  is  shown  in  Figure  4-11. 


Figure  4-11:  Test  and  calibration  setup 


The  blackbody  containing  the  test  samples  and  the  blackbody  used  for  two  point 
correction  were  mounted  on  a  rail  so  that  they  may  be  slid  in  and  out  of  the  optical  path  quickly. 
Three  types  of  metal  samples  were  mounted  to  the  extended  blackbody:  copper,  aluminum  and 
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steel.  The  camera  system  collected  numerous  temperature  data  points  for  each  of  the  samples  to 
insure  repeatability,  the  results  of  which  will  be  presented  later.  The  support  electronics  to  the 
Amber  camera  system  is  the  general  purpose  ProView  unit  (Figure  4-12).  The  ProView 
electronics  provided  camera  clocking  and  A/D  conversion  of  the  analog  data.  After  the  data  is 
acquired,  the  ProView  relays  the  image  information  to  the  system  computer  and  to  its  hard  drive. 
The  ProView  unit  comes  with  its  own  language  (TDF)  that  is  used  to  set  system  parameters  such 
as  timing,  gain,  offset,  etc.  The  code  for  which  is  located  in  appendix  B. 


Figure  4-12:  ProView  electronics  package 


4.4.  Validation 

There  were  two  means  used  to  validate  the  calculated  temperature.  First  the  temperature 
of  the  sample  on  the  blackbody  was  known  (internal  Thermocouple)  from  the  control  panel.  This 
temperature  was  checked  at  the  metal  surface  using  an  external  thermocouple.  In  addition  to  the 
thermocouples,  a  commercial  pyrometer  manufactured  by  Agema  Inc.,  was  used  to  compare 
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results.  This  type  of  pyrometer  uses  the  user  input  values  for  emissivity  to  compute  temperature 
(picture  shown  in  Figure  4-13).  This  camera  was  used  primarily  as  a  comparison  of  this  new 
technology  to  a  typical  commercial  product.  When  the  emissivity  calculated  by  the  high-speed 
pyrometer  was  used  in  the  AGEMA  camera  the  temperature  prediction  was  highly  accurate. 
However  this  system  relies  on  accurate  prediction  of  an  average  emissivity  over  the  3  to  5  band. 


Figure  4-13:  The  AGEMA  Thermo-vision  IR  camera  used  for 
validation. 


The  validation  for  the  emissivity  following  a  polynomial  dependence  was  presented 
earlier  and  the  purpose  of  the  above  instrumentation  was  to  show  that  the  system  could  predict 
temperature  accurately.  If  the  actual  model  for  emissivity  were  different  from  a  second-degree 
polynomial,  the  temperature  prediction  would  have  increasing  error  with  increasing  error 
between  the  two  models. 


4.5.  Results 
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After  many  measurements  were  made  to  insure  repeatability,  the  results  were  plotted  in 
Figure  4-14  through  4-17.  It  is  interesting  to  note  that  all  samples  followed  the  same  trend.  As 
temperature  was  increased,  the  curvature  for  the  emissivity  profile  flattened  out,  and  the  average 
emissivity  rose.  All  the  emissivity  profiles  were  concave.  The  results  for  temperature  prediction 
are  shown  in  tables  4-1  through  4-4  for  each  of  the  metals. 
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Figure  4-14:  Plot  of  the  emissivity  of  a  heated  piece  of  aluminum  at 
various  temperatures.  The  surface  of  the  sample  was 
prepared  using  200  grit  sandpaper. 


Figure  4-15:  Plot  of  the  emissivity  of  a  heated  piece  of  copper  at 
various  temperatures.  The  surface  of  the  sample  was 
prepared  using  200  grit  sandpaper. 
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Figure  4-16:  Plot  of  the  emissivity  of  a  heated  piece  of  steel  at  various 
temperatures.  The  surface  of  the  sample  was  prepared 
using  200  grit  sandpaper. 


Heavily  Rusted  Steel 


Wavelength  (>  m) 


120  degrees  centigrade 
160  degrees  centigrade 
200  degrees  centigrade 
250  degrees  centigrade 


Figure  4-17:  Plot  of  the  emissivity  of  a  heated  piece  of  heavily  oxidized 
steel  at  various  temperatures. 


Table  4-1:  Temperature  error  of  Aluminum  Sample 
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167 
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+4 
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+10 

278.1 
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-8.1 

Table  4-2:  Temperature  error  of  Copper  Sample 
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Table  4-3:  Temperature  Error  of  Steel 
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Nearly  all  the  calculation  errors  were  within  ±5  °C.  At  high  temperatures  there  was  a 
significant  error.  This  was  possibly  due  to  the  fact  that  the  heat  sink  compound  is  not  designed 
for  use  at  temperatures  in  excess  of  200  °C. 

The  singular  value  decomposition'  code  was  used  to  reconstruct  tempeiatuie  and 
essentially  the  same  results  were  attained  as  for  Gaussian  elimination.  There  was,  however,  a 
slight  degradation  in  predicting  temperature,  and  the  SVD  method  achieved  an  accuracy  of  about 
±8  °C.  This  additional  error  is  most  probably  due  to  rounding  errors  associated  with  the  more 
complex  algorithm  execution.  I  believe  with  some  careful  work,  the  code  could  be  improved  to 
attain  the  same  accuracy  as  Gaussian  elimination. 

For  the  Levenberg  and  Marquardt  nonlinear  optimization,  again,  the  accuracy  was 
somewhat  degraded  form  the  Gaussian  elimination  method,  seeing  an  accuracy  of  ±25  degrees 
error.  In  addition,  the  error  was  highly  subject  to  the  accuracy  of  the  seeding  algorithm.  It  is 
possible  that  the  code  is  becoming  hung  up  in  local  minimums.  Another  problem  is  that  with  flat 
emissivity  profiles,  the  algorithm  does  not  work  very  well  at  all.  This  was  seen  for  the  heavily 
oxidized  sample,  where  the  error  was  several  hundred  °C. 
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5.  Conclusions 


The  determination  of  the  temperature  of  an  object  through  the  use  of  only  passive 
radiation  has,  in  the  past,  been  met  with  ambiguous  results.  The  primary  problem  is  to  separate 
the  combined  effects  of  emissivity  and  temperature  on  the  spectral  radiance.  This  research  has 
shown  that  determining  the  temperature  of  metals  is  achieved  with  a  model  for  emissivity  as 
second  order  polynomial  dependence  with  wavelength.  This  model  was  developed  from  a 
physics-based  approached  with  the  understanding  that  the  free  electrons,  Fresnel  surface 
reflections  and  surface  topography  have  a  combined  contribution  to  emissivity.  This  model  was 
validated  and  shown  to  be  accurate  for  the  metal  found  in  weapon  testing.  This  physical  model 
allows  a  flexible  way  to  determine  a  metals  emissivity  ,  and  as  material  parameters  change,  the 
new  wavelength  dependent  model  can  be  quickly  determined. 

It  is  both  important  to  assume  an  accurate  functional  form  for  the  emissivity,  as  well  as 
develop  instrumentation  that  can  accurately  acquire  the  radiance  data.  Since  emissivity  and 
temperature  have  functional  dependencies  on  wavelength,  the  use  of  multi-spectral 
instrumentation  is  predicated.  To  measure  emissivity  and  temperature  to  the  required  accuracies, 
a  special  high-speed  IR  imaging  spectrometer  was  needed.  Using  the  knowledge  that  the 
emissivity  takes  the  form  of  a  second  degree  polynomial,  in  addition  to  assuming  that  the  object's 
temperature  is  in  steady  state  and  thus  its  temperature  dependent  spectral  behavior  follows  the 
Planck  distribution,  an  IR  spectrometer  capable  of  measuring  four  spectral  regions  was 
developed  and  tested.  Utilizing  a  common  set  of  optics  and  a  single  interference  filter  allowed 
efficient  simultaneous  collection  of  all  spectral  information  extremely  fast.  This  pyrometer  was 
shown  to  remain  calibrated  for  long  periods  of  time,  and  having  no  moving  parts,  this  system  is 
very  rugged  and  normal  handling  does  not  alter  calibration.  The  pyrometer  proved  to  be  accurate 
within  +5  degrees,  for  temperatures  ranging  from  100  to  200  degrees  centigrade.  Due  to  the 
laboratory  experimental  setup,  it  was  not  possible  to  test  the  system  at  higher  temperatures; 
however,  it  is  not  anticipated  that  the  system  would  degrade  much,  or  even  at  all,  until  the 
material  started  to  melt  or  burn. 
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Up  to  this  point,  experimental  verification  for  munition  models  has  only  been  performed 
for  object  deformation,  crystal  lattice  realignment,  and  total  heat  energy  conversion.  Additional 
experimental  measurements  of  object  radiation  profiles  must  be  performed  to  accurately  validate 
the  manner  in  which  temperature  is  incorporated  into  these  new  models.  Specifically,  point-to- 
point  positions  on  an  object  must  be  mapped.,  first  in  terms  of  that  object  s  radiance,  and 
subsequently  in  terms  of  the  object’s  emissivity  and  temperature.  Specific  characteiistics  of  the 
energetic  environment  give  rise  to  difficulties  when  attempting  to  deduce  temperatuie  and 
emissivity  from  measured  radiance.  These  difficulties  are  due  to  the  fast  motion  and  high  rate  of 
temperature  change  of  the  objects.  Temperature  is  typically  changing  on  the  order  of  hundieds  of 
degrees  per  millisecond,  while  motion  occurs  at  typical  velocities  of  1000  feet  per  second.  The 
accuracy  of  the  pyrometer  of  this  dissertation  will  not  suffer  from  the  dynamic  nature  of  the 
energetic  environment  and,  consequently,  is  the  prefect  tool  to  be  used  to  validate  the  empiiical 
calculations  of  the  munition  code  simulations. 

The  next  step  in  this  research  is  science  only  limited  detail  was  given  to  substantiate  the 
effect  of  surface  topography  on  emissivity  a  detailed  analysis  needs  to  be  done.  In  the  literature 
there  is  an  indication  that  diffraction  through  the  surface  may  provide  significant  insight  in 
describing  the  wavelength  dependent  emissivity  for  may  materials.  It  appears  that  incorporating 
this  effect  into  the  overall  model  for  emissivity  is  rich  for  development.  In  addition  to  a  more 
complete  model  for  emissivity  work  needs  to  be  done  to  simplify  the  lens  design  and  find  a  lower 
cost  substitution  for  the  Amber  Inc.  camera.  The  objective  here  is  to  lower  the  cost  of  the  overall 
system  to  compete  with  other  off-the-shelf  single  color  pyrometers.  In  addition,  further  work 
needs  to  be  done  to  characterize  the  emissivity  of  other  materials  of  interest.  After  emissivity  is 
characterized,  the  correct  algorithm  can  be  selected,  and  system  accuracy  analyzed.  When  the 
imaging  spectrometer  of  this  dissertation  is  used  with  an  accurate  emissivity  model  a  highly 
accurate  and  low  cost  optical  pyrometer  will  be  created,  suited  perfectly  for  high-speed  industrial 
applications. 
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Appendix  A-Temperature  Estimation  Software 


The  programs  contained  herein  are  those  written  to  process  the  raw  data  and  output  the 
estimated  temperature.  There  are  basically  two  programs,  a  gain  and  offset  procedure  and  the 
main  data  reduction  procedure.  The  main  data  reduction  procedure  consists  of  the  Gaussian 
reduction,  singular  value  decomposition  and  Levenberg-Marquet  solutions  to  fitting  the  four 
color  data. 

There  are  three  function  calls  that  accomplish  the  inversion  using  LU-decomposition  to 
implement  Gaussian  elimination;  MatbsD%,  MatInvD%  and  MatluD%.  An  LU  matrix  is  an 
upper  triangle  matrix  (one  with  all  non-zero  entries  on  or  above  the  main  diagonal)  and  can  be 
solved  immediately.  The  goal  of  Gaussian  elimination  is  to  transform  a  non  upper  triangle 
system  into  an  equivalent  triangular  one. 

MatbsD%  takes  a  matrix  in  LU  form,  found  by  matluD%,  and  a  vector  b  and  solves  the 
system  Ux=Lb  for  x.  matrices  A,b,x  are  double  precision.  Parameters:  LU  matrix  in  A, 
corresponding  pivot  vectors  in  rpvt  and  cpvt,  right  side  in  b  Returns:  solution  in  x,  b  is  modified, 
rest  unchanged 

MatInvD%  uses  the  matluD%  and  matbsD  procedures  to  invert  a  square,  double  precision 
matrix.  Let  e(N)  contain  all  zeroes  except  for  the  jth  position,  which  is  1.  Then  the  jth  column 
of  A is  x,  where  Ax=e.  Parameters:  A(n  x  n)  matrix  Returns:  A 

MatluD%  does  Gaussian  elimination  with  total  pivoting  to  put  a  square,  double 
precision  matrix  in  LU  form.  The  multipliers  used  in  the  row  operations  to  create  zeroes  below 
the  main  diagonal  are  saved  in  the  zero  spaces.  Parameters:  A(n  x  n)  matrix,  rpvt(n)  and  cpvt(n) 
permutation  vectors  used  to  index  the  row  and  column  pivots  Returns:  A  in  LU  form  with 
corresponding  pivot  vectors;  the  total  number  of  pivots  in  count,  which  is  used  to  find  the  sign  of 
the  determinant. 
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Al.  Gain  And  Offset  Acquired  Data  Program: 


Main  routine 

DECLARE  SUB  gainoff  (a#(),  BBmeslO,  MLbb'O,  MLmeslO,  tlow!,  tmed!,  thigh!,  temp#) 

DECLARE  SUB  calcML  (a#0,  MLbb!(),  MLsamp!(),  temp#0,  templ#()) 

DECLARE  SUB  getbbdata  (FILTERMAX!,  BBmes!(),  MLmes!()) 

DECLARE  SUB  Normalize  (a#0) 

DECLARE  SUB  PLOTDATA  (a#0) 

DIM  SHARED  o  AS  STRING  *  1 
DIM  SHARED  ML#(4,  10) 

COMMON  SHARED  DLAM#,  FILTERMAX,  LAMMAX#(),  LAMMIN#(),  EMISS#,  cl#,  c2#,  c3# 

SCREEN  9 

CLS 

cl#  =3.7415  *  10  A  4 
c2#=  1.43879*  10  A4 
c3#  =  1.88365  *  10  A  23 

DIM  a#(5,  500),  BBmes(4, 10),  MLbb(4,  10),  MLmes(4, 10),  MLsamp(4, 10),  tempi #(10),  temp#(10) 
DIM  LAMMIN#(5),  LAMMAX#(5) 

EMISS#  =  .9 
FILTERMAX  =  4 
DLAM#  =  .005 
LAMMIN#(1)  =  2.7 
LAMMAX#(1)  =  5.2  -  DLAM# 

LAMMIN#(2)  =  2.7 
LAMM  AX#  (2)  =  5.2  -  DLAM# 

LAMM1N#(3)  =  2.7 
LAMMAX#(3)  =  5.2  -  DLAM# 

LAMMIN#(4)  =  2.7 
LAMMAX#(4)  =  5.2  -  DLAM# 

LAMMIN#(5)  =  2.7 
LAMMAX#(5)  =  5.2  -  DLAM# 

CALL  Normalize(a#0) 

CALL  PLOTDATA(a#0) 

CALL  calcML(a#(),  MLbb(),  MLsampQ,  temp#(),  templ#()) 

CALL  getbbdata(FILTERMAX,  BBmesQ,  MLmcsQ) 


tlow  =  1 
tmed  =  3 
thigh  =  6 


CALL  gainoff(a#(),  BBmes(),  MLbbO,  MLmes(),  tlow,  tmed,  thigh,  templ#(tmed)) 
END 

SUB  calcML  (a#(),  MLbb(),  MLsampO,  temp#(),  templ#()) 

temp#(l)  =  98.6 +  273 
temp#(2)=  108.4  +  273 
temp#(3)=  117.8  +  273 
temp#(4)  =  129  +  273 
temp#(5)=  137.9  +  273 
temp#(6)=  149.4  +  273 
temp#(7)=  149.4  +  273 
temp#(8)=  1156  +  273 
temp#(9)  =  168.7  +  273 
temp#(10)  =  169.3  +  273 

FORt=l  TO  10 
FOR  filter  =  1  TO  FILTERMAX 
MLOLD#  =  0 
1  =  0 

FOR  lam#  =  LAMMIN#(filter)  TO  LAMMAX#(filter)  STEP  DLAM# 

1  =  1+1 

R#  =  a#  (filter,  1) 
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MLbb#  =  cl#  /  (lam#  A  5  *  (EXP(c2#  /  (lam#  *  temp#(t)))  - 1))  *  EMISS# 
MLbb(filter,  t)  -  MLbb#  *  R#  *  DLAM#  +  MLOLD# 

MLOLD#  =  MLbb(filter,  t) 

NEXT  lam# 

’PRINT  MLhigh#(filter) 

NEXT  filter 
NEXT  t 

templ#(l)  =  108  +  273 
templ#(2)  =119  +  273 
templ#(3)  =  140  +  273 
templ#(4)=  140  +  273 
templ#(5)=  159.5  +  273 
templ#(6)  =  159.6  +  273 
templ#(7)  =  164.9  +  273 
templ#(8)  =  168.8  +  273 
templ#(9)  =  179.7  +  273 
templ#(10)  —  207.1  +  273 

FOR  t  =  1  TO  10 
FOR  filter  =  1  TO  FILTERMAX 
MLOLD#  =  0 
1  =  0 

FOR  lam#  =  LAMMIN#(filter)  TO  LAMMAX# (filter)  STEP  DLAM# 

1  =  1+1 

R#  =  a# (filter,  I) 

MLbb#  =  cl#  /  (lam#  A  5  *  (EXP(c2#  /  (lam#  *  templ#(t)))  -  I))  *  1 
MLsamp(filter,  t)  =  MLbb#  *  R#  *  DLAM#  +  MLOLD# 

MLOLD#  =  MLsamp(fi!ter,  t) 

NEXT  lam# 

’PRINT  MLhigh#  (filter) 

NEXT  filter 
NEXT  t 

'FOR  filter  =  1  TO  4:  PRINT  MLh(filter,  1);  MLh(filter,  2);  MLh(filter,  3):  NEXT  filter 
END  SUB 


SUB  gainoff  (a#0>  BBmes0,  MLbbO?  MLmes(),  tlow,  tmed,  thigh,  temp#) 

DIM  MLtemp#(4) 

DIM  dfilter(4) 

FOR  filter  =  1  TO  FILTERMAX  'Gain  and  offset  the  data  using  low  and  high  blackbody  values 

ofcet  =  (MLbb(filter,  tlow)  *  BBmes(filter,  thigh)  -  MLbb(filter,  thigh)  *  BBmes(filter,  tlow))  /  (MLbb(filter,  tlow)  -  MLbb(filter,  thigh)) 
gall  =  (BBmes(filter,  thigh)  -  ofcet)  /  MLbb  (filter,  thigh) 

IF  filter  =  1  THEN  gall  =  gall  *  .93 
'IF  filter  =  4  THEN  gall  =  gall  *  .9 
FOR  t=  1  TO  10 

MLmes(filter,  t)  =  (MLmes(filter,  t)  -  ofcet)  /  gall 
'print  MLmes(filter,  t) 

NEXT  t 

BBmes(filter,  tmed)  =  (BBmes(filter,  tmed)  -  ofcet)  /  gall'Gain  and  offset  mid  bb  reading 
PRINT  ML;  ML#(4,  10),  gall;  ofcet 
NEXT  filter 

FOR  t  =  temp#  -  50  TO  temp#  +  50  STEP  5 
FOR  filter  =  1  TO  FILTERMAX 
MLOLD#  =  0 
1  =  0 

FOR  lam#  =  LAMMIN# (filter)  TO  LAMMAX# (filter)  STEP  DLAM# 

1  =  1+1 

R#  =  a#  (filter,  I) 

MLbb#  =  cl#  /  (lam#  A  5  *  (EXP(c2#  /  (lam#  *  t))  -  1))  *  1 
MLtemp#(filter)  =  MLbb#  *  R#  *  DLAM#  +  MLOLD# 
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MLOLD#  =  (MLtemp#(filter)) 

NEXT  lam# 

'PRINT  MLmed#(fiIter);  med(filter) 

NEXT  filter 

FOR  I  =  1  TO  FILTERMAX 

dfilter(I)  =  MLmes(l,  tmed)  /  MLtemp#(I)’  -  BBmes(i,  tmed) 

NEXT  I 

IF  t  =  temp#  THEN  COLOR  5  ELSE  COLOR  1 0 
PRINT  t;  dfilter(l);  dfilter(2);  dfilter(l)  -  dfilter(2);  dfilter(3);  dfilter(4) 

NEXT  t 

OPEN  "a:\DATA.FIL"  FOR  OUTPUT  AS#  1 

FOR  filter  =  1  TO  FILTERMAX  'Gain  and  offset  the  data  using  low  and  high  blackbody  values 

1  =  0 

FOR  lam#  =  LAMMIN#(filter)  TO  LAMMAX#(filter)  STEP  DLAM# 

1  =  1+1 

PRINT  #1,  a#(filter,  I) 

’  PRINT  a#(filter,  i) 

NEXT  lam# 

NEXT  filter 
CLOSE  #1 

OPEN  "a:\DATA.SYS”  FOR  OUTPUT  AS  #2 
PRINT  #2,  FILTERMAX 
PRINT  #2,  DLAM# 

FOR  I  =  1  TO  FILTERMAX:  PRINT  #2,  LAMMIN#(I):  NEXT  I 
FOR  I  =  1  TO  FILTERMAX:  PRINT  #2,  LAMMAX#(I):  NEXT  I 

FOR  t  =  1  TO  10 

FOR  I  =  1  TO  FILTERMAX:  PRINT  #2,  MLmes(I,  t):  NEXT  I 
NEXT  t 

CLOSE  #2 
END  SUB 

SUB  getbbdata  (FILTERMAX,  BBmes(),  MLmes()) 

OPEN  "a:\blackbol.arr"  FOR  BINARY  AS  #1 
GET  #1,  1,0 
a!  -  ASC(o) 

GET  #1,  a!,  o 
PRINT 

FOR  1=  1  TO  FILTERMAX 

FOR  t=  1  TO  10 
GET  #1, ,  BBmes(I,  t) 

1  PRINT  BBmes(I,  t); :  INPUT  a$ 

NEXT  t 
1  INPUT  aS 
NEXT  I 

CLOSE  #1 

OPEN  "a:\sample2. arr"  FOR  BINARY  AS  #1 
GET  #l,l,o 
a!  =  ASC(o) 

GET  #1,  a!,  o 
PRINT 

FOR  I  -  1  TO  FILTERMAX 

FOR  t=  1  TO  10 
GET  #1, ,  MLmes(I,  t) 

'  PRINT  MLmes(i,  t),  i,  t; :  INPUT  a$ 

NEXT  t 
1  INPUT  aS 
NEXT  I 


106 


CLOSE#! 


END  SUB 


SUB  Normalize  (a#Q) 


OPEN  "a:\filterl  .arrh  FOR  BINARY  AS  #2 
OPEN  "a:\temp.arr'5  FOR  BINARY  AS  #1 


GET  #2,  l,o 
b!  =  ASC(o) 
GET  #2,  b!,  o 


FOR  I  -  1  TO  499 
GET  #1, ,  m# 
a#(5, 1)  =  m#  *  1.5 
NEXT  I 

max  =  0 

FOR  filter  =  1  TO  4 
GET  #2, ,  0F1 
FOR  I  =  1  TO  499 
GET  #2, ,  Ml 
'  a#(filter,  I)  =  MI 

IF  a# (5, 1)  >  .01  AND  I  <  480  AND  I  >  10  THEN  a#(filter,  I)  =  Ml  /  a#(5, 1) 
IF  max  <  a#(filter,  I)  THEN  max  =  a#(filter,  I) 

NEXT  I 
NEXT  filter 

CLOSE  #1 
CLOSE  #2 

END  SUB 


SUB  PLOTDATA  (a#0) 

LINE  (5,  l)-(5,  100),  2 
LINE  (1,  100)-(400,  100),  2 

FOR  L  -  10  TO  100  STEP  10:  LINE  (1,  LH10,  L),  2:  NEXT  L 

FOR  filter  =  1  TO  FILTERMAX 
1  =  0 

PSET  (5,  100) 

FOR  lam#  -  LAMMIN#  (filter)  TO  LAMM  AX#  (filter)  STEP  DLAM#  'Calculate  and  store  emissivities  Plot  filter  values 
1  =  1+1 

X  =  (lam#  -2.7)  *200 
Y  =  (1  -  1  *  a#(filter,  I))  *  100 
LINE  -(X,  Y),  filter 

1  LOCATE  20,  1 :  PRINT  lam# ; :  INPUT  a$ 

NEXT  lam# 

*  LOCATE  20,  1:  INPUT  a$ 

NEXT  filter 


END  SUB 
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A2.  Process  Data  Program: 


Main  Subrotine 

DECLARE  SUB  FITDATUM  (TEMP#,  A#()) 

DECLARE  SUB  GETDATA  (samp!) 

DECLARE  FUNCTION  MATBSD%  (AO  AS  DOUBLE,  B()  AS  DOUBLE,  X()  AS  DOUBLE) 

DECLARE  FUNCTION  MATINVD%  (A()  AS  DOUBLE) 

DECLARE  FUNCTION  MATLUD%  (AO  AS  DOUBLE) 

DECLARE  SUB  SVBKSB  (U!0,  W!(),  V!0,  M!,N!,  MP!,  np!,  B!(),  X!()) 

DECLARE  SUB  SVDCMP  (A!0,  M!,  N!,  MP!,  np!,  W!(),  V!()) 

DECLARE  SUB  MRQCOF  (X!(),  Y!0,  SIG!(),  NDATA!,  A!(),  MA!,  LISTA!(),  MFIT!,  ALPHA!(),  BETA!(),  NALPL  CHISQ!,  DUM!) 
DECLARE  SUB  GAUSSJ  (A!(),  N!,  NP!,  B!(),  M!,  MP!) 

DECLARE  SUB  COVSRT  (COVAR!(),  NCVM!,  MA!,  LISTA!(),  MFIT!) 

COMMON  SHARED  BETAO 

COMMON  SHARED  cl#,  c2#,  C3#,  DLAM#,  filtermax,  LAMMAX#(),  lammin#() 

COMMON  SHARED  ML#0,  R#0,  filenames 

DIM  SHARED  LO  AS  INTEGER,  UP  AS  INTEGER 

DIM  SHARED  CONTINUE  AS  INTEGER,  COUNT  AS  INTEGER 

DIM  SHARED  RPVT(LO  TO  UP)  AS  INTEGER,  CPVT(LO  TO  UP)  AS  INTEGER 

DIM  LAMMAX#(10),  !ammin#(10),  ML#(10),  MQ#(10),  ATARGET#(10) 

DIM  R#(4,  500) 


'Decluration  of  constants 

SCREEN  9 

CLS 

cl#  =  3.7415  *  10  A4 
c2#=  1.43879*  10  A  4 
C3#  =  1.88365  *  10  A  23 

'CALL  ACAL(ATARGET#0) 

'CALL  FILTERS 

'CALL  FINDDATA(500,  ATARGET#(),  ML#(),  MQ#()) 


filenames  =  ”a:DATA" 
samp  =  7 

CALL  GETDAT  A(samp) 

DIM  A#(filtermax  -  1) 

FOR  TEMP#  =  350  TO  450  STEP  5  'search  temperature  to  find  solution 
CALL  FITDATUM  (TEMP#,  A #()) 

NEXT  TEMP# 


SUB  ACAL  (ATARGET#0) 

'Do  curvefit  of  data  points  as  selected 

DLAM#  =  .001 
filtermax  =  4 
lammin#(l)  =  3.4 
lammin#(2)  =  3.5 
lammin#(3)  =  3.6 
lammin#(4)  =  3.7 
lammin#(5)  =  3.8 
lammin#(6)  =  3.9 
lammin#(7)  =  4! 
lammin#(8)  =  4.1 
lammin#(9)  =  4.2 
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Iammin#(10)  =  4.3 
LAMMAX#(1)  -  3.5 
LAMM  AX#  (2)  =  3.6 
LAMMAX#(3)  =  3.7 
LAMM  AX# (4)  -  3.8 
LAMM  AX#  (5)  =  3.9 
LAMMAX#(6)  —  4! 
LAMMAX#(7)  =  4.1 
LAMMAX#(8)  =  4.2 
LAMM  AX#  (9)  =  4.3 
LAMMAX#(1 0)  =  4.4 


ATARGET#(1)  =  0 
ATARGET#(2)  -  0 
ATARGET#(3)  =  0 
ATARGET#(4)  =  -.1 
ATARGET#(5)  -  0 
ATARGET#(6)  =  2 
ATARGET#(7)  -  0 
ATARGET#(8)  =  0 
ATARGET#(9)  -  0 
ATARGET#(10)  -  0 

END  SUB 


SUB  FILTERS 


’Generate  the  filter  values.  The  filter  shapes  used  her  are  gaussian. 

OPEN  "DATA.FIL"  FOR  OUTPUT  AS  #1 

FOR  filter  =  1  TO  filtermax 
1  =  0 

FOR  lam#  =  lam  min#  (filter)  TO  LAMMAX#  (filter)  STEP  DLAM# 

1  =  1+1 

AVG  =  (lammin#(filter)  +  LAMMAX#  (filter))  /  2 

WID  =  (I ammin# (filter)  -  LAMMAX# (filter))  /  2 

R#  =  EXP(-3.141  *  ((lam#  -  AVG)  /  WID)  A  2)  'Gaussian  filter  shape 

PRINT  #1,R# 

'  PRINT  R# 

NEXT  lam# 

NEXT  filter 

CLOSE #1 

END  SUB 


SUB  FINDDATA  (TEMP#,  ATARGET#Q,  ML#(),  MQ#()) 


'Calculate  the  radiant  extdance  in  terms  of  pawer  (ML#)  and  photons 
'per  second  (MQ#) 

OPEN  "DATA.SYS"  FOR  OUTPUT  AS  #2 
OPEN  "DATA.FIL"  FOR  INPUT  AS  #  I 
FOR  filter  =1  TO  filtermax 
MLOLD#  =  0 
MQOLD#  =  0 
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FOR  lam#  =  lammin#(filter)  TO  LAMMAX# (filter)  STEP  DLAM# 

INPUT  #1,R# 

AOLD# = 0 

FOR  I  =  1  TO  filtermax  -  1 

emiss#  =  ATARGET#(I)  *  lam#  A  (filtermax  -  1  - 1)  +  AOLD# 

AOLD#  =  emiss# 

NEXT  I 

MLbb#  =  cl#  /  (lam#  A  5  *  (EXP(c2#  /  (lam#  *  TEMP#))  - 1))  *  emiss# 
MQBB#  =  C3#  /  (lam#  A  4  *  (EXP(c2#  /  (lam#  *  TEMP#))  -  1))  *  emiss# 
ML#(filter)  =  MLbb#  *  R#  *  DLAM#  +  MLOLD# 

MQ#(filter)  =  MQBB#  *  R#  *  DLAM#  +  MQOLD# 

MLOLD#  =  ML#(fiIter) 

MQOLD#  =  MQ#(fiIter) 

NEXT  lam# 

NEXT  filter 


PRINT  #2,  filtermax 
PRINT  #2,  DLAM# 

FOR  I  =  1  TO  filtermax:  PRINT  #2,  lammin#(I):  NEXT  I 
FOR  I  =  1  TO  filtermax:  PRINT  #2,  LAMMAX#(I):  NEXT  I 
FOR  1  =  1  TO  filtermax:  PRINT  #2,  ML#(I):  NEXT  I 

CLOSE #1 
CLOSE  #2 
END  SUB 


SUB  FITDATUM  (TEMP#,  A#0) 

'Find  coef.  at  a  given  temperature  and  compute  temperature 

DIM  MATRIXsvd#(l  TO  filtermax,  1  TO  filtermax  - 1) 

DIM  MATRIX#(I  TO  filtermax  -  1,  1  TO  filtermax  -  1) 

DIM  MLOLD#(l  TO  filtermax) 

DIM  U(4, 3),  W(4),  V(4, 3),  X(4),  ML(4),  ER(4) 

'OPEN  filenames  +  ".FIL”  FOR  INPUT  AS  #  1 

’  FOR  filter  =  1  TO  filtermax  - 1  ’Generate  the  matrix  of  simotanious  equations 
'  FOR  I  =  1  TO  filtermax  - 1 :  MLOLD#(I)  =  0:  NEXT  I 
'  FOR  lam#  =  lammin#(filter)  TO  LAMMAX#(filter)  STEP  DLAM# 

1  INPUT  #1,R# 

•  FOR  COL  =  1  TO  filtermax  - 1 

'  MLbb#  =  c  1  #  /  (lam#  A  (5  -  (filtermax  - 1  -  COL))  *  (EXP(c2#  /  (lam#  *  TEMP#))  - 1 )) 

•  MATRIX#(COL,  filter)  =  MLbb#  *  R#  *  DLAM#  +  MLOLD#(COL) 

•  MLOLD#  (COL)  =  MATRIX#(COL,  filter) 

’  NEXT  COL 

1  NEXT  lam# 

1  NEXT  filter 

'  ERRCODE%  =  MATTNVD%(MATRIX#0) 

'FOR  ROW  =  1  TO  filtermax  - 1  '  determine  emissivity  coef. 

'  AOLD#  =  0 

'  FOR  COL  =  1  TO  filtermax  - 1 

'  a# (ROW)  =  MATRIX#(COL,  ROW)  *  ML#(COL)  +  AOLD# 

•  AOLD#  =  a#  (ROW) 

'  NEXT  COL 
■NEXT  ROW 

'FOR  I  =  1  TO  FILTERMAX  - 1 :  a  =  a#(I):  PRINT  a; :  NEXT  I 

'  MLOLD#  =  0 
1  MLOLD  1#  =  0 

'  FOR  lam#  =  lammin#(filtermax)  TO  LAMMAX#(filtermax)  STEP  DLAM# 

'  INPUT  #1,R# 
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AOLD# = 0 


FOR  I  -  1  TO  filtermax  -  1  'determine  emissivity  of  forth  filter 
1  emiss#  =  a# (I)  *  lam#  A  (filtermax  -  1  - 1)  +  AOLD# 

'  AOLD#  =  emiss# 

'  NEXT I 

'  MLbb#  =  cl#  /  (lam#  A  5  *  (EXP(c2#  /  (lam#  *  TEMP#))  -  1))  *  emiss# 

1  ML#  =  MLbb#  *  R#  *  DLAM#  +  MLOLD# 

'  MLOLD#  =  ML# 

’  NEXT  lam# 

'  er  =  ML#  -  ML# (filtermax) 

1  AOLD#  =  0 

'  FOR  I  =  1  TO  filtermax  -  1 

1  em#  =  a#(l)  *  (lammin#(4)  +  450  *  DLAM#)  A  (filtermax  -  1  - 1)  +  AOLD# 

1  AOLD#  =  em# 

'  NEXT  I 

'  emi s  =  ML# (filtermax)/ ML  1# 

PRINT  emis;  emiss#;  TEMP#  -  273;  emis  -  em#,  emis  -  emiss# 

'CLOSE  #1 

FOR  filter  =  1  TO  filtermax  'Generate  the  matrix  for  svd 

L  =  0 

FOR  I  =  1  TO  filtermax  -  1:  MLOLD#(I)  =  0:  NEXT  I 

FOR  lam#  =  Iammin#(filter)  TO  LAMMAX# (filter)  STEP  DLAM# 

L  =  L  +  1 

R#  =  R#(filter,  L) 

FOR  COL  =  1  TO  filtermax  -  1 

MLbb#  =  cl#  /  (lam#  A  (5  -  (filtermax  -  1  -  COL))  *  (EXP(c2#  /  (lam#  *  TEMP#))  -  I)) 
M ATRIXsvd# (fi Iter,  COL)  =  MLbb#  *  R#  *  DLAM#  +  MLOLD#(COL) 

MLOLD# (COL)  =  MATRIXsvd#(filter,  COL) 

NEXT  COL 
NEXT  lam# 

NEXT  filter 

FOR  filter  =  1  TO  filtermax  'Copy  matrixsvd#  into  U# 

FOR  COL  =  1  TO  filtermax  -  I 
U(filter,  COL)  =  MATRIXsvd#(fiIter,  COL) 

NEXT  COL 
NEXT  filter 


CALL  SVDCMP(U(),  4,  3,  np,  np,  W(),  VO)  'Decompose  matrix  A 

FOR  filter  —  1  TO  filtermax 
ML(filter)  -  ML#(filter) 

NEXT  filter 


WMAX  =  0! 

FOR  filter  —  1  TO  filtermax  ’Find  maximum  singular  value 

IF  W(filter)  >  WMAX  THEN  WMAX  =  W(filter) 

NEXT  filter 

WMIN  =  WMAX  *  (.000001)  'Define  "small’’ 

FOR  filter  =  1  TO  filtermax  'Zero  the  "small"  singular  values 

IF  W(filter)  <  WMIN  THEN  W(filter)  =  0! 

NEXT  filter 

'Backsubstitute  for  each  right-hand  side  vector 
CALL  SVBKSB(U(),  W0,  V0,  4,  3,  np,  np,  ML(),  X()) 

FOR  filter  =  1  TO  filtermax 
MLOLD#  -  0 
MLOLD  1#  =  0 
L  =  0 

FOR  lam#  -  lammin#(filter)  TO  LAMMAX# (filter)  STEP  DLAM# 

L  =  L  +  1 
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R#  =  R#(filter,  L) 

AOLD#  =  0 

FOR  I  =  1  TO  filtermax  - 1  'determine  emissivity  of  forth  filter 
emiss#  =  X(I)  *  lam#  A  (filtermax  - 1  - 1)  +  AOLD# 

AOLD#  =  emiss# 

NEXT  I 

MLbb#  =  cl#  /  (lam#  A  5  *  (EXP(c2#  /  (lam#  *  TEMP#))  -  1))  *  emiss# 

ML#  =  MLbb#  *  R#  *  DLAM#  +  MLOLD# 

MLOLD#  =  ML# 

NEXT  lam# 

ER(filter)  =  ML#  -  ML#(fiIter) 

NEXT  filter 

PRINT 

PRINT  TEMP#  -  273;  ER(1)  +  ER(2)  +  ER(3)  +  ER(4); 

FOR  K  =  1  TO  filtermax  -  1 :  PRINT  USING  "#####.######";  X(K); :  NEXT  K 

END  SUB 


SUB  GETDATA  (samp) 

DIM  MLt(4,  10) 

OPEN  filenames  +  ".SYS”  FOR  INPUT  AS  #2 

INPUT  #2,  filtermax 
INPUT  #2,  DLAM# 

FOR  I  -  1  TO  filtermax:  INPUT  #2,  lammin#(I):  NEXT  I 
FOR  I  =  1  TO  filtermax:  INPUT  #2,  LAMM  AX#  (I):  NEXT  I 


t=l 

FOR  t  =  1  TO  10 
FOR  I  =  1  TO  filtermax 
INPUT  #2,  MLt(I,  t) 

’PRINT  MLt(I,t); 

NEXT  I 
PRINT 
NEXT  t 

FOR  I  ~  1  TO  filtermax 
ML# (I)  =  MLt(I,  samp) 

'  PRINT  MLt(I,  samp) 

NEXT  I 

CLOSE #2 

OPEN  filenames  +  ".FIL"  FOR  INPUT  AS  #1 

FOR  filter  -  1  TO  filtermax  ’Generate  the  matrix  of  filter  values 


L  =  0 

FOR  lam#  -  lammin#(filter)  TO  LAMM  AX#  (filter)  STEP  DLAM# 
L  =  L  +  1 
INPUT  #1,R# 

R# (filter,  L)  =  R# 

’PRINT  R# (FILTER,  L);  FILTER;  L'; :  INPUT  AS 
NEXT  lam# 

NEXT  filter 
CLOSE #1 


END  SUB 
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FUNCTION  MATBSD%  (AO  AS  DOUBLE,  B()  AS  DOUBLE,  X()  AS  DOUBLE) 

ON  LOCAL  ERROR  GOTO  DBSERR:  MATBSD%  =  0 
'do  row  operations  on  b  using  the  multipliers  in  L  to  find  Lb 
FOR  PVT%  =  LO  TO  (UP  - 1) 

C%  =  CPVT(PVT%) 

FOR  ROW%  =  (PVT%  +  1)  TO  UP 
R%  =  RPVT(ROW%) 

B(R%)  =  B(R%)  +  A(R%,  C%)  *  B(RPVT(PVT%)) 

NEXT  ROW% 

NEXT PVT% 

FOR  ROW%  =  UP  TO  LO  STEP  -1  'backsolve  UX-Lb  to  find  x 
C%  =  CPVT(ROW%) 

R%  =  RPVT(ROW%) 

X(C%)  =  B(R%) 

FOR  COL%  =  (ROW%  +  1)  TO  UP 
X(C%)  -  X(C%)  -  A(R%,  CPVT(COL%))  *  X(CPVT(COL%)) 

NEXT  COL% 

X(C%)  =  X(C%)  /  A(R%,  C%) 

NEXT  ROW% 

DBSEXIT: 

EXIT  FUNCTION 
DBSERR: 

MATBSD%  =  ERR 
RESUME  DBSEXIT 
END  FUNCTION 


SUB  MRQMIN  (X(),  YO,  SIGO,  NDATA,  A(),  MA,  LISTA(),  MFIT,  COVAR(),  ALPH A(), 
NCA,  CHISQ,  DUM,  ALAMDA) 

STATIC 

DIM  ATRY(MA),  DA(MA,  1) 

IF  ALAMDA  <0!  THEN 
KK  —  MFIT  +  1 
FOR  J  —  1  TO  MA 
IHIT  =  0 

FOR  K  -  1  TO  MFIT 
IF  LISTA(K)  =  J  THEN  IHIT  =  IHIT  +  1 
NEXT  K 

IF  IHIT  =  0  THEN 
LISTA(KK)  =  J 
KK  =  KK  +  1 
ELSEIF  IHIT  >  1  THEN 
PRINT  "Improper  permutation  in  LISTA" 

EXIT  SUB 
END  IF 
NEXT  J 

IF  KK  o  MA  +  1  THEN  PRINT  "Improper  permutation  in  LISTA":  EXIT  SUB 
ALAMDA  =.001 

CALL  MRQCOF(X0,  Y(),  SIGO,  NDATA,  AO,  MA,  LISTA(),  MFIT,  ALPHA(),  BETA(),  NCA,  CHISQ,  DUM) 

OCHISQ  =  CHISQ 
FOR  J  =  1  TO  MA 
ATRY(J)  =  A(J) 

NEXT  J 
END  IF 

FOR  J  =  1  TO  MFIT 
FORK=l  TO  MFIT 
COVAR(J,  K)  =  ALPHA(J,  K) 

NEXT  K 

COVAR(J,  J)  =  ALPHA(J,  J)  *  (1 !  +  ALAMDA) 

DA(J,  1)  =  BETA(J,  1) 
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NEXT  J 

CALL  GAUSS J(COVAR(),  MFIT,  NCA,  DAO,  1,1) 

IF  ALAMDA  =  0!  THEN 

CALL  COVSRT(COVAR(),  NCA,  MA,  LISTAO,  MFIT) 

ERASE  DA,  ATRY 
EXIT  SUB 
END  IF 

FOR  J  =  1  TO  MFIT 

ATRY(LISTA(J))  -  A(LISTA(J))  +  DA(J,  1) 

NEXT  J 

CALL  MRQCOF(X(),  YO,  SIGO,  NDATA,  ATRY(),  MA,  LISTAO,  MFIT,  COVAR(),  DA(),  NCA,  CHISQ,  DUM) 
IF  CHISQ  <OCHISQ  THEN 
ALAMDA  =  .  1  *  ALAMDA 
OCHISQ  =  CHISQ 
FOR  J  =  1  TO  MFIT 
FOR  K  =  1  TO  MFIT 
ALPHA(J,  K)  =  COVAR(J,  K) 

NEXT  K 

BETA(J,  I)  “  DA(J,  1) 

A(LISTA(J))  =  ATRY  (LISTA(J)) 

NEXT  J 
ELSE 

ALAMDA  =  10!  *  ALAMDA 
CHISQ  =  OCHISQ 
END  IF 

ERASE  DA,  ATRY 
END  SUB 


FUNCTION  MATINVD%  (AO  AS  DOUBLE) 

ON  LOCAL  ERROR  GOTO  DINVERR:  ERRCODE%  =  0 
LO  -  LBOUND(A,  1) 

UP  =  UBOUND(A,  1) 

DIM  AIN(LO  TO  UP,  LO  TO  UP)  AS  DOUBLE 

DIM  E(LO  TO  UP)  AS  DOUBLE,  X(LO  TO  UP)  AS  DOUBLE 

REDIM  RPVT(LO  TO  UP)  AS  INTEGER,  CPVT(LO  TO  UP)  AS  INTEGER 

ERRCODE%  =  MATLUD%(A0)  ’Get  LU  matrix 

IF  NOT  CONTINUE  THEN  ERROR  ERRCODE% 

FOR  COL%  -  LO  TO  UP  'Find  AM  one  column  at  a  time 

E(COL%)  =  1# 

BSERRCODE%  =  MATBSD%(A(),  E(),  X()) 

IF  BSERRCODE%  THEN  ERROR  BSERRCODE% 

FOR  ROW%  =  LO  TO  UP 
AIN(ROW%,  COL%)  =  X(ROW%) 

E(ROW%)  =  0# 

NEXT  ROW% 

NEXT  COL% 

FOR  COL%  =  LO  TO  UP  'Put  AA-1  in  A 

FOR  ROW%  =  LO  TO  UP 
A(ROW%,  COL%)  =  AIN(ROW%,  COL%) 

NEXT  ROW% 

NEXT  COL% 

IF  ERRCODE%  THEN  ERROR  ERRCODE% 

DINVEXIT: 

ERASE  E,  X,  AIN,  RPVT,  CPVT 
MATINVD%  =  ERRCODE% 

EXIT  FUNCTION 
DINVERR: 

ERRCODE%  =  (ERR  +  5)  MOD  200  -  5 
RESUME  DINVEXIT 
END  FUNCTION 
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FUNCTION  MATLUD%  (AO  AS  DOUBLE) 

ON  LOCAL  ERROR  GOTO  DLUERR:  ERRCODE%  -  0 
'Checks  if  A  is  square,  returns  error  code  if  not 

IF  NOT  (LO  =  LBOUND(A,  2)  AND  UP  *  UBOUND(A,  2))  THEN  ERROR  1 98 
DIM  ROWNORM(LO  TO  UP)  AS  DOUBLE 
COUNT  =  0  'initialize  count,  continue 

CONTINUE  =  -l 

FOR  ROW%  =  LO  TO  UP  'initialize  rpvt  and  cpvt 

RPVT(ROW%)  =  ROW% 

CPVT(ROW%)  =  ROW% 

ROWNORM(ROW%)  -  0#  'find  the  row  norms  of  A() 

FOR  COL%  =  LO  TO  UP 

ROWNORM(ROW%)  =  ROWNORM(ROW%)  +  ABS(A(ROW%  COL%)) 

NEXT  COL% 

IF  ROWNORM(ROW%)  =  0#  THEN  'if  any  rownorm  is  zero,  the  matrix 

CONTINUE  =  0  'is  singular,  set  error,  exit  and 

ERROR  1 99  'do  not  continue 

END  IF 
NEXT  ROW% 

FOR  PVT%  =  LO  TO  (UP  - 1) 

'Find  best  available  pivot 

MAX#  =  0#  'checks  all  values  in  rows  and  columns  not 

FOR  ROW%  =  PVT%  TO  UP  'already  used  for  pivoting  and  saves  the 
R%  =  RPVT(ROW%)  'largest  absolute  number  and  its  position 

FOR  COL%  -  PVT%  TO  UP 
C%  =  CPVT(COL%) 

TEMP#  =  ABS(A(R%,  C%))  /  ROWNORM(R%) 

IF  TEMP#  >  MAX#  THEN 
MAX#  =  TEMP# 

BESTROW%  =  ROW%  'save  the  position  of  new  max# 

BESTCOL%  =  COL% 

END  IF 
NEXT  COL% 

NEXT  ROW% 

IF  MAX#  =  0#  THEN  'if  no  nonzero  number  is  found,  A  is 

CONTINUE  =  0  'singular,  send  back  error,  do  not  continue 

ERROR  199 

ELSEIF  PVT%  >  1  THEN  'check  if  drop  in  pivots  is  too  much 

IF  MAX#  <  (deps#  *  oldmax#)  THEN  ERRCODE%  =  199 
END  IF 

oldmax#  =  MAX# 

IF  RPVT(PVT%)  o  RPVT (BESTROW%)  THEN 
COUNT  —  COUNT  +1  'if  a  row  or  column  pivot  is 

SWAP  RPVT(PVT%),  RPVT(BESTROW%)  'necessary,  count  it  and  permute 
END  IF  ’rpvt  or  cpvt.  Note:  the  rows  and 

IF  CPVT(PVT%)  o  CPVT(BESTCOL%)  THEN  'columns  are  not  actually  switched, 

COUNT  =  COUNT  +  1  'only  the  order  in  which  they  are 

SWAP  CPVT(PVT%),  CPVT (BESTCOL%)  'used. 

END  IF 

’Eliminate  all  values  below  the  pivot 
RP%  =  RPVT(PVT%) 

CP%  =  CPVT(PVT%) 

FOR  ROW%  =  (PVT%  +  1)  TO  UP 
R%  =  RPVT  (ROW%) 

A(R%,  CP%)  =  -A(R%,  CP%)  /  A(RP%,  CP%)  ’save  multipliers 
FOR  COL%  =  (PVT%  +  1)  TO  UP 

C%  =  CPVT(COL%)  ’complete  row  operations 

A(R%,  C%)  =  A(R%,  C%)  +  A(R%,  CP%)  *  A(RP%,  C%) 

NEXT  COL% 

NEXT  ROW% 

NEXT  PVT% 

IF  A(RPVT(UP),  CPVT(UP))  =  0#  THEN 
CONTINUE  =  0  'if  last  pivot  is  zero  or  pivot  drop  is 

ERROR  199  'too  large,  A  is  singular,  send  back  error 

ELSEIF  (ABS(A(RPVT(UP),  CPVT(UP)))  /  ROWNORM(RPVT(UP)))  <  (deps#  *  oldmax#)  THEN 
ERRCODE%  =  199  ’if  pivot  is  not  identically  zero  then 
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END  IF  ’continue  remains  TRUE 

IF  ERRCODE%  THEN  ERROR  ERRCODE% 
DLUEXIT: 

MATLUD%  =  ERRCODE% 

EXIT  FUNCTION 
DLUERR: 

IF  ERRCODE%  <  199  THEN  CONTINUE  =  0 
ERRCODE%  =  ERR 
RESUME  DLUEXIT 
END  FUNCTION 


SUB  PLOTDATA  (ATARGET#0) 


1  =  0 

LINE  (5,  1 )-(5,  100),  2 
LINE  (1,1 00)-(400,  100),  2 

FOR  L  =  10  TO  100  STEP  10:  LINE  (1,  L)-(10,  L),  2:  NEXT  L 
FOR  lam#  =  3  TO  5  STEP  DLAM#  'Draw/plot  emissivity  curve 

AOLD#  =  0 

FOR  I  =  1  TO  filtermax  - 1 

emiss#  =  ATARGET#(I)  *  lam#  A  (filtermax  -  1  - 1)  +  AOLD# 

AOLD#  =  emiss# 

NEXT  I 

X  =  (lam#  -3)  *300 
Y  =  (1  -  emiss#)  *  100 
PSET  (X,  Y),  3 
NEXT  lam# 


OPEN  filenames  +  ".FIL"  FOR  INPUT  AS  #  1 
’OPEN  "DATA”  FOR  INPUT  AS  #  1 
FOR  filter  =  I  TO  filtermax 

FOR  lam#  =  lamm in# (filter)  TO  LAMMAX# (filter)  STEP  DLAM#  'Calculate  and  store  emissivities  Plot  filter  values 
INPUT  #1,R# 

'  PRINT  R# 

AOLD#  =  0 

FOR  I  =  1  TO  filtermax  -  1 

emiss#  =  ATARGET#®  *  lam#  A  (filtermax  -  1  -  I)  +  AOLD# 

AOLD#  =  emiss# 

NEXT  I 
*  emiss#  =  1 

X  =  (lam#  -  3)  *  300 
Y  =  (1  -  emiss#  *  R#)  *  100 
PSET  (X,  Y) 

NEXT  lam# 

NEXT  filter 

CLOSE  #1 

END  SUB 


SUB  SVBKSB  (U0,  W0,  V0,  M,  N,  MP,  np,  B0,  XO) 

DIM  TMP(N) 

FOR  J  =  1  TO  N 
S  =  0! 

IF  W(J)  <>  0!  THEN 
FOR  I  =  1  TO  M 
S  =  S  +  U(I,  J)  *  B(I) 

NEXT  I 
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S  =  S/W(J) 

END  IF 
TMP(J)  =  S 
NEXT  J 

FOR  J  =  1  TON 
S  =  0! 

FOR  JJ  =  1  TON 
S  =  S  +  V(J,  JJ)  *  TMP(JJ) 
NEXT  JJ 
X(J)  =  S 
NEXT  J 
ERASE  TMP 
END  SUB 


SUB  SVDCMP  (AO,  M,  N,  MP,  iip,  WO,  VO) 

DIM  RV1(N) 

IF  M  <  N  THEN  PRINT  ’’You  must  augment  A  with  extra  zero  rows  EXIT  SUB 
G  =  0! 

SCALE  =  0! 

ANORM  =  0! 

FORI  =  I  TON 
L  =  I  +  1 

RV1(I)  =  SCALE  *  G 
G  =  0! 

S  =  0! 

SCALE  =  0! 

IF  I  <=  M  THEN 
FOR  K  =  I  TO  M 
SCALE  -  SCALE  +  ABS(A(K,  I)) 

NEXT  K 

IF  SCALE  <>  0!  THEN 
FOR  K  =  I  TO  M 

A(K,  I)  =  A(K,  I) /SCALE 
S  -  S  +  A(K,  I)  *  A(K,  I) 

NEXT  K 
F  =  A(I,  I) 

G  =  -ABS(SQR(S))  *  SGN(F) 

H  =  F  *  G  -  S 
A(I,  I)  —  F  -  G 
IF  I  o  N  THEN 

FOR  J-LTON 
S  =  0! 

FOR  K  =  I  TO  M 
S  =  S  +  A(K,  I)  *  A(K,  J) 

NEXT  K 
F  =  S/H 
FOR  K  =  I  TO  M 
A(K,  J)  =  A(K,  J)  +  F  *  A(K,  I) 

NEXT  K 
NEXT  J 
END  IF 

FOR  K  =  I  TO  M 

A(K,  I)  =  SCALE  *  A(K,  I) 

NEXT  K 
END  IF 
END  IF 

W(I)  =  SCALE  *  G 
G  =  0! 

S  =  0! 

SCALE  =  0! 

IF  I  <=  M  AND  I  ON  THEN 
FOR  K  =  L  TO  N 
SCALE  =  SCALE  +  ABS(A(I,  K)) 

NEXT  K 
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IF  SCALE  <>  0!  THEN 
FOR  K  =  L  TO  N 

A(I,  K)  -  A(I,  K)  /  SCALE 
S  =  S  +  A(I,  K)  *  A(I,  K) 

NEXT  K 
F  =  A(I,  L) 

G  =  -ABS(SQR(S))  *  SGN(F) 

H  =  F  *  G  -  S 
A(I,  L)  =  F  -  G 
FOR  K  =  L  TO  N 

R V 1  (K)  =  A(I,  K)/H 
NEXT  K 
IF  I  o  M  THEN 

FOR  J  =  LTOM 
S  —  0! 

FOR  K  =  L  TO  N 
S  =  S  +  A(J,  K)  *  A(I,  K) 

NEXT  K 
FOR  K  =  L  TO  N 
A(J,  K)  =  A(J,  K)  +  S  *  RV1(K) 

NEXT  K 
NEXT  J 
END  IF 

FOR  K  =  L  TO  N 

A(I,  K)  =  SCALE  *  A(I,  K) 

NEXTK 
END  IF 
END  IF 

IF  ABS(W(I))  +  ABS(RV1(I))  >  ANORM  THEN  ANORM  =  ABS(W(I))  +  ABS(RV]  (I)) 
NEXT  I 

FOR  I  =  N  TO  1  STEP-1 
IF  I  <  N  THEN 
IF  G  <>  0!  THEN 
FOR  J  =  L  TO  N 

V(J,  I)  =  (A(I,  J)  /  A(I,  L))/G 
NEXT  J 

FOR  J  =  L  TO  N 
S  =  0! 

FOR  K  =  LTON 
S  -  S  +  A(I,  K)  *  V(K,  J) 

NEXTK 
FOR  K  =  L  TO  N 
V(K,  J)  =  V(K,  J)+S*V(K,  I) 

NEXTK 
NEXT  J 
END  IF 

FOR  J  ~  L  TO  N 
V(I,  J)  =  0! 

V(J,  I)  -  0! 

NEXT  J 
END  IF 
V(I,I)  =  1! 

G  =  RV1(I) 

L  =  I 
NEXT  I 

FOR  I  =  N  TO  I  STEP -I 
L  =  I  +  1 
G  =  W(I) 

IF  I  <  N  THEN 
FOR  J  =  L  TO  N 
A(I,  J)  =  0! 

NEXT  J 
END  IF 

IF  G  o  0!  THEN 
G  =  1 !  /  G 
IF  I  o  N  THEN 
FOR  J  =  L  TO  N 
S  =  0? 
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FOR  K  =  LTOM 
S  =  S  +  A(K,  I)  *  A(K,  J) 

NEXT  K 

F  =  (S  /  A(I,  I))  *  G 
FOR  K  =  I  TO  M 
A(K,  J)  =  A(K,  J)  +  F  *  A(K,  I) 

NEXT  K 
NEXT  J 
END  IF 

FOR  J  =  I  TO  M 
A(J,  I)  =  A(J,  I)  *  G 
NEXT  J 
ELSE 

FOR  J  =  I  TO  M 
A(J,  I)  =  0! 

NEXT  J 
END  IF 

A  (I,  I)  =  A(I?  I)+  1! 

NEXT! 

FOR  K  =  N  TO  I  STEP-1 
FOR  ITS  -  1  TO  30 
FOR  L  =  K  TO  1  STEP-1 
NM  =  L  -  1 

IF-ABS(RV1(L))  +  ANORM  -  ANORM  THEN  EXIT  FOR 
IF  ABS(W(NM))  +  ANORM  =  ANORM  THEN  EXIT  FOR 
NEXT  L 

IF  ABS(RV1(L))  +  ANORM  <>  ANORM  THEN 
C  =  0! 

S  =  1 ! 

FOR  I  =  L  TO  K 

F  =  S  *  RV1(I) 

RVI(I)  =  C*  RV1(I) 

IF  ABS(F)  +  ANORM  <>  ANORM  THEN 
G  =  W(I) 

H  =  SQR(F  *  F  +  G  *  G) 

W(I)  =  H 
H  =  1!  /  H 
C  =  G*H 
S  =  -F  *  H 
FOR  J  =  1  TO  M 
Y  =  A(J,  NM) 

Z  =  A(J,  I) 

A(J,  NM)  =  Y  *  C  +  Z  *  S 
A(J,  I)  =  -Y*S  +  Z*C 
NEXT  J 
END  IF 
NEXT  I 
END  IF 
Z  =  W(K) 

IF  L  =  K  THEN 
IF  Z  <  01  THEN 
W(K)  =  -Z 
FOR  J  =  I  TO  N 
V(J,  K)  *  -V(J,  K) 

NEXT  J 
END  IF 
EXIT  FOR 
END  IF 

IF  ITS  =  30  THEN  PRINT  "No  convergence  in  30  iterations":  ERASE  RVI :  END 
X  =  W(L) 

NM  -  K  -  1 
Y  =  W(NM) 

G  =  RV1(NM) 

H  ~  R  V 1  (K) 

F  =  ((Y  -  Z)  *  (Y  +  Z)  +  (G  -  H)  *  (G  +  H))  /  (2!  *  H  *  Y) 

G  =  SQR(F  *  F  +  1!) 

F  -  ((X  -  Z)  *  (X  +  Z)  +  H  *  ((Y  /  (F  +  ABS(G)  *  SGN(F)))  -  H»  /  X 
C  =  1! 
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S  —  1 ! 

FORJ-LTONM 
I  =  J  +  1 
G  =  RV1(I) 

Y  =  W(I) 

H  =  S  *  G 
G  =  C  *  G 

Z  -  SQR(F  *  F  +  H  *  H) 

RV1(J)  =  Z 
C  =  F/Z 
S  =  H/Z 
F=X*C+G*S 
G  =  -X*S  +  G*C 
H  -  Y  *  S 

Y  =  Y  *  C 

FOR  JJ  =  1  TON 
X  =  V(JJ,  J) 

Z  =  V(JJ,  I) 

V(JJ,  J)  =  X*C  +  Z*S 
V(JJ,  I)  =  -X  *  s  +  z  *  c 
NEXT  JJ 

Z  =  SQR(F  *  F  +  H  *  H) 

W(J)  =  Z 
IF  Z  <>  0!  THEN 
Z=l!/Z 
C  =  F  *  Z 
S  =  H*Z 
END  IF 

F  =  C  *  G  +  S  *  Y 
X  =  -S*G  +  C*Y 
FOR  JJ  =  1  TO  M 
Y  =  A(JJ,  J) 

Z  =  A(JJ,  I) 

A(JJ,  J)  =  Y*C  +  Z*S 
A(JJ,  I)  =  -Y  *  S  +  Z  *  C 
NEXT  JJ 
NEXT  J 
RV1(L)  =  0! 

RV I  (K)  =  F 
W(K)  =  X 
NEXT  ITS 
NEXT  K 
ERASE  RV1 
END  SUB 


SUB  test  (TEMP#) 

OPEN  filenames  +  ”.FILM  FOR  INPUT  AS  #1 


FOR  filter  =1  TO  4 
MLOLD1#  =  0 

FOR  lam#  =  lammin#(filtermax)  TO  LAMMAX#(filtermax)  STEP  DLAM# 
INPUT  #1,R# 

AOLD# -  0 

MLbbl#  =  cl#  /  (lam#  A  5  *  (EXP(c2#  /  (lam#  *  TEMP#))  - 1))  *  1 
ML1#  =  MLbbl#  *  R#  *  DLAM#  +  MLOLD1# 

MLOLD1#  =  ML1# 

NEXT  lam# 

emis  -  ML#(filter)  /  ML  I# 

PRINT  emis; 
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NEXT  filter 
PRINT  TEMP#  -  273 

CLOSE  #1 

END  SUB 


SUB  MRQMIN  (X(),  YO,  SIGO,  NDATA,  AO,  MA,  LISTA(),  MFIT,  COVAR(),  ALPHA(), 
NCA,  CHISQ,  DUM,  ALAMDA) 

STATIC 

DIM  ATRY(MA),  DA(MA,  1) 

IF  ALAMDA  <  0!  THEN 
KK  =  MFIT  +  1 
FOR  J  =  1  TOMA 
IHIT  -  0 

FORK=l  TO  MFIT 
IF  LISTA(K)  =  J  THEN  IHIT  -  IHIT  +  1 
NEXT  K 

IF  IHIT  =  0  THEN 
LISTA(KK)  =  J 
KK  =  KK  +  1 
ELSEIF  IHIT  >  1  THEN 
PRINT  "Improper  permutation  in  LISTA" 

EXIT  SUB 
END  IF 
NEXT  I 

IF  KK  <>  MA  +  I  THEN  PRINT  "Improper  permutation  in  LISTA":  EXIT  SUB 
ALAMDA  =  .001 

CALL  MRQCOF(X(),  Y(),  SIGO,  NDATA,  A(),  MA,  LIST  AO,  MFIT,  ALPHAO,  BETAO,  NCA,  CHISQ,  DUM) 

OCH1SQ  =  CHISQ 
FOR  J  = I  TO  MA 
ATRY(J)  =  A(J) 

NEXT  J 
END  IF 

FOR  J  =  1  TO  MFIT 
FOR  K  =  1  TO  MFIT 
COVAR(J,  K)  =  ALPHA(J,  K) 

NEXT  K 

COVARO,  J)  =  ALPHAO,  J)  *  (1!  +  ALAMDA) 

DAO,  1)  =  BETAO,  1) 

NEXT  J 

CALL  GAUSSJ(COVAR0,  MFIT,  NCA,  DAO,  1,  Q 
IF  ALAMDA  =  0!  THEN 

CALL  COVSRT(COVAR0,  NCA,  MA,  LISTA(),  MFIT) 

ERASE  DA,  ATRY 
EXIT  SUB 
END  IF 

FORJ=l  TO  MFIT 

ATRY(LISTAO))  =  A(LISTA(J))  +  DAO,  1) 

NEXT  J 

CALL  MRQCOF(X0,  YO,  SIGO,  NDATA,  ATRYO,  MA,  LISTAO,  MFIT,  COVARO,  DAO,  NCA,  CHISQ,  DUM) 

IF  CHISQ  <  OCHISQ  THEN 
ALAMDA  =  .  1  *  ALAMDA 
OCHISQ  =  CHISQ 
FOR  J  =  1  TO  MFIT 
FOR  K  =  1  TO  MFIT 
ALPHAO,  K)  =  COVARO,  K) 

NEXTK 

BETAO,  1)  =  DAO,  1) 

A(LISTAQ))  =  ATRY (LISTAO)) 

NEXT  J 
ELSE 

ALAMDA  =  10!  *  ALAMDA 
CHISQ  =  OCHISQ 
END  IF 

ERASE  DA,  ATRY 
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END  SUB 


A3.  Reconstruct  Program: 


DECLARE  SUB  Normalize  (a#0) 

DECLARE  SUB  PLOTDATA  (ATARGET#()) 

DIM  SHARED  o  AS  STRING  *  1 
DIM  SHARED  ML#(10) 

COMMON  SHARED  DLAM#,  FILTERMAX,  LAMMAX#(),  LAMMIN#() 

SCREEN  9 

CLS 

DIM  a# (5,  500),  lo(4),  med(4),  hi(4),  MLhigh#(4),  MLIow#(4),  MLmid#(4) 
DIM  LAMMIN#(5),  LAMM  AX#  (5) 

EMISS#  =  I 
FILTERMAX  =  4 
DLAM#  =  .005 
LAMMIN#(1)  =  2.7 
LAMMAX#(1)  =  5.2  -  DLAM# 

LAMMIN#(2)  =  2.7 
LAMM  AX#  (2)  =  5.2  -  DLAM# 

LAMMIN#(3)  =  2.7 
LAMM  AX#  (3)  =  5.2  -  DLAM# 

LAMMIN#(4)  =  2.7 
LAMM  AX#  (4)  =  5.2  -  DLAM# 

LAMMIN#(5)  =  2.7 
LAMM  AX#  (5)  =  5.2  -  DLAM# 

x$  =  STR$(1) 

x$  =  MID$(x$,  2) 

f$  =  "a:\monocr"  +  x$  +  ".arr" 

PRINT  f$ 

OPEN  f$  FOR  BINARY  AS  #1 
GET #1,  l,o 
a!  =  ASC(o) 

GET  #1,  a!,  o 
GET#1,,  OF 
FOR  I  =  1  TO  499 
GET  #  1 , ,  m 
a#(l,  I)  =  m 
NEXT  I 
CLOSE  #1 

LINE  (5,  l)-(5,  100),  2 
LINE  (1,  100M400,  100),  2 

FOR  1  =  10  TO  100  STEP  10:  LINE  (1, 1)-(10, 1),  2:  NEXT  1 


1  =  0 

PSET(5,  100) 

FOR  lam#  =  LAMMIN#(1)  TO  LAMMAX#(1)  STEP  DLAM#  ’Calculate  and  store  emissivities  Plot  filter  values 
1  =  1+  1 

IF  I  <  55  THEN  a#(l,  I)  =  a#(l,  55)  *  I  /  55 

IF  I  >  294  AND  I  <  336  THEN  a#(l,  I)  =  a#(l,  294)  -  (a#(l,  294)  -  a#(l,  336))  *  (I  -  294)  /  40 
IF  I  >  435  THEN  a#(l,  I)  =  a#(l,  435)  -  a#(l,  435)  *  (I  -  435)  /  64 
x  =  (lam#  -2.7)  *200 
Y  =  (1  -  .1  *a#(I,I)>*  100 
LINE  -(x,  Y),  5 

’  LOCATE  20,  1 :  PRINT  lam#;  i; :  INPUT  a$ 

NEXT  lam# 

LOCATE  20,  1:  INPUT  a$ 

OPEN  f$  FOR  BINARY  AS  #  1 
OPEN  "a:\temp. arr"  FOR  BINARY  AS  #2 
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GET  #1,  l,o 
a!  =  ASC(o) 

GET  #1,  a!,  o 
GET #1,,  OF 
FOR  1=1  TO  499 
PUT  #2, ,  a#(l,  1) 

NEXT  I 
CLOSE  #1 
CLOSE  #2 

SUB  calcML  (a#0,  MLhigh#(),  MLlow#()) 

Cl#  =  3.7415  *  10  A  4 
C2#  =  1.43879*  10  A  4 
C3#  =  1.88365  *  10  A  23 

INPUT  "input  low  temp'1,  temp# 

FOR  FILTER  =  1  TO  F1LTERMAX 
MLOLD#  =  0 
1  =  0 

FOR  lam#  =  LAMMIN#(FILTER)  TO  LAMMAX#(FILTER)  STEP  DLAM# 
1  =  1  +  1 

R#  =  a#  (FILTER,  I) 

MLBB#  =  Cl#  /  (lam#  A  5  *  (EXP(C2#  /  (lam#  *  temp#))  - 1))  *  EMISS# 
MLlow#(FILTER)  =  MLBB#  *  R#  *  DLAM#  +  MLOLD# 

MLOLD#  =  MLIow#(FILTER) 

NEXT  lam# 

NEXT  FILTER 

INPUT  "input  high  temp”,  temp# 

FOR  FILTER  =  1  TO  FILTERMAX 
MLOLD#  =  0 
1  =  0 

FOR  lam#  =  LAMMIN#(F1LTER)  TO  LAMMAX#(FILTER)  STEP  DLAM# 
1  =  1+1 

R#  =  a#  (FILTER,  I) 

MLBB#  =  Cl#  /  (lam#  A  5  *  (EXP(C2#  /  (lam#  *  temp#))  - 1))  *  EMISS# 
MLhigh#(FILTER)  =  MLBB#  *  R#  *  DLAM#  +  MLOLD# 

MLOLD#  =  MLhigh#(FILTER) 

NEXT  lam# 

NEXT  FILTER 

END  SUB 

SUB  gainoff(a#0,  Io(),  mid(),  hit),  MLIow#0,  MLhigh#0,  MLmid#()) 


'INPUT  "input  data  file  name  (samplel)  and  no.  in  file  (0-2)'',  filenames,  no! 
filenames  =  "samplel.arr":  no!  =  2 
DIM  dfilter(4) 

OPEN  filenames  FOR  BINARY  AS  #3 

GET  #1,  l,o 

a!  =  ASC(o)  +  4!  *  no! 

GET  #  I ,  a!,  o 


OPEN  "a:\DATA.FIL"  FOR  OUTPUT  AS  #1 

FOR  FILTER  -  1  TO  FILTERMAX  'Gain  and  offset  the  data  using  low  and  high  blackbody  values 
Gall  -  (hi (FILTER)  -  lo(FILTER))  /  (MLhigh#  (FILTER)  -  MLlow# (FILTER)) 
ofcet  =  lo(FILTER)  -  MLlow#  (FILTER)  *  Gall 
GET  #3, ,  ML 

ML#(FILTER)  =  (ML  -  ofcet)  /  Gall 
1  =  0 

FOR  lam#  -  LAMMIN# (FILTER)  TO  LAMMAX# (FILTER)  STEP  DLAM# 

1  =  1+  1 

PRINT  #1,  a#  (FILTER,  I) 

NEXT  lam# 

NEXT  FILTER 
CLOSE  #1 
CLOSE  #3 
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INPUT  "input  mid  temp",  t  'Calculate  errors  using  3rd  blackbody  value 
FOR  temp#  =  t  -  1 0  TO  t  +  1 0 
FOR  FILTER  =  1  TO  FILTERMAX 
MLOLD# = 0 
1  =  0 

FOR  lam#  =  LAMMIN#  (FILTER)  TO  LAMM  AX#  (FILTER)  STEP  DLAM# 
1  =  1+1 

R#  =  a#(FILTER,  I) 

MLBB#  =  Cl#  /  (lam#  A  5  *  (EXP(C2#  /  (lam#  *  temp#))  -  1))  *  EMISS# 
MLmid# (FILTER)  =  MLBB#  *  R#  *  DLAM#  +  MLOLD# 

MLOLD#  =  (MLmid# (FILTER)  -  ofcet)  /  Gall 
NEXT  lam# 

NEXT  FILTER 
FOR  I  =  1  TO  FILTERMAX 
dfilter(I)  =  MLmid#(I)  -  mid(I) 

NEXT  I 

PRINT  temp#,  dfilter(l),  dfilter(2),  dftlter(3),  dfilter(4) 

NEXT  temp# 

OPEN  "DATA.SYS"  FOR  OUTPUT  AS  #2 
PRINT  #2,  FILTERMAX 
PRINT  #2,  DLAM# 

FOR  I  =  1  TO  FILTERMAX:  PRINT  #2,  LAMMIN#(I):  NEXT  I 
FOR  I  =  1  TO  FILTERMAX:  PRINT  #2,  LAMM  AX#  (I):  NEXT  I 
FOR  I  =  1  TO  FILTERMAX:  PRINT  #2,  ML#(I):  NEXT  I 
CLOSE  #2 
END  SUB 

SUB  getbbdata  (FILTERMAX,  lo(),  med(),  hi()) 

OPEN  "a:\blackbol. arr"  FOR  BINARY  AS  #1 
GET  #  1 ,  l,o 
a!  =  ASC(o) 

GET  #1,  a!,  o 

FOR  I  =  1  TO  FILTERMAX 
GET  Io(I),  med(I),  hi(I) 

NEXT  I 
CLOSE #1 

END  SUB 

SUB  Normalize  (a#()) 


END  SUB 

SUB  PLOTDATA  (a#Q) 


END  SUB 


Appendix  B- Amber  Proview  Software 


The  programs  contained  in  this  appendix  have  been  developed  in  the  Amber  ProView 
TDF  language.  There  were  used  to  interface  the  Amber  256  x  256  snapshot  camera  and  there 
purpose  is  basic  data  collection. 


Bl.  Calibrate  Monochromator  Light  Source: 


settimeout(120) 

arr  discard(alO) 

arr  discard(al  1) 

arr  discard(al2) 

arr  discard(al3) 

arr  create(al3,2,  500) 

arr  create(alO,  256,  256) 

arr  create(all,  5,  500) 

arrcreate(al2,  5, 3) 

all  =  0.0 

al2  =  0.0 

al3  =  0.0 

i5=0 


loop  (500) 


i5“  i5  +  1 
s5  =  chr(13) 

sO  =  7700.0  GOTO'  +  s5 
gpib  write(l  l,s0) 
delay(  15000) 


r5=1.0 
11=0 
loop  (500) 
acquire(l,l,N) 
alO  =  up!oad(cl,  0,  65536) 
al3[l,Il]  =  al0[l,  127] 
al3[l,Il]  =  al3[l,Il]  +  al0[l,  128] 
al3[l,Il]  =  a!3[l,Il]  +  al0[l,  129] 
al3[l,Il]  =  al3[l,Il]  +  a!0[2,  127] 
al3[l,Il]  =  al3[l,Il]  +  al0[2,  128] 
al3[l,Il]  =  a!3[l,Il]  +  al0[2,  129] 
al3[l,Il]  =  al3[l,Il]  +  al0[3,  127] 
al3[l,Il]  =  al3[l,Il]  +  al0[3,  128] 
al3[l,Il]  =  al3[l,Il]  +  al0[3,  129] 
if  (il>  0  ) 

r5=al3[l,  il]/al3[l,0] 
endif 

message(l,30,r5) 
al3[0,ll]  =  al0[127,  222] 
al3[0,Il]  =  al3[0,Il]  +  al0[127,  223) 
al3[0,Il]  =  al3[0,Il]  +  a!0[127, 224] 
al3[0,Il]  =  al3[0,Il]  +  a!0[128,  222] 
al3[0,Il]  =  a!3[0,Il]  +  al0[128,  223] 
al3[0,Il]  =  al3[0,Il]  +  al0[128,  224] 
al3[0,Il]  =  al3[0,Il]  +  al0[126,  222] 
al3[0,Il]  =  al3[0,Il]  +  al0[126,  223] 
al3[0,Il]  =  al3[0,Il]  +  a!0[126,  224] 
if  (il>0) 
r6=al3[0,0]  *  r5 
al3[0,Il]  =  al3[0,il]  -  r6 
endif 
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11=11  +  1 

12  =  II  *  5 

13  =  i2  +  2700 
s2  =  str(i3) 

sO  =  s2  +  '.0  GOTO' 
sO  =  sO  +  s5 
gpib  write(l  l,sO) 

14  =  II  - 1 
message  (l,l,i4) 
message  (2,l,i3) 
message  (4,l,al3[0,i4]) 
message  (1, 20,  i5) 

endloop 

si  =  str(i5) 

s9  =  'b:\monocr'  +  si 

arr  save  file(al3,s9) 

endloop 


B2.  Get  Data  For  Filter  Bandpasses: 


set  timeout(120) 

arr  discard(alO) 

arr  discard(all) 

arr  discard(al2) 

arr  discard(al3) 

arr  create(a!3,2,  500) 

arr  create(al0,  256,  256) 

an-create(all,  5,  500) 

arr  create(al2,  5 ,  3) 

all  =0.0 

al2  =  0.0 

a!3  =  0.0 


s5  =  chr(13) 

sO  =  '2700.0  GOTO'  +  s5 
gpib  write(ll,sO) 
si -Wait  for  Monocrometer' 
message  (l,l,sl) 
read(sl„) 

;echo  (off) 


II  =0 
r5  =  1.0 
loop  (500) 

acquire(l,l,N) 
alO  =  upload(cl,  0,  65536) 
all[4,Il]  =  al0[l,  127] 
all[4,Il]  =  all[4,Il]  +  al0[l,  128] 
all [4,  II]  =  all[4,  II]  +  al0[l,  129] 
al 1[4,  II]  =  all [4,  II]  +  al0[2,  127] 
al  1[4,  II]  =  al  1[4,  II]  +  al0[2,  128] 
al  1  [4,  II]  =  al  1[4,  II]  +  al0[2,  129] 
all  [4,  II]  =  all  [4,  II]  +  al0[3, 127] 
al  1  [4,  II]  =  al  1[4,  II]  +  al0[3,  128] 
al  1  [4,  II]  =  al  1  [4,  II]  +  a!0[3, 129] 
if  (il>  0  ) 

r5=al  1  [4,  i  1  ]/a  1 1  [4,0] 
endif 
rl0=l-r5 
message(l,30,r5) 

all[0,  II]  =  al0[127,  30] 

al  1[0,  II]  =  all[0,  II]  +  al0[127,  31] 

al  1[0,  II]  =  al  1[0,  II]  +  al0[127,  32] 
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all[0,  II]  -  al  I [0,  II]  +  al0[128,  30] 
all[0,  II]  -  all [0,  II]  +  al0[128,  31] 
all[0,  II]  =  all [0,  II]  +a!0[128,  32] 
all[0,  II]  =  al  1[0,  II]  +  al0[126>  30] 
al  1[0,  II]  =  all[0,  II]  +  al0[I26,  31] 
al  1  [0,  II]  =  al  1  [0,  II]  +  al0[126,  32] 
if  (il>0) 

r6 -  all [0,0]  *  r5 
all[0,Il]  =  all[0,Il]-r6 
endif 

al  1  [  1 ,  II]  =  al0[127,  94] 

all[l,  II]  =  all[l,  II]  +  al0[127,  95] 

all[l,Il]  =  all[l,Il]+alO[127,  96] 

all[l,Il]  =  all[l,II]  +  alO[128,  94] 

all[l,  Il]  =  all[l,  II]  +  al0[128,  95] 

al 1[1,  II]  =  all[l,  II]  +  al0[128,  96] 

all[l,  II]  =  all[l,  II]  +  al0[126,  94] 

all[l,  II]  —  all [1,  II]  +  al0[126,  95] 

all[l,Il]  =  all[I,Il]  +  al0[126,  96] 

if  (i  1  >0  ) 

r7=al  1  [1,0]  *  r5 

al  1[1,  II]  -  al  1  [1 ,  II]  -  r7 

endif 

al  1  [2,  II]  =  al 0[1 27,  158] 

al  I  [2,  II]  =  al  1  [2,  II]  +  al0[127,  159] 

all [2,  II]  -  al  1[2,  II]  +  al0[127,  160] 

all [2,  II]  =  all [2,  II]  +  al0[128,  158] 

al  1  [2,  II]  =  al 1[2,  II]  +  al0[128,  159] 

al  1  [2,  II]  -  al  1  [2,  II]  +  al0[128,  160] 

all [2,  II]  =  all [2,  II]  +  a!0[126,  158] 

al  1[2,  II]  =  all[2,  II]  +  al0[126,  159] 

al  1[2,  II]  =  al  1[2,  II]  +  a!0[126,  160] 

if  (il>0) 

r8=al  1  [2}0]  *  r5 

al  1  [2, 11]  -  al  1[2,  II]  -  r8 

endif 

all [3,  II]  -  al 0[1 27,  222] 

al  1[3,  II]  =  al  1[3,  II]  +  al0[127,  223] 

al  1[3,  II]  =  al  1[3, 11]  +  al0[127,  224] 

all [3,  II]  =  all [3,  II]  +  al0[128,  222] 

al 1[3,  II]  =  al  1[3,  II]  +  al0[128,  223] 

al  1  [3,  II]  =  al  1  [3,  II]  +  al0[128,  224] 

all [3?  II]  =  al 1[3,  II]  +  a!0[126,  222] 

al  1[3,  II]  =  al  1[3,  II]  +  a!0[126,  223] 

al  1  [3,  II]  ==  all [3,  II]  +  al0[126,  224] 

if  (i  1  >0) 

r9=al  1[3,0]  *  r5 

al 1[3, II]  =  all[3, II]  -  r9 

endif 

11  =11  +  1 

12  =  II  *  5 

13  =  12  +  2700 
s2  =  str(i3) 

s0  =  s2 +  ’.0  0010’ 
sO  =  sO  +  s5 
gpib  write(ll,sO) 

14  =  II  - 1 
message  (l,l,i4) 
message  (2,l,i3) 
message  (3,l,all[0,  i4]) 
message  (4,l,all[l,  i4]) 
message  (5,l,all[2,  i4]) 
message  (6,l,all[3,  i4]) 
message  (7,l,all[4,  i4]) 
rl=al  1[0,  i4]-al  1[1,  i4] 
r2=all[l,  i4]-all[2,  i4] 
r3=all[2,  i4]-all[3,  i4] 
message  (9,l3rl) 


message  (10,1,  r2) 
message  (ll,l,r3) 
endloop 
s9  -  'b:\filterr 
arr  save  file(al  l,s9) 

B3.  Calibrate  Monochromator 

set  timeout(120) 

arr  discard(alO) 

arr  discard(all) 

arr  discard(al2) 

arr  discard(a!3) 

arr  create(al3,2,  500) 

arr  create(alO,  256,  256) 

arr  create(all,  5,  500) 

arr  create(al2,  5,  3) 

all  =0.0 

al2  =  0.0 

al3  —  0.0 


s5  =  chr(13) 

sO  =  '2700.0  GOTO'  +  s5 
gpib  write(ll,s0) 
s  1  -  Wait  for  Monocrometer' 
message  (l,l,sl) 
read(sl„) 

;echo  (off) 


r5=1.0 
II  =0 
loop  (500) 
acquire(l,l,N) 
a!0  =  upload(cl,  0,65536) 
al3[l,Il]  =  a!0[l,  127] 
al3[l,Il]=al3[l,Il]+alO[l,  128] 
al3[l,Il]  =  al3[l,Il]  +  al0[l,  129] 
al3[l,Il]  =  al3[l,Il]  +  al0[2,  127] 
al3[l,Il]  =  al3[l,Il]  +  al0[2,  128] 
al3[l,Il]  =  al3[l,Il]  +  al0[2,  129] 
al3[l,Il]  =  al3[l,Il]  +  al0[3,  127] 
al3[l,Il]  =  al3[l,Il]  +  al0[3,  128] 
a!3[l,Il]  =  al3[l,Il]  +  a!0[3,  129] 
if  (il>  0  ) 

r5— a  1 3 [  1 ,  il]/al3[l,0] 
endif 

message(l,30,r5) 
al3[0,Il]  “  al0[127,  222] 
al3[0,Il]  =  al3[0,Il]  +  al0[127,  223] 
a!3[0,Il]  =  al3[0,Il]  +  al0[127, 224] 
al3[0,Il]  =  a  1 3  [0,1 1  ]  +  al0[128,  222] 
al3[0,Il]  -  al3[0,Il]  +  al0[128,  223] 
al3[0,Il]  =  al3[0,Il]  +  al0[128,  224] 
al3[0,Il]  =  al3[0,Il]  +  al0[126, 222] 
al3[0,Il]  -  al3[0,Il]  +  al0[126,  223] 
al3[0,Il]  =  al3[0,Il]  +  al0[126,  224] 
if  (il>0) 
r6=al3[0,0]  *  r5 
al3[0,Il]=al3[0,il]-r6 
endif 

11  =11  +  1 

12  =  II  *  5 

13  =  i2  +  2700 
s2  =  str(i3) 

s0  =  s2 +  ’.0  00X0' 
sO  =  sO  +  s5 
gpib  write(l  l,sO) 
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i4  =  Il  -1 
message  (l,l,il) 
message  (2,l,i3) 
message  (4,l,al3[0,i4]) 
endloop 

s9  =  'b:\monociT 
arr  save  file(al3,s9) 

□ 

i  1=0 

set  timeout  (60) 
arr  discard  (alO) 
arr  create  (al  0,256,256) 

loop  (100) 


il=il+l 
acquire(l,l,N) 
al0=upload(cl,0, 65536) 
s8=str(il) 
s9='b:\test'  +  s8 
;  arr  save  fi!e(a!0,s9) 
greyscale(al0,6,nowait,noclear) 
read  (sl„) 
if  (si  in  ('q')) 
exitloop 
endif 

message  (l,l,s9) 
endloop  □ 


B4.  Get  Blackbody  And  Sample  Data: 


set  timeout(120) 
arr  discard(alO) 
arr  discard(a!2) 
arr  discard(al3) 
arr  discard(al4) 

arr  create(al0, 256,  256) 
arrcreate(al2,  5,  10) 
arr  create(al3,  5,  10) 
arr  create(al4,  5,  10) 
al2  =  0.0 
al3  =  0.0 
al4  =  0.0 


i0=0 

loop(10) 

si- Set  Blackbody  at  temp' 
message  (l,l,sl) 
read(sl„) 
acquire(l,l,N) 
alO  =  upload(cl,  0,  65536) 
al2[0,  iO]  -  al0[127,  31] 
al2[l,  iO]  =  al0[127,  95] 
al2[2,  iO]  ~  a!0[127,  159] 
al2[3,  iOj  =  al0[127,  223] 
al2[4,  iO]  =  alO[l,  1] 

si  - Set  Sample  at  temp' 
message  (I,l,sl) 
read(sl„) 
acquire(l,l,N) 
alO  =  upload(cl,  0,  65536) 
al3[0,  iO]  =  al0[127,  31] 
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al3[l,i0]  =  a!0[127,  95] 
al3[2,  iO]  =  a!0[127,  159] 
al3[3,  i0]  =  al0[127,  223] 
a!3[4,  iO]  =  alO[l,  1] 


si  - Set  Sample  at  temp' 
message  (1,1, si) 
read(sl„) 
acquire(l,l,N) 
alO  =  upload(cl,  0,  65536) 
al4[0,  iO]  =  alO[I27,  31] 
al4[l,i0]  =  al0[127,  95] 
al4[2,  iO]  =  al0[127,  159] 
al4[3,  iO]  ~  al0[127, 223] 
al4[4,  iO]  =  alO[l,  1] 
iO=iO+l 


endloop 

s9  =  'b:\bIackbor 
arr  save  file(al2,s9) 

s9  =  'b:\samplel' 
arr  save  file(al3,s9) 

s9  =  'b:\sample2' 
arr  save  file(a!4,s9) 
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ABSTRACT 

Real  time  analysis  of  the  IR  spectrum  is  becoming  increasingly  important  for  the  accurate  measurement  of  the 
emissive  characteristics  of  materials.  The  objective  is  once  the  emissivity  is  defined  the  temperature  across  a  surface  can  be 
observed  using  passive  IR  techniques.  It  is  both  important  to  create  accurate  emissivity  models  as  well  as  develop 
instrumentation  that  can  acquire  data  accurately.  With  the  advent  of  large  high  speed  infrared  detector  arrays  it  has  become 
possible  to  realize  high  resolution  imaging  spectrometers.  To  develop  a  model  for  materials  such  as  metals  we  observe  that 
the  primary  contribution  to  absorption  characteristics  can  be  attributed  to  the  free  electron  plasma.  It  can  be  shown  that  for 
heated  metals  in  thermal  equilibrium  the  emissivity  can  be  modeled  as  a  second  order  polynomial.  Consequently,  a  four  color 
pyrometer  can  be  used  to  accurately  determine  emissivity  profile  so  that  the  true  temperature  can  be  determined. 

1.  INTRODUCTION 


Over  the  past  few  years  I  have  been  working  in  the  area  of  energetic  materials  pyrometry.  I  have  found  a  great  need 
for  accurate  temperature  measurement  for  events  such  as  the  internal  combustion  in  engines;  materials  deformation  in  metal 
processing;  and  detonating  or  defligrating  explosives.  For  these  environments  the  objects  that  we  are  interest  in  are  moving  at 
high  rates.  In  addition  the  temperature  may  also  be  moving  at  high  rates;  on  the  order  of  hundreds  of  degrees  per  millisecond. 
After  an  extensive  survey  of  the  technologies  that  could  possibly  solve  some  of  this  need  I  found  little  that  is  applicable  to  the 
energetic  environment. 

^  In  recent  years  the  inclusion  of  temperature  effects  in  high  speed  deformation  models  has  become  more  critical.  The 
Wright  laboratory  armament  directorate  has  been  working  to  develop  accurate  models  that  predict  penetrator  type  weapons 
deformation  as  the  projectile  encounters  harden  structures.  The  work  presented  in  this  paper  is  in  the  support  of  temperature 
profilometry  of  the  metal  casing  material  used  in  penetrator  weapons.  It  is  the  intent  of  this  paper  to  develop  the  theory 
leading  to  an  accurate  model  of  emissivity  and  to  substantiate  this  model  with  experimental  data. 

2.  BACKGROUND 

As  energy  is  delivered  to  a  stricture  at  first  it  undergoes  elastic  deformation  and  will  return  to  its  original  condition 
after  the  force  is  released.  When  enough  force  has  been  applied  the  material  it  reaches  its  “yield  point”  and  the  structure  will 
deform  plastically  until  it  breaks.  After  the  yield  point  is  reached  the  stricture  will  not  return  to  its  original  condition  after  the 
force  is  released.  Consequently,  work  has  been  done  to  the  structure  in  the  form  of  material  realignment  and  the  generation  of 
heat  energy.  In  fact  most  of  the  deformation  energy  gets  turned  into  heat  and  this  heat  energy  modifies  how  the  material  will 
plastically  deform.  If  the  emissivity  of  the  material  is  known  radiometric  techniques  can  be  used  to  determine  the  surface 
temperature  passively. 

Methods  for  determining  the  surface  temperature  of  heated  material  using  optical  pyrometers  has  often  been  based  on 
assumptions  that  lead  to  erroneous  results.  The  primary  problems  here  is  the  uncertainties  in  the  collected  optical  energy. 
Emissivity,  atmospheric  transmission,  detector  non-linearities  etc.  can  all  affect  the  amount  of  energy  collected  depending  on 
the  wavelength  and  bandwidth  the  measurement  was  made.  Often  the  problem  can  be  reduced  to  an  uncertainty  in  Emissivity1. 
Emissivity  is  often  characterized  with  a  M---order  polynomial-with - - - - 


eefffy)~aM+i +  a m\m  1  +aM_x\M  2+ +av 
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where  the  a  s  are  unknown  constants  that  describe  emissivity  behavior.  If  for  instance  a  given  material  had  a  constant 
emissivity  we  would  need  to  determine  ax  where  all  the  other  a’s  would  be  identically  zero.  Consequently,  we  would  require 
two  spectral  measurements  to  determine  temperature.  Accordingly,  if  an  Mth  order  polynomial  provides  a  given  materials 
emissivity  profile  we  will  require  a  minimum  of  M+2  spectral  measurements.  However,  the  instrumentation  errors  increase 
strongly  as  M  increases  ,  consequently,  it  is  very  important  to  limit  the  number  of  spectral  measurements  by  use  of  simple 
models. 


3J3LACKBODY  RADIATION 

Light  emitted  from  solids  under  conditions  of  thermal  equilibrium  and  in  the  absence  of  any  other  external  sources,  is 
known  as  thermal  light.  In  this  section  I  will  attempt  to  validate  the  thermal  equilibrium  assumption  for  the  high  speed 
application.  The  thermal  light  derivation  done  here  will  closely  follow  that  which  is  done  in  Fundamentals  of  Photonics3. 
Consider  a  cavity  of  unit  volume  whose  walls  have  many  atoms  that  have  two  energy  states  that  are  separated  by  energy 
difference  hv.  The  cavity  will  support  braudband  radiation.  Now  let  Nx  and  N2  representing  the  number  of  atoms  per  unit 
volume  occupying  energy  levels  one  and  two.  Spontaneous  emission  creates  radiation  in  the  cavity,  assuming  one  atom  is  at 
energy  level  two.  Subsequently  absorption  and  stimulated  emission  can  occur.  The  three  states  exist  in  equilibrium  when 
steady  state  is  reached.  We  assume  that  there  are  an  average  number  of  photons  per  radiation  mode  ft  . 

We  first  consider  the  spontaneous  emission  rate  and  then  we  will  move  on  to  absorption  and  stimulated  emission. 
The  probability  that  an  atom  in  the  upper  energy  state  undergoes  spontaneous  emission  into  any  of  the  modes  within  the  time 

increment  from  t  to  t  +  At  is  P^At  =  At/ .  There  are  N2  such  atoms  and  the  average  number  of  photons  emitted  in  At 

seconds  is  therefor  N^t/t^  .  This  is  also  the  number  of  atoms  that  leave  level  2  during  the  time  interval.  Consequently  the 
rate  that  atoms  increase  in  level  2  due  to  spontaneous  emission  is  given  by, 


dN^ 

dt 


The  solution  is  a  exponential  decaying  function  as  a  function  of  initial  state  at  time  equals  zero. 

f~] 

N1=N0/'*) 


In  the  presence  of  radiation  absorption  and  stimulated  emission  will  occur.  This  changes  the  populations  of  N}  and 
N2  ■  Lets  consider  the  absorption  first.  Since  there  are  jVj  in  level  one  capable  of  absorbing,  the  rate  of  increase  of  the 
population  of  atoms  in  the  upper  level  due  to  absorption  is, 


dN2  Nxn 


dt 


In  a  similar  way  the  stimulated  emission  can  be  written  as, 


sp 


dN 2 
dt 


N2n 


sp 


As  you  can  see  the  absorption  and  spontaneous  rates  are  proportional  to  the  average  number  of  photons  in  each  mode. 


We  can  now  write  the  rate  equation  for  the  overall  process  as. 
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dN2 

~dt 


N2  N^n*  N2n 
t  t  t 

sp  sp  1  sp 

Notice  that  this  equation  does  not  include  transitions  to  other  than  the  two  levels,  however,  adding  more  levels  is  a  trivial 
extension  when  we  are  in  thermal  equilibrium.  In  steady  state  the  rate  of  change  the  levels  will  be  zero.  Consequently,  we 
have 


Nl  l  +  o 

We  know  that  we  are  in  thermal  equilibrium  and  therefore  we  know  that  the  populations  obey  a  Boltzman 
distribution,  i.e.,  1 


*1 


and  we  can  now  write, 


1 


The  average  energy  E  of  each  mode  in  this  situation  will  be  simply  nhv  ,  so  that 


E  = 


hv 


Multiplying  this  expression  by  the  modal  density  87tV  2/c3  we  have  the  blackbody  radiation  law 


Mv  = 


Snhv 2 


r  hv  ) 

Sw) 


-1 


We  have  derived  the  relation  above  under  the  assumption  that  we  could  assume  thermal  equilibrium.  Let  me  now  see 
is  this  assumption  is  a  good  one.  If  the  temperature  changes  from  T}  to  T2  During  the  time  interval  At  the  change  in  the 
populations  will  be  related  as, 


AN2  =  -ANl 

.  If  ten_atom  move  to  level .  two L_then  level  one  looses  ten  atoms. _We  can  also  write 


AN, 


n 


1  +  n 


— — — 
l  +  «  1 
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For  our  purposes  we  will  be  looking  at  objects  at  about  500  K.  4f  the  temperature  raises  by  100°  and  our  operational 
wavelength  is  at  three  microns  we  can  write 


^“KL-foo-Kikio. 

We  can’t  solve  this  equation  without  knowledge  of  the  population  of  level  one,  however,  we  can  ballpark  the  range  of  A N2 . 

In  the  equation  above  we  know  that  the  first  term  is  larger  than  the  second  term  and  that  they  are  both  positive.  Dropping  the 
second  term  we  can  write 


We  can  also  write  the  population  of  level  two  in  terms  of  level  one  as 


and 


an2  <n2 


Consequently,  if  the  time  of  our  observation  is  one  microsecond  and  the  relaxation  time  is  10  14  seconds  then  the  changing 
term  is  negligible,  i.e.. 


and 


dN 2  _  N2  N{n  N2n  N2 

dt  tsp  tsp  tsp  ~  tsp 


dN ,  =  10~9  Nn 


~sp 


Therefore, 


dN2 

dt 


=  0 


and  the  steady  state  assumption  is  validated.  It  appears  as  though  the  blackbody  radiation  law  will  hold  true  until  the 
relaxation  time  approaches  the  camera  integration  time. 

4.  EMISSIVITY  USING  A  FREE  F.T.F.CTRON  MODEL 

We  have  looked  the  characteristic  manner  that  a  ideal  solid  generates  its  IR  spectrum  and  now  I  divert  our  attention  to 
the  characterization  of  a  metals  absorption  spectrum  and  consequently  its  emissivity.  Here  I  will  model  the  metal  as  an 
electron  plasma  with  characteristic  relaxation  time  and  effective  mass.  This  model  does  show  promise  to  accurately  model 
emissivity.  In  equilibrium  according  to  Kirkoffs  law  the  absorption  of  a  material  is  equivalent  to  the  emissivity. 
Consequently  if  we  can  determine  the  absorption  characteristics  of  a  materiel  we  know- the  emissivity  characteristics,  atleast  if 
the  temperature  is  not  changing  quickly.  Since  much  of  our  application  is  considering  the  properties  of  heated  metals  I  will 
obtain  the  absorption  relationship  for  a  simple  metal.  Test  results  will  support  the  accuracy  of  this  technique. 

The  earliest  attempts  to  predict  the  optical  properties  of  metals  were  made  by  Lorentz  and  Durde4,  who  assumed  the 
metal  to  contain  electrons  which  were  essentially  free  to  move  under  the  influence  of  the  electric  field  induced  by  the  incident 
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electromagnetic  wave.  These  free  electrons  are  the  valence  electrons  in  the  outer  shell  of  the  atoms  constituting  the  metal. 
When  the  wave  is  incident  upon  its  surface,  an  oscillating  electric  field  parallel  to  the  surface  is  induced  in  the  metal  and  the 
free  electrons  will  oscillate  under  the  influence  of  this  field  at  the  frequency  of  the  incident  wave.  To  describe  the  optical 
behavior  requires  two  parameters:  the  number  density  of  free  electrons,  n,  being  excited  by  the  induced  field,  and  the  average 
time  (relaxation  time  t)  between  collisions  of  the  electron  with  the  atomic  lattice.  This  is  called  the  Drude  Free  Electron  model. 
If  the  assumption  can  be  made  that  cot«1  then  a  simplified  model  can  be  generated. 

For  many  metals,  particularly  at  higher  temperatures  and  for  relatively  high  purity,  the  electrical  receptivity  is  proportional 
to  temperature;  through  integration  of  the  emissivity  equation  it  can  be  shown  that  the  normal  total  emissivity  is  proportional  to 
the  first  power  of  the  temperature5.  To  obtain  the  absorption  characteristic  we  need  to  determine  the  polarization  of  the 
material  from  which  we  can  derive  the  dielectric  function.  Using  some  well  known  relations  for  coefficient  of  absorption  in 
terms  of  dielectric  function  we  can  derive  the  dependence  with  wavelength.  Let  us  first  assume  that  the  metal  crystal  structure 
consists  of  charges  (electrons)  that  respond  to  an  electric  field,  E(t),  polarized  in  the  x  direction.  We  can  write  the  equations  of 
motion  for  each  electron  in  the  field  using  Newton’s  second  law 

dt  x  dt 

where  x(t)  is  the  electron  displacement  from  equilibrium,  X  is  the  relaxation  time  and  me  is  the  effective  mass  of  an  electron. 
Notice  we  have  neglected  the  nucleus  motion  effect  known  as  the  dipole  approximation.  The  second  term  of  the  above 
equation  represents  the  damping  due  to  collations  with  ions.  Notice  that  since  the  electrons  in  a  metal  are  free  to  conduct  there 
is  no  restoring  force  in  the  above  equation. 

If  we  assume  am  single  wavelength  interacting  at  frequency,©  we  have 

E(t)  =  E0e-ie“ . 


Accordingly  the  position  of  the  electron  will  follow  the  field  occultation  and  can  be  written  as, 

x(t)  =  xQe~t(d* 

Substituting  into  the  first  equation  and  deviling  both  sides  by  e'mi ,  we  have 

2  to  „ 

-me<o  xQ~me  — x0  =  -eE0 


x0 


ex  ^  cox  -  i 

E0  2  2  7 
meco  co  x  + 1 


This  is  a  solution  to  the  equation  of  motion  differential  equation  for  which  a  general  oscillator  model  can  be  derived. 

,  x  4nne2x  cox  -i 

e(co)  =s,  ~ 2^j"—  •_ 

The  term  is  referred  to  as  the  background  dielectric  function  and  is  due  to  the  displacement  of  bound  electrons.  This 
formula  is  known  as  the  Drude  oscillator  model  dielectric  function  and  can  also  be  written  as 
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e(©)=eK 


Anne  x 


2 


me  <d2t2+1  mm  coV  +1 


1  *  .  4nne2T  1 

+ 1  - 


Now  we  can  go  about  solving  for  the  coefficient  of  absorption  in  terms  of  dielectric  function.  This  relationship  is  derived  in 
many  texts  an  is 


a(co)  = 


co 


cn(m ) 


and 


”(g))= vKe,+^)2+(^'^1)’ 

Accomplishing  the  appropriate  substitutions  we  have 

V2co  6  .._V2o>Wt  1 


a(co)  = 


mm 


m  2t  2  + 1 


1 


Anne % 


2 


ma 


coV  4-1 


4- « 


,  47tfte2T2  1  .y,r  „  Anne1  X  1 


v 


CO  \ 2  +  1, 


meco  cd  2t  2  + 1 ) 


In  this  form  it  seems  rather  formidable,  however  we  can  look  at  it  a  little  more  simply,  or  at  least  how  we  may  go 
about  solving  for  the  unknowns.  There  are  five  unknown  constants  in  the  system,  ea)',  eM",  x,  me  andn.  Consequently, 
with  unknown  temperature  there  will  be  six  unknowns  and  a  six  color  pyrometer  will  be  required.  In  addition  the  system  of  six 
equations  will  be  nonlinear  and  require  a  nonlinear  optimization  technique  such  as  proposed  by  Levenberg  and  Marquard. 
These  types  of  algorithms  tent  to  take  a  long  time  to  optimize  and  often  get  caught  in  local  minim’s.  There  is,  however,  some 
simplifications  that  can  be  done  for  our  application.  Since  polished  metals  have  very  few  bound  electrons  that  interact  in  the  3 
to  5pm  energy  band  we  can  assume  that  is  small,  therefore 


v/2co4jtne2i 


a(co)  =- 


mm 


m  \ 2  + 1 


47t?7e  T 


2_  2 


mc  co  2t  2  + 1 


■  +  . 


or  by  rearranging  terms  can  be  written  as 


4nne2x 2 

1  ) 

™e 

co2t2  +1> 

8nne2 

1 

4nne2x 


V  me&  co2t2  +1J 


a(co)  = 


m,  m  2t  2  + 1 


1  + 
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Conductivity  measurements  have  shown  that  the  relaxation  time  for  metals  is  in  the  range  of  10  13  to  10  15  seconds6.  Since 

1 

is  small  and  with  using  the  binomial  expansion 


our  sensor  is  operating  at  3  to  5  pm  thenT  (D  >  36  when 
approximation  we  have, 


^  2^  2 
CO  T 


oc(cg)  = 


16nne2  od2t2 
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m„ 


co  2x 2  +1 


or 


a(oo)  = . 


1 1 6nne 2 


m„ 


1  + 


2  2 
CO  X 


While  we  have  simplified  the  equation  it  is  still  nonlinear  and  suffers  from  the  same  problem  of  the  previous  equation.  We  can 
expand  this  equation  in  a  binomial  series  expansion  and  absorption  can  be  written  as. 


\6nne2  1  31 

me  _  2co  2x 2  8  co  4x 4 

or. 


a(co) 


1 1671772  2 
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1- 
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8tc  2c2t2 


X1 
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1 


128  ti4c4t  4 
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The  third  term  is  very  small  and  can  be  left  out.  Therefore  a  second  order  polynomial  should  provide  a  good  model  for  our 
calculations  in  temperature  emissivity  decoupling.  As  you  can  see  there  is  no  X  term  in  the  equation  and  experimental  results 
in  the  past  have  shown  that  there  is  a  small  amount  of  dependence  on  the  first  order  term.  However,  the  third  order  terms  and 
higher  have  been  essentially  nonexistent  in  the  data  previously  collected.  While  I  these  results  will  again  be  encountered  in 
this  research  I  will  later  address  the  impact  on  accuracy  from  the  neglecting  of  these  terms. 

5.  EXPERIMENTAL  RESULTS 


To  validate  the  emissivity  model  proposed  in  this  paper  I  set  up  an  experiment  using  a  heated  plate  of  aluminum  and  copper 
(T=500°K)  and  a  fourier  transform  IR  spectrometer  manufactured  by  Boman.  Let  me  first  start  with  a  discussion  of  the 
technique  used  for  the  solution  of  the  system  of  equations.  The  emissivity  can  be  expanded  to  as  many  orders  as  we  have 
spectral  bins  in  our  spectrometer.  However,  due  to  computational  rounding  errors  the  calculation  of  coefficients  beyond  the 
seventh  order  are  erroneous.  This  should  work  fine  for  the  validation  here  because  the  series  should  converge  after  the  second 
term.  The  experiment  calculated  the  eight  coefficients  for  the  emiossivity  model  for  the  aluminum  and  copper  samples, 
however,  to  simplify  the  discussion  I  will  use  a  second  order  model.  The  four  equations  of  the  system  are 


®PUeifyi,T)=  | +^A+.a3), Re^ yk _  for i=l,  2, 3  & 4 


X, 


ex,T  -l 
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where: 

cx  -  const. 
c2  =  const. 

^^Temperature  of  target 
X  =  Wavelength  in  microns 
Re(A,)=  Response  of  detector/filter 
dk  =  Differential  wavelength 

Rearranging  the  terms  we  have 


®P<Xeifyi>T)=al 
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Since  the  temperature  is  held  constant  then  any  three  of  the  four  equations  has  a  linear  solution.  The  first  three  equations  can 
be  written 


®,/„! (x,,r)=o,  Jci(r)Re^',^f''1  +a,  | 


e.OT«?^+ 


a3  JC3  (  T) 


Re(A,.)/A, 


A; 


Each  of  the  integrals  is  evaluatible  and  can  be  written  as 

®  J)=  Vi,  CD  +  V2,  (r)  +  a,/,,  (T) 

or  can  be  rewritten  in  linear  algebra  form  as 


®x. 

= 

h,_ 

CD  /,,(D  C,(D 
C,(D  4(D  /,,(D 
4(D  4  CD  I  AT) 


i 

• 

*2 

_ i 

The  solution  for  the  unknown  a1?  a2  &  a3  coefficients  is 


V 

/,.(D  4  CD  4( T) 

-1 

X" 

a2 

= 

4<D  4(D  4(D 

• 

a3_ 

4  CD  4(D  4(D 

.^3. 

The  procedure  of  inverting  a  matrix  is  fairly  complicated.  There  many  techniques  for  inverting  a  matrix  and 
there  usage  depends  on  the  specific  applicationT^Thavechosen  the  Gaussian  elimination  method  because- if  provides  the 
greatest  versatility  and  it  is  easily  implemented.  The  data  from  the  spectrometer  was  substituted  into  the  above  equation  the 
solution  matrix  was  attained.  What  we  saw  was  that  coefficients  after  the  square  term  were  negligible.  The  inclusion  of  the 
other  terms  constitute  a  ±7%  variation  in  the  calculated  emissivity.  Keeping  the  last  three  terms  and  calculating  the 
temperature  from  a  ninth  equation  I  arrived  at  a  575°K  prediction  and  seventy  five  degrees  off  the  true  temperature. 
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6.  CONCLUSIONS 


The  fundamental  process  of  passively  predicting  the  temperature  of  heated  materials  inevitably  will  be  limited  by  the 
accuracy  of  the  model  for  emissivity  and  instrumentation  noise.  The  theory  has  shown  that  metal  materials  have  an  emissivity 
that  follow  a  second  order  polynomial  with  an  absent  linear  term.  However,  the  experimental  data  has  shown  that  the  linear 
term  is  significant  and  can’t  be  left  out  of  the  model.  It  seems  likely  that  a  four  color  pyrometer  will  be  able  to  collect  all  the 
basic  data  needed  to  estimate  temperature  accurately.  The  source  of  the  error  in  temperature  prediction  of  the  experiment 
performed  for  this  paper  is  primarily  due  to  instrumentation  calibration  inaccuracies.  To  eliminate  this  source  of  error  new 
instrumentation  is  required  to  be  developed.  This  will  be  the  subject  of  further  work. 
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Introduction: 

This  report  documents  the  research  methods  used  and  data  compiled  during  this  project. 
The  development  of  algorithms  to  extract  spectral  emissivity  and  temperature  values  from 
radiometric  data  taken  in  several  bands  in  the  3-5  micron  region  is  described.  In  addition, 
progress  toward  an  instrument  capable  of  providing  the  raw  data  for  these  algorithms  for 
each  pixel  in  a  single  snapshot  image  is  detailed. 


I.  Tasks: 


A)  Search  available  literature  for  relevant  spectral  emissivity  studies. 

The  literature  search  was  completed.  Our  findings  indicated  that  the  spectral  emissivity  of 
many  metals  in  the  3-5  micron  region  could  be  accurately  represented  by  a  third  order 
polynomial  in  wavelength: 


e(X)  -a-t-  b-X  +-  c-X 


2 


As  shown  in  Fig.  1,  Aluminum,  and  Steel  Alloy  show  a  parabolic  dependency  of 
emissivity  with  wavelength .  Ceramics  and  dielectric  materials  however,  exhibit  more 
complicated  spectral  behavior  as  shown  for  brick  and  fireclay  in  Fig.  2.  Since  composite 
armor  can  produce  projectile  debris  consisting  of  ceramic  and  metal,  a  higher  order 
polynomial  approximation  of  emissivity  must  be  used.  In  addition,  a  higher  order 
expression  allows  the  identification  of  emission  lines  from  combustion  processes  which 
could  perturb  the  temperature  solution.  To  do  otherwise  would  be  to  adapt  the  problem  to 
an  existing  technology  when  competent  engineering  practice  demands  the  opposite 
approach. 


B)  Determine  the  feasibility  of  using  spectral  radiance  data  to  find  absolute 
temperature  without  prior  knowledge  of  spectral  emissivity. 

Originally,  a  polynomial  fit  of  spectral  emissivity  was  combined  with  the 
Levenberg-Marquardt  algorithm1  to  solve  a  system  of  equations  of  the  form: 
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Where  Kcal  is  a  known  radiometric  calibration  constant  and  Data;  is  the  photon 
irradiance  ( photons/cmA2-second  )  in  each  band.  The  integrals  on  the  left-hand  side 


Fig.  1  Spectral  emissivity  of  Aluminum  (Top)  and  Steel  Alloy2. 


Fig.  2  Spectral  emissivity  of  brick  (Top)  and  fireclay  100  -  300  C  2. 


represent  the  calculated  band  inadiance,  which  is  ideally  equal  to  the  measured  band 
irradiance  on  the  right.  For  a  third  order  polynomial  expression  of  emissivity,  at  least 
four  such  equations  must  be  solved  simultaneously  for  a,  b,  c  and  T. 

In  practice,  the  system  can  be  solved  only  approximately.  It  was  found  during 
tests  of  the  algorithm  on  simulated  data  with  additive  Gaussian  noise  that  higher  order 
polynomial  expressions  for  e(X)  resulted  in  unacceptable  susceptibility  to  noise.  Since 
we  were  aiming  for  field  images  of  projectile  debris  at  small  integration  times,  we 
anticipated  noisy  data.  Therefore  the  polynomial  fit  approach  was  abandoned  in  favor  of 
independent  band  averaged  emissivities  for  each  band.  This  resulted  in  a  system  of 
equations  of  the  following  form: 
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Effectively,  we  chose  to  model  e(X)  by  a  rectangle  function  in  each  band  and  the  entire 
spectral  profile  by  several  such  bands  as  shown  in  Fig.  3.  The  astute  reader  will  note  that 
this  system  is  inherently  underdetermined  due  to  the  presence  of  one  more  variable  ( T  ) 
than  the  number  of  equations.  Accordingly,  additional  constraints  such  as  positivity  on 
si,  e2  etc.  had  to  be  applied  and  several  different  initial  values  of  temperature  had  to  be 
used  as  starting  points.  Further  computer  simulations  indicated  that  Levenberg- 
Marquardt,  applied  to  such  a  system,  tended  to  produce  inaccurate  results  owing  to  the 
algorithm’s  tendency  to  converge  to  local,  rather  than  absolute,  minima.  Accordingly, 
globally  convergent  algorithms  were  investigated  and  tested.  After  considerable  trial  and 
error,  Broyden’s  Method  proved  to  give  good  noise  immunity  and  reasonably  accurate 
results  when  applied  to  simulated  data  with  noise3.  As  with  Levenberg-Marquardt, 
several  different  seed  values  of  T  over  the  interval  of  interest  had  to  used,  then  the 
algorithm  automatically  chose  the  best  solution  based  on  a  least-squares  merit  function. 
An  attempt  to  fit  one  fewer  emissivity  value  than  the  number  of  data  points,  thereby 
producing  an  adequately  determined  system,  proved  unsuccessful. 

Having  obtained  good  results  with  this  technique  on  simulated  noisy  data,  we 
tested  the  program  on  real  data  generated  with  the  experimental  setup  shown  in  Fig.  4. 

We  simulated  a  target  of  unknown  emissivity  and  temperature  by  placing  a  known 
spectral  filter  in  the  radiation  path  from  a  blackbody  set  at  a  known  temperature.  Since 
we  had  not  yet  manufactured  an  IR  disperser,  the  output  of  the  laboratory  blackbody  was 
passed  through  a  monochromator  to  produce  radiation  in  six  contiguous  bands  with  center 
wavelengths  of 3.20, 3.44, 3.68, 3.92 , 4.16  and  4.40  microns.  The  monochromator  output 
was  first  chopped  and  then  sent  to  a  256  x  256  pixel  PtSi  camera  operating  at  78  K. 

This  was  the  imager  we  intended  to  use  in  the  actual  IR  snapshot  imaging  spectrometer. 

A  magnified  image  of  the  exit  slit  of  the  monochromator  was  sent  via  RS-170  to  a  video 
monitor  and  subsequently  to  a  lock-in  amplifier  to  extract  the  signal  magnitude.  Since  we 
were  unable  to  capture  video  frames  from  the  camera  with  frame  grabbers  currently  in  our 
lab,  the  use  of  a  lock-in  was  necessary  to  separate  die  signal  of  interest  from  the  various 
components  of  the  PtSi  imager’s  composite  video  output  4 
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Fig.  3  Rectangle  approximation  of  spectral  emissivity. 


Fig.  4  Experimental  setup  for  unknown  target  measurements. 


Readings  were  taken  in  each  band  with  die  blackbody  set  for  873  K  to  obtain  a 
calibration  of  the  system  as  a  whole.  Then  the  blackbody  temperature  was  increased  to 
973  K ;  a  previously  scanned  bandpass  filter  was  placed  in  the  radiation  path  after  the 
chopper  and  the  data  collection  process  was  repeated.  The  filter  transmission  versus 
wavelength  was  used  to  simulate  a  target  with  unknown  spectral  emissivity  characteristics 
at  an  unknown  temperature.  Both  the  system  calibration  data  and  the  target  data  were 
processed  with  a  program  that  incorporated  Broyden’s  Method.  The  results  are  shown  in 
Fig.  5  below. 


FILTER  TRANSMISSION  VERSUS  WAVELENGTH  IN  MICRONS 
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973  K  Actual  Temp. 

(B<mm>)0  =  962.461  Solved  Temp. 

Fig.  5  Results  of  the  emissivity-temperature  algorithm.  The  bandpass 
filter  transmission  is  the  continuous  dashed  line;  the  measured  data  for  the  six 
spectral  bands  is  the  upper  trace  with  circles;  the  output  of  the  program  is  the 
lower  trace  with  squares. 


The  algorithm  converged  to  a  temperature  of  962.5  K  which  is  10.5  K  lower  than  the 
actual  source  temperature,  thereby  proving  the  ability  to  recover  temperature  with  good 
accuracy.  The  lack  of  a  good  emissivity  fit  was  probably  due  to  inadequate  wavelength 
calibration  of  the  monochromator. 


C)  Design  and  manufacture  diffractive  elements  for  the  3  to  5  micron  region. 


This  task  has  been  completed.  Prior  to  the  design  and  manufacture  of  dispersive 
elements,  it  was  strongly  suggested,  by  the  contract  monitor,  that  we  spend  several 
months  developing  a  model  of  the  IR  Computed  Tomography  Imaging  Spectrometer.  Our 
extensive  experience  with  the  visible  version  of  this  instrument  was  apparently  not 
applicable  to  the  3  -  5  micron  band.  We  then  used  this  model  as  a  starting  point  to 
optimize  the  disperser  design  while  recognizing  that  trial  and  error  experiments  were 
inevitable.  The  optical  layout  upon  which  this  model  is  based  appears  below. 


OPTICAL  SYSTEM  LAYOUT 
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Fig.  6.  Simplified  optical  layout  of  the  IR  CTIS  system  showing 
the  focal  lengths  of  the  objective  lens,  the  collimating  lens  and  the  final 
imaging  lens  which  is  actually  attached  to  the  PtSi  camera  dewar.  The 
field  stop  (F.S.)  is  located  at  the  focal  plane  of  both  the  objective  lens 
and  the  collimating  lens  to  reduce  stray  light. 


A  Fourier  optics  approach  was  chosen  for  the  system  model  to  account  for  the 
diffractive  effects  of  the  grating.  The  method  is  an  extension  of  a  incoherent  diffraction 
calculation  presented  by  Goodman  to  two-dimensional  binary  phase  gratings  .  A  two 
dimensional  binary  phase  grating  is  one  which  is  etched  into  a  transmissive  material  and 
has  only  two  levels  as  shown  schematically  below  in  Fig.  7. 
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Fig.  7.  Two-dimensional  binary  phase  grating  in  cross  section.  The 
grating  period  is  given  by  A  and  the  grating  modulation  thickness 
is  given  by  t.  From  above,  the  grating  looks  like  a  pattern  of 
square  islands  with  the  “moats”  separating  them  as  wide  as  the  squares 
themselves.  This  particular  pattern  was  chosen  because  it  could  be 
fabricated  at  the  Optical  Sciences  Center  in  GaAs  and  required  only 
a  single  mask  and  etch  step. 


Following  and  generalizing  Goodman’s  model,  we  found  for  a  one-dimensional 
sinusoidal  phase  grating  coincident  with  a  square  aperture, 
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where 

t(xiyi)  =  The  pupil  transmission  function  of  the  grating. 

m(/t)  =  2 a{n - 1  )t!  X  is  the  phase  modulation  factor,  n  =  3.4  for  GaAs. 

fo  =  1  /  A  is  the  grating  frequency  in  line  pairs  per  mm. 

1  =  The  lateral  dimensions  of  the  square  aperture. 

and  rect(xi/l)rect(yi/l)  is  a  function  which  is  1  inside  the  aperture  and  0  outside. 


At  this  point,  we  chose  only  to  model  a  one  dimensional  grating  coincident  with  a 
square  aperture  and  calculate  the  resulting  diffraction  pattern  a  distance  z  away.  Effects 
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of  the  other  optical  components  in  the  IR  CHS  would  be  included  later.  To  apply  this 
result  to  a  binary  squarewave  phase  grating,  we  expanded  the  squarewave  phase  profile 
in  a  Fourier  series  to  obtain: 
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for  the  pupil  transmittance  function  of  a  one-dimensional  binary  squarewave  phase 
grating.  The  grating  is  a  squarewave  along  X|  and  constant  along  yt.  To  find  the 
Fraunhofer  diffraction  pattern  a  distance  z  away,  we  took  the  Fourier  Transform  of  this 
expression,  noting  that  multiplications  become  convolutions  in  the  transform  domain. 


155 


F  t 


Vyt 


00 

V 


.  \  x  ^  0  y/ 


q  = 


-0- 


00 

V 


-5 'f 
3-Jt  ;  x 


•0 


.  q  = 


L2sinc  ^L-fJ-sinc  (L-fJ 


(6) 

where  we  truncated  the  expansion  of  the  sine  function  to  four  terms  for  simplicity.  For  a 
unit  amplitude  plane  wave  incident  upon  the  aperture,  this  transmission  function  is  also 
the  amplitude  of  the  transmitted  wave  immediately  behind  the  aperture. 
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We  took  advantage  of  the  properties  of  convolution  applied  to  delta  functions,  namely  if 
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and  arrived  at  the  following  result. 
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Here  the  asterisk  denotes  multiplication. 

The  Fraunhofer  diffraction  pattern  at  a  distance  z  from  the  aperture  was  found  by  taking 
the  modulus  squared  of  this  expression  and  scaling  the  appropriate  variables.  The 
coordinate  system  in  the  Fraunhofer  plane  is  denoted  by  xo,  yo. 
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This  result  was  applied  to  the  situation  of  a  unit  amplitude  plane  wave  incident  on 
a  two-dimensional  binary  phase  grating  coincident  with  a  lens  of  focal  length  B  to 
calculate  the  diffraction  pattern  at  the  focal  plane. 
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Equation  12  describes  the  idealized  case  of  a  unity  amplitude  plane  wave  incident 
on  a  grating-lens  system  similar  to  the  grating-imaging  lens  combination  in  the  IR  CTIS. 
By  considering  the  diffraction  properties  of  the  other  lenses  in  the  system,  we  were  able 
to  modify  (12)  to  characterize  the  whole  system.  To  do  this,  we  considered  the  case  of  an 
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extended  object  which  produces  a  uniform  irradiance  distribution  in  the  field  stop  of  Fig. 
6.  We  noted  that  an  ideal  objective  lens  would  merely  scale  the  object  exitance  by 
l/(4(F/#)2+l)  and  we  asked  about  the  nature  of  the  wavefront  incident  upon  the  grating¬ 
imaging  lens  combination  just  in  front  of  the  PtSi  FPA.  We  were  able  to  show  that  in 
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collimated  space  just  prior  to  the  grating,  the  plane  wavefront  amplitude  is 
given  by 


where 

F  num  =  Objective  lens  F/#. 

I0  =  Generalized  “Intensity”  leaving  the  object. 

fj  =  focal  length  of  the  collimating  lens. 

D  =  Lateral  dimension  of  the  square  field  stop. 

To  simplify  the  derivation,  we  took  into  account  the  finite  extent  of  the  field  stop 
only,  not  the  objective  or  collimating  lenses.  To  further  simplify  the  derivation  from  this 
point  on,  we  chose  not  include  the  finite  extent  of  the  grating-imaging  lens  combination. 
These  approximations  result  in  the  following  expression  for  the  focal  plane  irradiance  due 
to  an  extended  object  of  exitance  Mobj« 
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Here, 

m  =  m(X)  =  2n{n-X)tl  X 
A  =  The  grating  period. 


Since  the  measurable  quantity  we  sought  was  the  number  of  electrons  collected 
from  each  pixel  of  the  PtSi  camera,  we  needed  to  take  into  account  the  spectral  response 
of  the  FPA  and  the  300  K  background  noise.  Fortunately,  both  quantities  had  been 
previously  measured  for  this  particular  camera5.  The  spectral  quantum  efficiency  of  the 
array  can  be  characterized  by5: 
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with  Cl  =  quantum  yield  =  0.190  ev'1  and  v|/  =  barrier  height  =  0.240  ev  by  previous 
measurement6.  This  gave  for  the  number  of  electrons  produced  by  a  pixel  during  one 
integration  period  when  the  system  was  viewing  an  object  of  temperature  T  after  a  300  K 
background  offset  correction: 
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where 

Topt  =  System  optics  transmission, 
tint —  FPA  integration  time. 

Apucd  =  Pixel  area. 

A.2  =  Longest  wavelength  of  interest 
A.1  =  Shortest  wavelength  of  interest. 

and  the  973  K  object  exitance  (photons/cm2-sec.)  with  the  300  K  background  subtracted 
in  equation  15  was  given  by 


159 


2-it  -c 


4  /  h-c 

•exp - ; 

\  U-k-T/ 


1: ; 

‘  1=973 


2-it-c 


.4  h-c 
A,  -  exp 


X-k-T/ 


li 


1=300 


(18) 


For  tractable  computer  implementation,  the  summations  from  -qo  to  qo  in  (15)  were 
reduced  to  summations  from  -2  to  2.  In  other  words. 
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The  model  was  implemented  in  MathCad  and  run  for  various  combinations  of  etch  depth 
( t  in  Fig.  7)  and  wavelengths  of  interest  in  the  3  to  5  pm  band  using  a  constant  grating 
period  of  A  =  1/1 1  mm.  A  run  for  an  etch  depth  of  0.85  pm  with  XI  =  3.4  pm  and  X2  = 
3.6  pm  and  a  blackbody  temperature  of  T=  973  K  is  shown  in  Fig.  8. 
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Fig.  8.  Simulation  run  Of  IR  CTIS  physical  system  including  PtSi  quantum 
efficiency.  Shown  in  this  figure  is  the  number  of  electrons  collected  per  pixel 
over  the  surface  of  the  256x256  array  during  one  integration  time.  Pixels  have 
been  burned  into  8x8  groups  to  increase  computational  speed. 
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Our  goal  was  to  achieve  uniformity  of  response  across  the  surface  of  the  PtSi 
FPA  and  to  compensate  for  the  characteristic  falloff  of  quantum  efficiency  at  longer 
wavelengths  in  this  type  of  photodiode  array.  Numerous  simulations  were  run  for  etch 
depths  between  0.3  pm  and  1.0  pm  for  various  wavelength  intervals  between  3.4  pm  and 
4.9  pm.  Based  on  these  simulation  runs,  it  was  determined  that  an  etch  depth  near  0.7  pm 
was  optimal.  Due  to  the  expense  and  delays  associated  with  outside  vendors,  we 
eventually  chose  to  develop  grating  fabrication  capability  in-house.  A  series  of  small 
prototype  gratings  were  fabricated  on  optical  quality  Gallium  Arsenide  wafers  using  a 
reactive  ion  etch  process.  Depths  ranging  from  0.6  pm  to  0.9  pm  were  chosen  and 
evaluated  in  the  breadboard  prototype.  These  results  will  be  discussed  in  the  next  section. 

Initially,  problems  were  encountered  with  maintaining  a  uniform  etch  depth  across 
a  sufficiently  large  area  (25mm  X  25  mm)  to  be  useful  in  our  instrument.  These 
difficulties  were  surmounted  by  increasing  the  flow  rate  of  the  methane  and  chlorine  gas 
mixture  used  in  the  etching  chamber  and  a  grating  with  a  1  micron  nominal  etch  depth 
was  fabricated.  A  surface  profile  scan  of  this  grating  was  taken  with  a  Wyko 
interferometer  and  the  results  appear  in  Fig.  9. 
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Fig.  9.  Wyko  interferometer  surface  profile  of  the  large  area  crossed 
binary  phase  grating  fabricated  using  Reactive  Ion  Etching  at  the  Optical 
Sciences  Center  by  Michael  Descour  and  Daniel  Simon.  Average  etch 
depth  came  in  at  1.1  pm.  (Courtesy  M.  Descour  and  D.  Simon). 


161 


r 


D)  Breadboard  a  prototype  snapshot  imaging  spectrometer. 

Having  been  graciously  warned  that  the  AFOSR  does  not  like  hardware  built 
under  6.1  funding,  we  nevertheless  proceeded  to  breadboard  a  prototype  snapshot 
imaging  infrared  spectrometer  system,  as  specified  in  our  original  contract,  without 
hesitation.  Clearly,  any  perceived  agenda  underlying  such  an  inconsistency  could  only 
have  been  the  product  of  an  overactive  imagination.  During  this  phase  of  the  project  we 
gained  many  optical  insights  into  the  system  shown  below  in  Fig.  10. 


Fig.  10.  Protoytpe  snapshot  IR  imaging  spectrometer.  In  the  foreground,  a 
commercial  blackbody  (white  cylindrical  object )  supplies  IR  radiation  to  a 
Perkin-Elmer  monochromator.  The  exit  slit  of  the  of  the  monochromator  is 
reflected  by  a  folding  mirror  into  the  spectrometer  optical  train  in  the  right 
background.  From  left  to  right  we  have  1)  a  MOO  mm  objective  lens  followed 
by  2)  a  2  mm  square  field  stop.  The  black  conical  object  with  it’s  apex  to  the  left 
is  3)  a  100  mm  collimating  lens  followed  by  the  disperser  mount. 

Finally,  on  the  far  left  is  the  50  mm  re-imaging  lens  attached  to  the  liquid 
nitrogen  dewar  which  houses  the  PtSi  array.  See  Fig.  6.  for  further  details.  Not 
shown  is  the  computer  and  associated  frame-grabber  card  used  to  digitize  the 
image. 


Using  this  setup,  we  evaluated  small  prototype  gratings  of  varying  etch  depths  to 
determine  the  optimum  depth  as  mentioned  in  section  C)  above.  Once  we  eventually 
found  a  way  for  a  commercial  frame-grabber  board  to  reliably  accept  the  non-standard 
vidio  output  from  the  PtSi  camera,  we  were  able  to  acquire  foil  images.  The 
monochromator  was  used  to  restrict  the  spectral  bandwidth  so  that  the  diffraction 
efficiency  as  a  function  of  wavelength  could  be  gauged.  The  result  for  one  grating 
appears  in  Fig.  11. 
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Fig.  11.  Grayleve!  profile  of  a  0.85  pm  deep  grating  with  a  period  of  A=  1/1 1  mm 
at  a  center  wavelength  of  =  3.5  pm.  The  vertical  axis  is  proportional  to  the 
number  of  electrons  collected  by  each  pixel  of  the  focal  plane  array  during  one 
one  integration  period.  The  horizontal  axes  are  proportional  to  the  x  and  y  axis  of 
the  focal  plane.  The  blackbody  temperature  was  973  K. 


A  comparison  with  Fig.  8  reveals  significant  discrepancies  with  the  mathematical 
model.  In  particular,  the  center  order  in  the  real  data  shown  above  is  far  greater  than 
predicted  by  the  model.  We  believe  this  is  due  the  requirement  of  having  to  truncate  the 
quadruple  summation  in  (19)  for  tractable  (12  hrs.  on  a  266  MHz  Pentium  II)  computer 
simulation.  However,  photons  from  higher  order  summation  terms  contribute  directly  to 
the  zero  and  first  order  diffraction  patterns  in  a  complex  fashion  as  indicated  in  (15). 
More  work  needs  to  be  done  in  this  area  if  modeling  is  determined  to  be  actually 
important  But  as  far  as  disperser  design  is  concerned,  we  believe  we  have  a  workable,  if 
not  optimal,  large  area  candidate. 
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E)  Measure  high  temperature  laboratory  targets  which  approximate  actual 
penetrator  debris. 

Actual  behind  panel  debris  from  Eglin  AFB  (Tests  at  Range  C-64A)  would  have 
been  analyzed,  had  it  been  made  available. 


F)  Extract  temperature  and  emissivity  values  to  compare  with  true  parameters. 

This  task  has  been  completed  for  the  purposes  of  software  validation.  However  a 
serious  IR  camera  malfunction  prevented  us  from  completing  this  phase  of  the  project 
with  an  actual  snapshot  imaging  spectrometer.  To  date  we  have  been  unable  to  fix  the 
problem. 


III.  Publications:  None. 


IV.  Technical  Exchange  Meetings: 

Telephone  calls  and  in-person  meetings  with  the  technical  contract  monitor  plus 
several  visits  by  the  technical  monitor  to  Tucson. 


V.  Travel:  None. 


VI.  Problem  Areas: 

A  clock  driver  coaxial  cable  in  the  IR  camera  was  damaged  during  routine 
maintenance.  All  attempts  at  repair  have  failed  to  date.  The  camera  probably  requires 
major  rewiring  inside  the  dewar,  but  hopefully  not  a  new  PtSi  array. 
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Behind  Panel  Debris 
Temperature  Profilometry 


Objective 

The  objective  of  this  research  is  to  extend  the  techniques  developed  during  the 
first  two  years  of  this  project.  Namely,  we  will  improve  the  emissivity-temperature 
recovery  algorithm  already  developed  and  we  will  improve  the  instrumentation  that 
provides  the  raw  data  to  this  algorithm.  As  in  the  first  phase  of  this  research,  the 
fundamental  goal  is  to  investigate  ways  of  providing  an  absolute  temperature  map  of 
behind  panel  debris  from  penetrating  type  weapons  and  to  understand  the  radiative 
process  of  the  debris  itself. 


Background 

Under  the  present  contract,  a  unique  iterative  approach  was  shown  to  obtain 
accurate  temperature  of  a  remote  object  with  a  complex  emissivity  profile  from 
measurements  of  a  test  object’s  radiant  exitance  in  six  contiguous  spectral  sub-bands 
covering  the  3.4  to  4.9  micron  region.  This  method  does  not  expand  the  object’s 
emissivity  profile  in  a  polynomial  function  of  wavelength  such  as 

£(X)  ■—  a  ■+•  bX  +■  cX  "t"  dX 

which  is  valid  over  the  entire  spectral  region.  Instead  our  method  expresses  emissivity  as 
an  independent  constant  within  each  spectral  sub-band.  The  emissivity  profile  for  the 
entire  spectral  region  then  is  simply  the  juxtaposition  of  all  the  sub-band  emissivities. 
Simulations  of  noisy  data  have  shown  this  technique  to  be  effective  in  recovering  target 
temperatures  even  when  the  data  signal-to-noise  ratio  is  low.  An  experimental  setup 
consisting  of  a  hot  target,  a  monochromator  and  a  PtSi  camera  operating  in  the  3.4  -  4.9 
micron  region  was  used  to  simulate  an  infrared  imaging  spectrometer1.  The  radiation 
from  the  target  was  separated  in  to  six  bands  with  center  wavelengths  of  3.20, 3.44, 3.68, 
3.92, 4.16  and  4.40  microns.  In  Fig.  1,  the  continuous  dashed  curve  shows  the  target’s 
true  emissivity  as  measured  with  a  high  resolution  non-imaging  spectrometer.  The  broad 
upper  trace  with  circles  is  the  simulated  imaging  spectrometer  output  with  the  instrument 
spectral  response  divided  out.  The  lower  broad  trace  with  squares  is  the  fit  given  by  our 
iterative  algorithm  which  converged  to  within  1 1  K  of  the  actual  target  temperature  of 
973  K,  despite  less  than  ideal  recovery  of  the  true  emissivity  profile. 

In  a  parallel  effort,  an  imaging  spectrometer  using  computed  tomography  was 
developed  for  the  visible  spectrum  under  another  contract  The  software  will  be  readily 
applicable,  with  minor  modifications,  to  the  infrared  version  of  this  instrument  which  is 
currently  under  construction.  This  will  enable  us  to  apply  the  temperature  recovery 
algorithm  on  a  pixel  by  pixel  basis. 
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FILTER  TRANSMISSION  VERSUS 
WAVELENGTH  IN  MICRONS 


Actual  Temperature  =  973  K 
Solved  Temperature  =  962.5  K 


Fig.  1  Results  of  the  emissivity-temperature  algorithm.  The  true  target 
emissivity  is  the  continuous  dashed  line;  the  measured  data  for  the  six  spectral 
bands  is  the  broad  upper  trace  with  circles;  the  output  of  the  program  is  the 
broad  lower  trace  with  squares. 


168 


Approach 


Our  approach  to  the  next  phase  of  this  research  will  focus  on  the  fundamental 
physics  of  radiating  surfaces,  the  data  acquisition  system  as  a  whole,  the  infrared 
disperser  in  particular  and  further  evolution  of  the  temperature  recovery  algorithm.  The 

block  diagram  of  the  overall  system  in  Fig.  2  indicates  the  interrelationship  of  these 
components. 
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Fig.  2  Overall  system  block  diagram 

In  the  second  phase  of  the  research,  we  will  continue  to  pursue  increased 
understanding  of  both  the  measurement  technology  and  the  physics  of  IR  radiation 
production.  We  will  investigate  the  limitations  of  the  Planck  model  itself  that  can 
influence  accuracy.  Studies  in  the  400  nm  to  800  nm  wavelength  range  of  behind  armor 
debris  indicate  that  emission  lines  from  reacting  Fe02  can  perturb  the  otherwise  T  = 
1050  K  greybody  spectrum  of  the  debris  2.  If  similar  emission  lines  exist  in  the  3  -  5 
micron  region,  they  could  seriously  degrade  the  accuracy  of  the  temperature  recovery 
algorithm.  We  need  to  identify  potentially  reactive  species  and  find  ways  of  -• 
incorporating  them  into  the  model. 

The  data  acquisition  system  as  a  whole  can  benefit  from  improved  spatial  and 
spectral  resolution.  For  a  given  focal  plane  format,  both  spatial  and  spectral  resolution 
are  primarily  determined  by  the  tomographic  reconstruction  software.  By  modifying  this 
software  we  can  will  greater  spatial  and  spectral  resolution  than  the  current  goal  of  10  x 
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10  reconstructed  pixels  with  10  spectral  bands  p*er  pixel.  In  addition,  the  overall 
throughput  of  the  system  can  be  improved  for  true  IR  snapshot  imagery  to  be  feasible  . 
Modeling  of  the  present  optical  system  including  the  current  PtSi  camera  indicates  that 
low  SNR  will  occur  at  wavelengths  beyond  4.5  microns  making  integration  times  less 
than  1/60  of  a  second  a  practical  impossibility.  Since  integration  times  below  100 
microseconds  are  needed  to  freeze  rapidly  moving  behind  panel  debris,  increased  optical 
throughput  and  a  higher  quantum  efficiency  FPA  are  required. 

For  similar  reasons,  more  efficient  holographic  dispersers  must  be  designed  and 
tested.  The  current  system  will  employ  a  simple  binary  crossed  phase  grating  which  will 
not  make  efficient  use  of  the  available  IR  radiation  from  the  target.  In  conjunction  with 
design  and  fabrication  techniques  being  developed  in  parallel  for  the  visible  version  of 
this  instrument,  we  will  fabricate  better  dispersers. 

Finally,  additional  work  needs  to  be  done  on  the  emissivity-temperature  algorithm 
both  to  improve  its  performance  and  to  understand  its  limitations.  As  it  stands,  the 
software  is  able  to  take  six  sub-band  data  points  and  solve  for  six  sub-band  emissivity 
values  only  by  imposing  numerous  constraints  upon  the  solution.  While  the  technique  is 
entirely  automatic  and  needs  no  user  intervention,  it  is  slow  and  reduces  the  temperature 
recovery  accuracy.  We  will  investigate  methods  such  as  Simulated  Annealing  which  has 
proven  capable  of  dealing  with  previously  insoluble  problems  of  this  type 3. 
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DEPARTMENT  OF  THE  AIR  FORCE 

AIR  FORCE  RESEARCH  LABORATORY  (AFMC) 
EGLIN  AIR  FORCE  BASE,  FLORIDA 


MEMORANDUM  FOR  MARK  F.  HOPKINS 

430  W.  NEW  ENGLAND  AVE 
SUITE  10 

WINTER  PARK,  FL  32789 
(407)539-1639 


FROM:  AFRL/MNMF 


SUBJECT:  Patent  Award  and  Final  Review 


1 .  Congratulations  on  your  being  awarded  a  patent  for  your  outstanding  work  on  the  development  of  a 
new  and  novel  spectroscopic  imaging  method.  The  patent  number  5,708,504  for  an  Interfering  Imaging 
Spectrometer  was  issued  on  13  Jan  98.  A  copy  of  this  award  will  be  placed  in  the  program  files  prior  to 
retirement  and  forwarded  to  Civilian  personnel  for  inclusion  in  your  records.  Your  contributions  have 
been  recognized  as  key  to  future  capabilities  for  research  in  energetic  materials.  Recent  discussions  have 
indicated  interest  in  further  advancing  the  work  you  initiated  here.  Innovative  Research  investments  may 
be  available  at  some  date  in  the  future  to  support  University  Research  and  Post  Doctoral  Research  in 
development  and  validation  of  the  empirical  models  for  surface  emissivity  you  discussed  in  your  final 
report.  Future  sponsored  research  would  be  dependent  on  this  work,  based  on  input  from  key  scientists  in 
various  discipline  areas. 

2.  In  order  to  complete  and  close  out  the  prior  research  in  this  area,  it  is  anticipated  that  a  review  of  the 
reports  and  research  files  be  performed  during  the  week  of  12-18  April  at  your  facility  and  University  of 
Central  Florida.  A  research  plan  from  our  discussion  will  be  sponsored  for  briefing  to  the  MN  associate 
chief  scientist  and  to  AFOSR  (Capt.  Brian  Sanders). 

3.  Mr  Donald  R.  Snyder  is  the  point  of  contact  for  this  effort  and  can  be  reached  at  (850)  882-2005  ext. 
257  or  (850)  882-5810. 


Kenneth  L.  Echtemacht,  Jr.,  Major,  USAF 
Chief,  Fuzes  Branch 


Attachment: 

Notification  of  Patent  Award 
cc: 

CA-N  (Dr.  Sierakows.ki) 
MNOR  (Kay  Cushing) 
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WRIGHT-PATTERSON  AIR  FORCE  BASE  OHIO 


30-Dec-9  7 


MEMORANDUM  FOR  WL/MNSI  (MARK  F.  HOPKINS) 

101  W.  EGLINBLVD,  STE309 
EGLINAFB  FL  32542 

FROM:  AFMCLO/JAZ 

2240  B  Street,  Room  100 

Wright- Patterson  AFB,  OH  45433-7109 

SUBJECT:  S/N  08/736,690  (AF  Inv  D00103)  -  M.  Hopkins  -  INTERFERING 
IMAGING  SPECTROMETER 


1.  Attached  is  a  copy  of  the  Issue  Notification  in  the  matter  of  the  above-identified 
case.  As  indicated,  the  patent  issue  date  is  13  Jan  98  and  the  patent -number  is 
5,708,504. 


2. 


A  copy  of  the  issued  patent  will  be  forwarded  to  you  as  sbon  as  we  receive  it. 

/ 


BOBBY  D.  SCEARCE 


Chief,  Physics  Section 

New  Technology  Division 

Directorate  of  Intellectual  Property  Law 


Attachment: 
Issue  Notification 
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[57]  ABSTRACT 

An  interfering  imaging  spectrometer  is  described  wherein  a 
single  interference  filter  disposed  in  preselected  angular 
relationship  with  one  or  more  mirrors  separates  an  incoming 
optical  signal  or  field  of  view  into  multiple  contiguous 
images  separated  spatially  and  spectrally,  each  image  having 
a  unique  spectral  content. 

2  Claims,  3  Drawing  Sheets 
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INTERFERING  IMAGING  SPECTROMETER 
RIGHTS  OF  THE  GOVERNMENT 

The  invention  described  herein  may  be  manufactured  and 
used  by  or  for  the  Government  of  the  United  States  for  all 
governmental  purposes  without  the  payment  of  any  royalty. 

BACKGROUND  OF  THE  INVENTION 

The  present  invention  relates  generally  to  imaging  spec¬ 
trometer  systems  and  methods,  and  more  particularly  to  an 
imaging  spectrometer  system  having  no  moving  parts  and 
rapid  response  time. 

In  the  operation  of  existing  commercially  available  imag¬ 
ing  spectrometers,  spectral  scanning  is  generally  performed 
using  a  prism,  diffraction  grating  or  filter  which  is  moved 
into  and  out  of  the  optical  path  of  an  incident  optical  signal. 
The  time  required  to  scan  the  spectral  energy  of  the  optical 
signal  onto  an  optical  detector,  such  as  a  charge-coupled- 
device.  can  be  at  least  several  milliseconds  and  is  therefore 
undesirably  long.  Because  certain  spectral  analysis 
applications,  such  as  the  temperature  analysis  of  projectiles, 
require  fast  (of  the  order  of  a  microsecond)  response  time, 
prior  art  spectrometer  systems  and  methods  are  inadequate. 

The  invention  solves  or  substantially  reduces  in  critical 
importance  problems  with  prior  ait  spectrometer  systems 
and  methods  by  providing  an  improved  imaging  spectrom¬ 
eter  including  a  single  interference  filter  for  spectral  analysis 
of  the  incoming  signal  and  takes  advantage  of  the  spectral 
transmission  and  reflection  variations  with  incident  angle  of 
the  incoming  signal  within  the  interference  filter.  The  inven¬ 
tion  exhibits  a  rapid  response  time  in  the  microsecond  range, 
and  multiple  images  can  be  formed  from  a  single  incoming 
signal  so  that  a  single  snapshot  will  acquire  all  spectral 
information  of  the  signal.  The  spectrometer  of  the  invention 
has  no  moving  parts  and  can  therefore  have  a  rugged  simple 
structure  which  is  easily  and  accurately  calibrated. 

It  is  therefore  a  principal  object  of  the  invention  to  provide 
an  improved  spectrometer  system  and  method. 

It  is  a  further  object  of  the  invention  to  provide  an 
interfering  imaging  spectrometer  system  and  method. 

It  is  yet  another  object  of  the  invention  to  provide  an 
imaging  spectrometer  having  very  fast  response  time. 

It  is  another  object  of  the  invention  to  provide  an  inter¬ 
fering  imaging  spectrometer  having  no  moving  parts. 

It  is  yet  another  object  of  the  invention  to  provide  an 
interfering  imaging  spectrometer  which  is  rugged  in  con¬ 
struction  and  which  can  be  easily  and  accurately  calibrated. 

It  is  another  object  of  tile  invention  to  provide  an  inter¬ 
fering  imaging  spectrometer  having  a  single  spectral  filter. 

These  and  other  objects  of  the  invention  will  become 
apparent  as  a  detailed  description  of  representative  embodi¬ 
ments  thereof  proceeds. 

SUMMARY  OF  THE  INVENTION 

In  accordance  with  the  foregoing  principles  and  objects  of 
the  invention,  an  interfering  imaging  spectrometer  is 
described  wherein  a  single 'interference  filter  disposed  in 
preselected  angular  reiationship  with  one  or  more  mirrors 
separates  an  incoming  optical  signal  or  field  of  view  into 
multiple  contiguous  images  separated  spatially  and 
spectrally,  each  image  having  a  unique  spectral  content. 

r  9  ' 

DESCRIPTION  OF  THE  DRAWINGS 

The  invention  will  be  more  clearly  understood  from  the 
following  detailed  description  of  representative  embodi- 
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ments  thereof  read  in  conjunction  with  the  accompanying 
drawings  wherein: 

FIG.  1  is  a  schematic  block  diagram  illustrating  the 
essential  components  of  a  spectrometer  system  according  to 
5  the  invention; 

FIG  la  is  a  schematic  block  diagram  of  another  embodi¬ 
ment  of  the  invention; 

FIG.  2  is  a  schematic  illustration  of  the  interaction  of 
10  incoming  optical  signals  with  the  interference  filter  of  the 
KEG.  1  system; 

FIG.  3  is  a  schematic  illustration  of  the  interaction  of 
incoming  optical  signals  with  the  interference  filter  of  FIG. 
2  tilted  with  respect  to  the  optical  axis  of  the  FIG.  1  system; 
15  FIG.  4  is  a  schematic  illustration  of  the  interaction  of  an 
off-axis  incoming  optical  signal  with  the  interference  filter 
of  the  FIG.  1  system;  and 

FIG.  5  is  a  schematic  illustration  of  the  interaction  of  an 
off-axis  incoming  optical  signal  with  the  interference  filter 
20  tilted  with  respect  to  the  optical  axis  of  the  FIG  1  system  and 
with  the  addition  of  a  minor  for  providing  multiple  images 
of  the  incoming  signal. 

DETAILED  DESCRIPTION 

25 

Referring  now  to  the  drawings.  FIG  1  illustrates  in 
schematic  block  diagram  form  the  essential  components  of 
a  representative  spectrometer  system  10  according  to  the 
invention.  In  system  10.  any  suitable  optical  train  11 
3Q  (minors,  etc)  may  be  disposed  for  directing  an  incoming 
optical  signal  12  along  optical  axis  O  onto  collimating  lens 
13.  Collimated  beam  C  is  then  directed  onto  interference 
filter  15  having  transmittance  in  a  preselected  passband. 
Hlter  15  have  any  desired  passband  in  any  selected  spectral 
35  region,  ultraviolet  to  infrared,  depending  on  the  intended  use 
of  the  invention,  the  particular  selected  passband  not  con¬ 
sidered  limiting  of  the  invention.  For  example,  in  selecting 
a  passband  for  an  incoming  signal  in  the  range  of  about  3 
microns  (jim)  to  5  pm.  it  may  be  appropriate  to  select  a  filter 
40  having  four  filter  bands,  viz..  3-3.4  pm.  3.5-3.9  pm. 
4.04-4.5  pm  and  4.6-5.0  pm.  Spectral  energy  within  the 
transmission  passband  of  filter  15  is  transmitted  as  trans¬ 
mitted  beam  T  onto  imaging  lens  17  and  imaged  onto  optical 
detector  18  as  suggested  in  the  schematic  of  FIG  1.  Optical 
45  energy  within  beam  C  which  is  not  within  the  passband  of 
filter  15  is  substantially  totally  reflected  from  filter  15  at 
incident  surface  14  thereof. 

Any  suitable  optical  detector  18  sensitive  to  die  selected 
passband  of  filter  15  may  be  used,  within  scope  of  the 
50  invention,  as  would  occur  to  the  skilled  artisan  practicing  the 
invention  and  guided  by  these  teachings,  such  as  a  charge- 
couple  device,  photodiode  array  or  infrared  array,  detector 
selection  therefore  not  being  considered  limiting  of  the 
invention. 

55  Referring  now  to  FIG.  leu  shown  therein  is  a  schematic 
block  diagram  of  another  embodiment  of  the  invention  in  the 
form  of  representative  system  20.  wherein  incoming  signal 
12  is  formed  into  multiple  images  separated  spatially  and 
spectrally.  In  accordance  with  a  principle  of  operation  for  the 
60  FIG.  la  embodiment  filter  15  is  disposed  at  preselected 
angle  a  because  the  transmission  characteristics  of  interfer¬ 
ence  filter  15  having  a  specific  passband  is  a  function  of  the 
angle  of  incidence  of  the  incident  optical  energy.  Selection 
of  the  angle  a  of  filter  15  with  respect  to  axis  O  may  be  made 
65  by  the  skilled  artisan  practicing  the  invention  and  will  the 
material  used  and  the  selected  band  pass  for  each  of  the 
subsequent  filters.  Accordingly,  one  or  more  mirrors  21  may 
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be  disposed  between  collimating  lens  13  and  filter  15  at 
respective  separation  angles  p  and  at  different  planes  cor¬ 
responding  to  different  azimuthal  angles  of  incidence  of 
signal  12  relative  to  optical  axis  O  (the  separation  angle  p  of 
one  such  mirror  21  is  illustrated  in  FIG.  la)  to  reflect  the 
beam  R  reflected  from  surface  14  onto  filter  15  at  an  angle 
of  incidence  different  from  the  angle  of  incidence  of  colli¬ 
mated  beam  C.  The  transmission  character  of  filter  15  is 
dependent  on  the  cosine  of  the  angle  of  incidence  of  the 
incident  beam,  and  filter  15  therefore  transmits  differently 
for  the  reflected  beam  R  as  redirected  by  mirror  21  than  for 
collimated  beam  C.  The  selected  angle  p  between  mirror  21 
and  filter  15  will  correspond  to  a  particular  field  angle  for 
incoming  signal  12.  and  to  its  unique  cosine  dependent  filter 
value  (passband).  Multiple  (depending  on  the  number  of 
mirrors  21  included  in  system  20)  contiguous  images  are 
therefore  formed  at  detector  18  from  a  single  field  of  view 
or  incident  signal  12.  and  each  such  image  has  its  own 
unique  spectral  content. 

FIG.  2  shows  a  simplified  schematic  illustration  of  the 
interaction  of  incoming  optical  signal  12  with  filter  15  in 
system  10  of  FIG.  1.  Consider  a  filter  15  having  a  bandpass 
in  the  red  end  of  the  optical  spectrum.  For  an  optical  signal 
12  comprising  substantially  white  light,  that  portion  25  of 
signal  12  corresponding  to  a  field  angle  along  optical  axis  O 
would  result  in  a  transmitted  beam  26  in  the  red  correspond¬ 
ing  to  the  filter  15  passband.  For  any  portion  28  of  signal  12 
corresponding  to  a  field  of  view  at  successively  greater 
incidence  angles  on  filter  15.  the  corresponding  transmitted 
beam  29  would  be  shifter  to  a  shorter  (green)  wavelength  as  30 
a  function  of  the  cosine  of  the  angle  of  incidence.  For 
example,  filter  15  has  a  bandpass  at  normal  incidence  of  one 
micron,  then  for  a  portion  28  in  the  field  of  view  incident  at 
60°.  the  bandpass  would  be  about  one-half  micron  and 
somewhat  more  narrow  than  the  passband  at  normal  inci-  35 
dence.  Referring  now  to  FIG.  3.  it  is  noted  that  if  filter  15  is 
inclined  with  respect  to  optical  axis  O.  as.  for  example, 
shown  in  FIG  la.  portion  25  of  signal  12  incident  along  axis 
O  has  a  corresponding  transmitted  beam  31  shifted  to  the 
green,  and  portion  28  has  a  corresponding  transmitted  beam  40 
32  shifted  further  into  the  blue.  Accordingly,  and  with 
reference  now  to  FIG.  4.  wherein  filter  15  is  disposed  for 
normal  Incidence  of  signal  42  along  optical  axis  O.  object  41 
comprising  a  white  source  would  result  in  an  image  44  on 
detector  18  having  an  end  45  nearest  axis  O  shifted  in  45 
spectral  content  to  the  yellow  and  an  end  46  furthest  from 
axis  O  shifted  in  spectral  content  further  into  the  blue. 

Referring  now  to  FIG.  5.  shown  therein  is  a  schematic 
illustration  of  the  interaction  of  a  white  signal  52  from  object 
51  with  inclined  (with  respect  to  axis  O)  filter  15  and  mirror  50 
21  as  in  the  FIG  la  system.  That  portion  of  signal  52  which 
passes  through  filter  15  without  reflection  is  indicated  by 
solid  lines  and  results  in  a  first  image  54  having,  similarly  to 
the  illustration  of  FIG.  4.  spectral  content  shifted  toward  the 
yellow  nearest  axis  O  and  further  into  the  blue  furthest  from  55 
axis  O.  In  addition,  reflected  portion  R  (shown  by  broken 
lines)  of  signal  52  is  redirected  toward  filter  15  by  reflection 
off  mirror  21.  where  the,  total  deviation  is  now  twice  the 
separation  angle  between  filter  15  and  mirror  21.  which 
results  in  a  second  transmitted  beam  57  (also  shown  by 
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broken  lines)  and  a  second  image  58  shifted  in  spectral 
content  further  into  the  red  with  distance  from  axis  O.  Each 
additional  mirror  21  (not  shown)  may  be  disposed  at  a 
respective  preselected  angles  with  respect  to  filter  15  to 
5  provide  an  additional  image  contiguous  with  images  derived 
from  reflections  from  the  other  mirrors,  each  such  image 
containing  its  unique  spectral  content.  Contiguous  images 
will  result  where  toe  separation  angle  between  filter  15  and 
a  respective  mirror  21  is  one-half  toe  maximum  field  angle 
10  subtended  by  the  object 

The  invention  therefore  provides  an  interfering  imaging 
spectrometer  wherein  a  single  interference  filter  and  detector 
separates  an  image  into  its  spectral  components,  and  a  single 
set  of  optics  collects  all  optical  energy  needed  to  form  the 
is  desired  images.  The  invention  has  no  moving  parts,  is  simple 
and  accurate  in  its  calibration,  and  is  highly  light  efficient 
which  is  important  for  high  speed  applications.  It  is  under¬ 
stood  that  modifications  to  the  invention  may  be  made  as 
might  occur  to  one  with  skill  in  the  field  of  the  invention 
20  within  the  scope  of  the  teachings  herein  and  of  toe  appended 
claims.  All  embodiments  contemplated  hereunder  which 
achieve  toe  objects  of  the  invention  have  therefore  not  been 
shown  in  complete  detail.  Other  embodiments  may  be 
developed  without  departing  from  the  spirit  of  the  invention 
25  or  from  the  scope  of  the  appended  claims. 

I  claim; 

1.  An  interfering  imaging  spectrometer  system,  compris¬ 
ing: 

(a)  optical  means  for  directing  an  optical  signal  along  a 
preselected  optical  axis; 

(b)  a  collimating  lens  disposed  along  said  optical  axis  for 
collimating  said  optical  signal  and  directing  a  colli¬ 
mated  optical  signal  along  said  optical  axis; 

(c)  an  interference  filter  disposed  along  said  optical  axis 
for  selectively  filtering  said  collimated  optical  signal, 
said  interference  filter  having  transmittance  in  a  prese¬ 
lected  passband  and  disposed  at  a  preselected  angle  to 
said  optical  axis  for  transmitting  a  first  transmitted 
signal  along  said  optical  axis  and  reflecting  a  reflected 
signal  at  an  angle  corresponding  to  said  preselected 
angle  of  said  interference  filter; 

(d)  a  mirror  disposed  near  said  interference  filter  in 
preselected  angular  relationship  to  said  interference 
filter  and  said  optical  axis  for  directing  said  reflected 
signal  onto  said  interference  filter  and  providing  a 
second  transmitted  signal  along  said  optical  axis  and  a 
second  reflected  signal  from  said  interference  filter, 
said  second  transmitted  signal  having  spectral  content 
dependent  on  toe  angular  relationship  between  said 
interference  filter  and  said  mirror  and  defining  an  image 
separated  spatially  and  spectrally  from  toe  image 
defined  by  said  first  transmitted  signal;  and 

(e)  an  imaging  lens  and  an  optical  detector  disposed  along 
said  optical  axis,  said  imaging  lens  disposed  for  imag¬ 
ing  the  transmitted  signals  onto  said  optical  detector. 

2.  The  spectrometer  system  of  claim  1  wherein  said 
optical  detector  is  a  charge-couple  device,  photodiode  array 
or  infrared  array. 
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